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Abstract 



This is an introduction to geometrie algebra, an alternative to traditional vector algebra that expands 
on it in two ways: 

1. In addition to scalars and vectors, it defines new objects representing subspaces of any dimension. 

2. It defines a product that's strongly motivated by geometry and can be taken between any two 
objects. For example, the product of two vectors taken in a certain way represents their common 
plane. 

This system was invented by William Clifford and is more commonly known as Clifford algebra. It's 
actually older than the vector algebra that we use today (due to Gibbs) and includes it as a subset. Over 
the years, various parts of Clifford algebra have been reinvented independently by many people who 
found they needed it, often not realizing that all those parts belonged in one system. This suggests that 
Clifford had the right idea, and that geometrie algebra, not the reduced version we use today, deserves 
to be the Standard "vector algebra." My goal in these notes is to describe geometrie algebra from that 
standpoint and illustrate its usefulness. The notes are work in progress; Pil keep adding new topics as I 
learn them myself. 
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1. Introduction 
1.1. Motivation 

I'd say the best intuitive definition of a vector is "anything that can be represented by arrows that add 
head-to-tail." Such objects have magnitude (how long is the arrow) and direction (which way does it point). 
Real numbers have two analogous properties: a magnitude (absolute value) and a sign (plus or minus). 
Higher-dimensional objects in real vector spaces also have these properties: for example, a surface element 
is a plane with a magnitude (area) and an orientation (clockwise or counterclockwise). If we associate real 
scalars with zero-dimensional spaces, then we can say that scalars, vectors, planes, etc. have three features 
in common: 

1. An attitude: exactly which subspace is represented. 

2. A weight: an amount, or a length, area, volume, etc. 

3. An orientation: positive or negative, forward or backward, clockwise or counterclockwise. No matter 
what the dimcnsion of the space, there are always only two orientations. 

If spaces of any dimension have these features, and we have algebraic objects representing the zero- and 
one-dimensional cases, then maybe we could make objects representing the other cases too. This is exactly 
what geometrie algebra gives us; in fact, it goes farther by including all of these objects on equal footing in 
a single system, in which anything can be added to or multiplied by anything else. I'll illustrate by starting 
in three-dimensional Euclidean space. 

My goal is to create a product of vectors, called the geometrie product, which will allow me to build up 
objects that represent all the higher-dimensional subspaces. Given two vectors u and v, traditional vector 
algebra Iets us perform two operations on them: the dot product (or inner product) and the cross product. 
The dot product is used to project one vector along another; the projection of v along u is 

Pu{v) = 1 

\u\ 2 

where u ■ v is the inner product and \u\ 2 = u ■ u is the square of the length of u. The cross product represents 
the oriented plane defined by u and v; it points along the normal to the plane and its direction indicates 
orientation. This has two limitations: 

1. It works only in three dimensions, because only there does every plane have a unique normal. 

2. Even where it works, it depends on an arbitrarily chosen convention: whether to use the right or left 
hand to convert orientations to directions. So the resulting vector does not simply represent the plane 
itself. 

Because of this, I'll replace the cross product with a new object that represents the plane directly, and it 
will generalize beyond three dimensions as easily as vectors themselves do. 

I begin with a formal product of vectors uv that obeys the usual rules for multiplication; for example, 
it's associative and distributive over addition. Given these rules I can write 

uv = — (uv + vu) + — (uv — vu). (2) 
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The first term is symmetrie and bilinear, just like a generic inner product; therefore I set it equal to the 
Euclidean inner product, or 

— (uv + vu):=u-v. (3) 

I can immediately do something interesting with this: notice that u 2 = u-u = |u| 2 , so the square of any 
vector is just its squared length. Therefore, the vector 

u- 1 := — 2 (4) 

is the multiplicative inverse of u, since obviously ra -1 — u 2 /u 2 = 1. So in a certain sense we can divide by 
vectors. That's neat. By the way, the projection of v along u from Eq. ([T]) can now be written 

P u (v) = (vu)u- 1 . (5) 

In non-Euclidean spaces, some vectors are null, so they aren't invertible. That means that this projection 
operator won't be defined. As it turns out, projection along noninvertible vectors doesn't make sense 
geometrically; IT1 explain why in Section 17.11 Thus we come for the first time to a consistent theme in 
geometrie algebra: algebraic properties of objects frequently have direct geometrie meaning. 

What about the second term in Eq. @? I call it the outer product or wedge product and represent it 
with the symbol A, so now the geometrie product can be decomposed as 

uv = u-v + uAv. (6) 

To get some idea of what u Av is, I'H use the fact that it's antisymmetric in u and v, while u ■ v is symmetrie, 
to modify Eq. © and get 

vu ~ u ■ v — uAv. (7) 

Multiplying these equations together I flnd 

uvvu = (u ■ v) 2 — (u A v) 2 . (8) 
Now vv = \v\ 2 , and the same is true for u, while u ■ v = \u\\v\ cos ö, so 

{uAv) 2 = -\u\ 2 \v\ 2 sm 2 9. (9) 
So whatever tiA» is, its square has two properties: 

1. It's a negative scalar. (Just like an imaginary number, without my having to introducé them separately. 
Hmm.) 

2. Aside from the minus sign, it is the square of the magnitude of the cross product. 

The first property means that u A v is neither scalar nor vector, while the second property makes it look like 
a good candidate for the plane spanned by the vectors. u Av will turn out to be something called a simple 
bivector or 2-blade, so 2-blades represent planes with an area and an orientation (interchange u and v and 
you change the sign oï uAv). There's no unique parallelogram associated with uAv because for any A, 

u A(v + Am) = uAv. (10) 

So sliding the tip of one side along the direction of the other side changes the parallelogram but not the outer 
product. It is the plane (attitude), area (weight), and orientation that the outer product defines uniquely. 
With these definitions, the product of two vectors turns out to be the sum of two very different objects: a 
scalar and a bivector. For the moment think of such a sum as purely formal, like the sum of a real and an 
imaginary number. 

Later Hl define the outer product of any number of vectors, and this product will be associative: 

(u A v) A w = u A(v A w) = u A v A w. (11) 
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This guy is called a simple trivector or 3-blade, and it represents the three-dimensional space spanned by 
its factors, again with a weight (volume) and orientation. We can also form 4-blades, 5-blades, and so on 
up to the dimension of whatever vector space we're in. Each of these represents a subspace with the three 
attributes of attitude, weight, and orientation. These r-blades and their sums, called multivectors, make up 
the entire geometrie algebra. (Even scalars are included as 0-vectors.) The geometrie product of vectors can 
be extended to the whole algebra; you can multiply any two objects together, which Iets you do all sorts of 
useful things. Just multiplying vectors already Iets us do a lot, as I'll show now. 



1.2. Simple applications 

I'll start by solving two Standard linear algebra problems. Let's suppose a plane is spanned by vectors a and 
b, and you have a known vector x in the plane that you want to expand in terms of a and b. Therefore you 
want scalars a and j3 such that 

x = aa + /3b. (12) 
To solve this, take the outer product of both sides with a; since aAo = 0, you get 

aAx = l3aAb. (13) 

It will turn out in Euclidean space that every nonzero vector, 2-blade, and so on is invertible, so this can be 
solved to get 

/3 = (oAijfaA^ 1 . (14) 

This makes sense geometrically: both a Ax and a A b are bivectors in the same plane, so one should be a 
scalar multiple of the other. Since (3 is effectively a ratio of areas, I'm going to write instead 

aAb 

The problem with this is that it could mean either (aAi)(flAi) _1 or (a A &)~ 1 (a A x); but in this case 

they're the same, so there's no harm. Taking the outer product of both sides with b similarly gets you 
a = (x A b)/(a A &), so now we know that 

xAb\ /aAx\ , ,„„. 

r a + ] b. (16 

aAbJ KaAbJ y ' 

This expression is called Cramer's Rule. Here I've derived it much more quickly than is done in regular vector 
algebra, it's expressed directly in terms of the vectors instead of in components, and the geometrie meaning 
of the cocfBcicnts (ratios of areas in the plane) is immcdiatcly apparent. Also note that this expression is 
dcfined iff a A b ^ 0, which is exactly the condition that a and b span the plane. 

The generalization from planes to volumes is straightforward; if a, b, and c span the space then 

xAbAc\ /sAiAc\ , faAbAx 



AbAc \a Ab A 



\ , / aAóAx \ 

) i+ U^ c (l7) 



and so on for higher dimensions. 

When you have a linear equation like this, taking the outer product with one of the terms, and thus 
removing that term, is often a handy trick. Here's another example. Suppose I have two lines that lie in a 
plane: The fust passes through point p and points in direction a, while the second passes through point q 
and points in direction b. Assuming the lines aren't parallel, at what point x do they cross? 

If the lines aren't parallel then a and b aren't parallel, so they span the plane. Therefore x is a linear 
combination of a and b as given by Eq. ()16[) . That's nice but unhelpful, because this time x is unknown and 
we're trying to solve for it. But wait; x lies on the line through p pointing along a, or 

x = p + Xa (18) 
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for some A. That means that aAx = aAp. And the fact that x lies on the line through q pointing along b 
tells me that x A b — q A b, so when I put all this in Eq. (I16|) I find that the intersection point x is 




expressing the unknown x in terms of the four known vectors defining the two lines. 

The solutions to these last two exercises are expressed in an entirely intrinsic, coordinate-free way, which 
means that the results of this calculation can be used as inputs in any further calculations. Once you get 
to the end, of course, you can certainly use coordinates to perform the final computations. To do all this, 
though, you have to be comfortable with these new kinds of products and their inverses. Fm here to help 
with that. 

Now for a little geometry. I'U start by looking at reflections, like the operation performed by a mirror. 
How do we perform a mirror refiection on a vector? Well, we often think of a reflection as happening in one 
of two complementary ways: either through a plane (components in the plane are left alone, the remaining 
component gets a minus sign) or along an axis (the component along the axis gets a minus sign and the othcr 
components are left alone). However, these ways of thinking are interchangeable only in three dimensions, 
because only there does any plane have a unique normal. I want a picture that works in any number of 
dimensions, and only the second one does that, because it works even in one dimension. So I'll use it from 
now on. 

Let v be the vector we want to rcflect and let n be a vector along the reflection axis. Then 

v = ti(nn _1 ) 

= (mi)h _1 + (uAii)ii _1 . (20) 

The first term looks like the right hand side of Eq. ([5]), so it represents the orthogonal projection of v along 
n. That means the other term is the component of v perpendicular to n, also called the orthogonal rejection 
of v from n. (I'll bet you've never heard that term before.) Now let v' be the refiected vector; its component 
along n has the opposite sign, while its perpendicular component is the same, so it is given by 

v' = —(v ■ + (v A n)n~ x . (21) 

Using the symmetry and antisymmetry of the inner and outer products respectively, I can recast this as 

V 1 = —(n-v)n^ 1 — (nAv)n^ 1 
= — (n ■ v + n A v)n~ x 

= -nvn~ x . (22) 

This is a nifty little result; the equation for rcflccting a vector along an axis is very tidy. Compare that to 

n ' v 

v' = v-2—n, (23) 
\n\ 

which is the simplest one can do with traditional vector algebra. 

The appearance of both n and n^ 1 in Eq. (f22)) guarantees that the result depends neither on the weight 
(length) nor the orientation of n, only its attitude (the axis it represents), as it should. 

The next operation I'll describe is rotation. First note that the usual way one thinks of rotations, 
as being performed around an axis, works only in three dimensions. In general, it is better to think of 
a rotation of a vector as being performed in a plane; the component in the plane is rotated while the 
components perpendicular to the plane are left alone. Again, this picture works perfectly well in any number 
of dimensions. 

Hamilton discovered a great way to perform rotations: To rotate through angle 9 in a plane, perform two 
reflections in succession along any two axes in the plane, as long as (a) the angle between the axes is Ö/2 
and (b) a rotation from the first axis to the second is in the same direction as the rotation to be performed. 
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Figure 1: The vector v makes an angle <f> with axis n. It is then reflected along n, producing vector v', then 
along axis m, producing vector v". Notice that the angle between vectors v and v" is 6, twice the angle 
between n and m, regardless of the value of <f>. 



This is shown in Figure [TJ So if I want to rotate vector v, then I let m and n be vectors along axes satisfying 
the conditions, and the result of the rotation is 

v' = —m(—nvn~ 1 )m~ 1 
— (ran)i)(ran) _1 

= RvR- 1 (24) 

where R = mn. R is an object called a rotor. Typically a rotor is the product of unit vectors, in which case 
m = Tïh~ and n — n" 1 , which means R^ 1 = nm. 

These two examples, reflections and rotations, introducé a second theme in geometrie algebra: elemcnts 
of the algebra represent both geometrie objects (vectors, subspaces) and operations on those objects. There's 
no need to introducé any new elements (e.g. matrices) to represent operators. 

Til finish this section by looking at the rotor more closely. 

R = mn 

— m ■ n + m An 

= m-n — nAm. (25) 

I reversed the order of m and n because the sense of rotation of this rotor is from n to m (reflection n 
was applicd first). Since m and n are now unit vectors and the angle between the two axes of reflection is 
0/2, m-n = cos(ö/2) and (nAm) 2 = — sin 2 (ö/2). Therefore the bivector B — (nAm)/ sin(ö/2) is a unit 
bivector: B 2 = — 1. So now 

R = cos(0/2) - sin(Ö/2)B 

= exp(-B6/2) (26) 
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where the exponential is defined by its power series; the scalar terms and the terms proportional to B can 
be grouped and summed separately. Now I have a rotation operator that explicitly displays the rotation's 
angle, plane (attitude of B), and direction (orientation of B), and all without coordinates. 

Recall how little I started with: a product of vectors with the minimal algebraic properties to be useful, 
plus the extra bit that the symmetrie part equals the inner product. From only that much, Fve gotten a 
formula for rotating a vector that looks a lot like the formula for rotating a complex number, z' = e lS z, 
except that it's double-sided and uses half of the rotation angle. The resemblance to complex numbers is 
no accident; as I will show later on, the complex numbers are contained in the geometrie algebra of the 
real Euclidean plane. Therefore, all of complex algebra and analysis is subsumed into and generalized to 
arbitrary dimensions by geometrie algebra. As for the half angles, in physics they normally show up in the 
quantum theory of half-integer spin particles, but this suggests that there's nothing particularly quantum 
about them; they arise simply because a rotation equals two reflections. 

1.3. Where now? 

When I first read about geometrie algebra, examples like these immediately made me think it might have 
a lot to offer in terms of conceptual simplicity, unifying potential, and computational power. I was on the 
lookout for something like this because I had found Standard mathematical physics dissatisfying in two main 
ways: 

1. We use a hodgepodge of different techniques, particularly in more advanced work, each of which 
seems to find a home only in one or two specialized branches of theory. It seems like an unnecessary 
fragmentation of what should be a more unihed subject. 

2. As long as you stay in three dimensions and work only with vectors, everything is very concrete, 
geometrical, and easy to express in intrinsic form without coordinates. However, none of these desirable 
features seem to survive in more general situations (say, special relativity) . Either you end up expressing 
everything in coordinates from the start, as in classical tensor analysis, or you use coordinate-free forms 
like those found in modern differential geometry, which I hnd hard to calculate with and which seem 
to leave their geometrical roots in favor of some abstract analytic form. (I put differential forms in this 
category.) 

Despite signs of promise, however, I also have to admit I was taken aback by what looked like an enormous 
prolifcration of new objects. After all, it seems like geometrie algebra Iets you multiply vectors all day long 
and keep getting new things, and I had very little sense of how they all related to each other. (I imagine this 
is why geometrie algebra lost out to Gibbs' vector algebra in the first place.) I was also puzzled about how 
the rules for using these objects really worked. For example, if I had just read the previous two sections, Fd 
have questions like these. 

1. In my first two examples in Section [F2l I used the inverse of a 2-blade, (aAb)^ 1 , and I mentioned 
that in Euclidean space every nonzero r-blade has an inverse. Fve shown how to calculate the 2-blade 
itself: it's the antisymmetrized product. But how do you calculate the inverse? 

2. In Eq. (|22|) . I multiplied three vectors, nvn~ x , and the result was also a vector. However, you could 
teil that only by following the derivation. What's the product of three vectors in general? Is it always 
a vector? Is it something else? How can you teil? 

3. Then I multiplied a bivector by a vector, (u A What's that? In this case the result was anothcr 
vector, but again you had to follow the derivation to know that. In addition, I also said it was 
perpendicular to n. How do I check that? Presumably I should show that 

[{vAn)n- 1 ]-n = 7 (27) 

but that looks scary. 

To answer these and other questions for myself, I wrote these notes. I suspect Fm not the only one who 
reacted this way on seeing geometrie algebra for the first time, so I hope the notes can help others understand 
geometrie algebra and decide for themselves whether it's as good as advertised. 

The structure of the notes reflects the best way Fve found to explain geometrie algebra to myself. 
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• In Section [3] I lay out a set of axioms. I find it helpful to present axioms first, so we can have all of our 
basic rules in place immediately, knowing we won't be surprised later by having to take anything else 
into account. 

• With the axioms in hand, in Section [3] I answer the first question I asked myself when I saw all this: 
exactly what's in here? I describe both algebraically and geometrically what a generic multivector 
looks like, and I justify the claims I made at the end of Section [TTT1 in terms of the axioms. By this 
point, a multivector should seem a lot more concrete than just "lots of vectors multiplicd togethcr." 

• Having explained what a general multivector looks like, in Section|3]I explain what a general product of 
multivectors looks like. I also explain how to take the inner and outer products of any two multivectors, 
and I explain what they mean geometrically; this is a natural continuation of the geometrical discussion 
in Section[3J I claimed earlier that geometrie algebra Iets you take coordinate-free, intrinsic calculations 
much farther than Standard methods; it does this because it has a large number of algebraic identities, 
which I'll start to derive here. These identities make vector algebra start to look a lot more like algebra 
with real numbers. 

• A handful of additional operations are used all the time in calculations, and I collect them in Section 
[5j By describing them all together, I can show the relationships between them more easily. 

• At this point even I think the reader needs a break, so I pause in Section [üj for an "application" by 
describing what our favorite vector spaces, two- and three-dimensional real Euclidean space, look like 
in these terms. I show how the complex numbers pop up all by themselves in the two-dimensional 
algebra, and in three dimensions I show how to convert back and forth from geometrie algebra to the 
traditional language of cross products, triple products, and so on. 

• With the full algebra at my disposal, in Section [7] I return with a vengeance to my initial examples: 
orthogonal projection, reflections, and rotations. Now I really get to show you why this isn't your 
grandpa's vector algebra. We can project vectors into subspaces, and even subspaces into other sub- 
spaces, far more easily than traditional methods ever made you think was possible. And wait till you 
see what rotations look like. Ever tried to rotate a plane? In geometrie algebra, it's easy. 

• Coordinates do have a role in geometrie algebra, although it's vastly reduced, and I describe it in 
Section [3 

• Linear algebra looks very different when it's done not just on vector spaces but on geometrie algebras; 
that's the subject of Section [9] I'll review the basics, but even familiar subjects like adjoints and skew 
symmetrie operators take on a new flavor and signihcance in this system. And eigenvectors will be 
joined by eigenspaces of any dimension. I'll even show how to act with a linear operator on the whole 
vector space at once, and the eigenvalue of that operation will be our friend the determinant. 

• Right now, the notes are very light on applications to physics; so far I have included only a brief 
discussion of classical angular momentum (which is no longer a vector, by the way) and the Kepler 
problem, which gets a pretty snazzy treatment. I'll add more applications soon. 

All the important dcfinitions and relations are listed together in Appendix [XJ and the topics I plan to include 
in future versions are listed in Appendix [Bj 

I. 4. References and comments 

Although geometrie algebra dates from the 19th century, it was recovered in the form described here only in 
the 20th century by David Hestenes |1I2I3| . and it is slowly gaining popularity in various math and applied 
math communities. My primary sources are Hestenes' books; Doran and Lasenby [1], written for physicists; 
Dorst, Fontijne, and Mann 5 , written for computer scientists; and the introductory linear algebra text by 
Macdonald [B], which includes geometrie algebra alongside traditional linear algebra topics. You'11 see their 
influence everywhere; for example, my axioms were inspired by [2], Section fl.ll and the second half of Section 

II. 21 come from [4], and the first half of Section Tl. 21 is lifted from [5]. I'll mention other areas where I'm 
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particularly indebted my sources as I come to them. I follow [5] in defining two inner products, instead of 
Hestenes' one, but I continue to refer to them as inner products instead of "contractions" as Dorst et al. 
do. Finally, this approach to geometrie algebra is far from the only one: Lounesto [7] describes this one and 
several others, and he gives a great overview of the history that has brought Clifford algebra to this point. 

Given all the other introductions to geometrie algebra out there, I hope this treatment is made distinctive 
by two elements. First, I have worked everything out in more detail than I've seen anywhere else, which 
I think is very helpful for getting one's initial bearings in the subject. Second, I don't believe this way of 
organizing the material is found in other sources either, and as I said in the previous section, this is the 
way I've found easiest to understand. I try to convey Hestenes' attitude toward Clifford algebra as not just 
anothcr algebraic system but the natural extension of real numbers to include the geometrie idea of direction, 
which I find very attractive. 

I also prefer a more general treatment over a more specific one when the theory seems to be equally easy 
in either case. For example, all applications of geometrie algebra Fm familiar with take the scalars to be M, 
the real numbers, and an important part of Hestenes' view is that many of the other number systems used in 
mathematics are best understood not separately but as subsets of certain real Clifford algebras. (I dropped 
a hint about this regarding complex numbers in Section 11.11 to which I'll return in Section [6j where I'll 
handle the quaternions too.) However, I don't force the scalars to be R here, because the majority of results 
don't actually depend on what the scalars are. One thing that does change a bit, however, is the geometrical 
interpretation. For example, suppose the scalars are complex; how does the orientation of a vector change 
when you multiply by il In fact, the two notions of weight and orientation make sense only for real vector 
spaces, and as a result they won't have a place in a general geometrie algebra. They're still important for all 
those applications, however, so I'll make sure to explain them at the right time. And whenever the scalars 
have to be real for something to be true, I'll say so. 

As part of my goal to work everything out in detail but keep the notes easy to follow, I've set the theorem 
proofs off from the rest of the text so they can be easily skipped. Nonetheless, I urge you to take a look at the 
shorter proofs; I tried to motivate them well and convey some useful insights. Even some of the long proofs 
consist of more than just turning the algebra crank. I like proofs that do more than show that something is 
true; they give a sense of why. I have tried to write those sorts of proofs here. 

Because geometrie algebra has found its way into most applied mathematics, albeit in a very fragmented 
way, everything I describe in these notes can be done using some other system: matrices, Grassmann 
algebras, complex numbers, and so on. The advantage that I see here is that one system, a natural extension 
of elementary vector algebra, can do all these things, and so far I've always found I can better understand 
what's going on when all the different results are relatcd through a unified perspective. 

2. Definitions and axioms 

The purpose of this section is to define a geometrie algebra completely and unambiguously. This is the 
rigorous version of the discussion from Section 11.11 and you'11 see all of the basic ideas from that section 
reintroduced more precisely here. 

A geometrie algebra is a set Q with two composition laws, addition and multiplication (also called the 
geometrie product), that obey these axioms. 

Axiom 1. Q is a ring with unit. The additive identity is called and the multiplicative identity is called 1. 

Axiom [T] is the short way to say that (a) addition and multiplication in Q are both associative, (b) 
both operations have identities, (c) every element has an additive inverse, (d) addition commutes, and (e) 
multiplication is left and right distributive over addition. So now I've said it the long way too. 

A generic element of Q is denoted by a capital Roman letter (A, B, etc.) and is called a multivector. 
Notice that a geometrie algebra is one big system from the get-go: all multivectors, which will eventually 
include scalars, vectors, and much more, are part of the same set, and addition and multiplication are equally 
available to all. I'll continue to follow this philosophy as I introducé new operations by defining them for all 
multivectors. Also, this axiom formalizes the first requirement I made of the geometrie product in Section 
11.11 it gives addition and multiplication the minimal properties needed to be useful. 

Axiom 2. Q contains a field Q$ of characteristic zero which includes and 1. 
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A member of Go is called a O-vector, a homogeneous multivector of grade 0, or a scalar. Scalars are 
denoted by lower case Greek letters (A, fi, etc.). Being a field means that Go is closed under addition and 
multiplication, it contains all inverses of its elements, and it obeys all the rules that Q obeys from Axiom 
Q] plus the additional rules that (a) everything except has a multiplicative inverse and (b) multiplication 
commutes. The rational numbers, real numbers, and complex numbers are all fields. The property of having 
characteristic zero saves me from getting in trouble in the following way Since Go doesn't have to be R, 
the integers aren't actually the usual integers, but sums of terms all equaling 1, the multiplicative identity 
of Go- (So, for example, by 2 I literally mean 1 + 1.) If I don't specify any further properties of Go, then I 
haven't ruled out 1 + 1 = 0, which would be bad when I try to divide by 2. (Which I'll be doing frequently; 
see Eq. ([2]).) Having characteristic zero means that no finite sum of terms all equaling 1 will ever add up to 
0, so I can divide by integers to my heart's content. 

Axiom 3. G contains a subset Gi closed under addition, and X € Go, v € Gi implies Xv = vX G Gi- 

A member of Gi is called a 1-vector, a homogeneous multivector of grade 1, or just a vector. Vectors are 
denoted by lower case Roman letters (a, b, u, v, etc). The axioms imply that Gi obeys all the rules of a 
vector space with Go as scalars, justifying their names. However, all is not the same as what you're used 
to. In Standard vector algebra, the scalars and vectors are usually separate sets. For example, consider the 
vector space R 3 with the real numbers as scalars; the zero scalar is the number 0, but the zero vector is the 
ordered triple (0,0,0). In geometrie algebra this is not the case, and here's why 

1. is a scalar by Axiom [2] 

2. 0i> = for any vector v by Axiom [T] 

3. A scalar times a vector is a vector by Axiom [3] 

4. Therefore, is also a vector. 

It will turn out that is a whole lot of other things too. 

So far the axioms have told us how to add scalars, add vectors, multiply scalars, and multiply a scalar 
and a vector. Multiplying vectors is next. 

Axiom 4. The square of every vector is a scalar. 

As it was in Section fl.il this is the most important axiom of the bunch. Here's the first consequence: for 
any vectors u and v, 

l i 

- {uv + vu) = - [{u + vf - u 2 - v 2 } (28) 

(you can easily prove this by expanding out the right hand side), and the right side is a scalar thanks to 
Axiom[4l so it follows that the symmetrized product of any two vectors is a scalar. In fact, this is not merely 
implied by Axiom |4j it's equivalent to it. (Assume the statement is true. Since the square of a vector is its 
symmetrized product with itself, Axiom 2] follows.) The symmetrized product of two vectors defined above 
is called their inner product and is denoted either u\ v or u \_v. It is symmetrie and linear in both terms, 
thus obeying the usual rules for an inner product on a real vector space (but not a complex vector space). 
Vectors u and v are said to be orthogonal iï u \ v = 0, u is & unit vector if u 2 = ±1, and u is null if u = 0. 
Notice that vectors are orthogonal iff they anticommute. This turns out to be handy. Recall my earlier 
comment that if v is non- null, then v is invertible and v^ 1 = v/v 2 . 

I have two inner products, J and [, instead of just the usual •, for reasons that won't be clear until 
Section POl However, the two products are equal when both factors are vectors, so I can continue to use the 
Standard terminology of inner products as I please. The next axiom is an example. 

Axiom |4]by itself is a little too general; for instance, it would allow the product of any two vectors to be 
zero. That seems pointless. To prevent that, I'll add another axiom. 

Axiom 5. The inner product is nondegenerate. 
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This means that the only vector orthogonal to all vectors, including itself, is 0. This axiom is true in every 
application I can imagine, and I use it to prove some useful results in Section [5.5l However, it is possible to 
replace it with a weaker axiom that accomplishes most of the same things; I discuss that in Section [5.5l too. 
So if you ever find yourself reading other treatments of Clifford algebras, watch out to see whether they use 
this axiom or put something else in its place. 

Now I'll name other elements of Q. Let r > 1; then an r-blade or simple r-vector is a product of r 
orthogonal (thus anticommuting) vectors. A finite sum of r-blades is called an r-vector or homogeneous 
multivector of grade r. (I'll bet you didn't see that coming.) 2-vectors are also called bivectors, 3-vectors 
trivectors. The set of r- vectors is called Gr- Notice that this definition of simple r- vectors uses the geometrie 
product of orthogonal vectors, not the outer product of arbitrary vectors as I did in Section 11.11 The 
dehnitions are equivalent, as I'll show later. 

Products of vectors play an important role, so they get their own name. An r-versor is a product of r 
vectors. So far we've seen two types of versor: blades (where the vectors in the product are orthogonal) and 
rotors, introduced in Section [1.2l A rotor was defined to be a product of two invertible vectors, so a rotor is 
an invertible biversor. Later, a rotor will be any invertible even versor. 

From these dehnitions and Axiom [3] it follows that that if A £ Gr and A is a scalar, 

XA = AX G Gr- (29) 

So multiplication by a scalar doesn't change the grade of an r-vector. This in turn implies that (a) each Gr 
is a vector space with Go as scalars and (b) € Gr for every r. So all the results I gave right after Axiom [3] 
generalize fully. 

Now we know that G contains all the Q r , and we know a few things about how different Gr are related. 
For example, suppose u and v are orthogonal and consider the 2-blade uv. It anticommutes with both u and 
v, which means that it can't have a scalar part, because that part would have commuted with all vectors. 
In fact, for this same reason no even blade can have a scalar part; and no odd blade can either, as long as 
there's another vector orthogonal to all the factors in the blade. You can continue on this line and deduce 
a few more results, but it's not clear to me that you can use only the axioms so far to show that all the Gr 
are completely independent of each other. So I add one hnal axiom for cleaning up. 

Axiom 6. If Ga — Gi, then G = Go- Otherwise, G is the direct sum of all the Gr- 

The hrst part of the axiom covers a special case: a field by itself, without any vectors, can be a geometrie 
algebra. When there are vectors around, the axiom says that every A £ G may be expressed one and only 
one way as A — ^2 r A r where A r G Gr and all but finitely many A r vanish. Therefore, every A ^ is either 
an r-vector for only one r or is of mixed grade. 

For each r, let the grade operator () : G — > Gr project each A € G onto its unique grade-r component. 
Then 

(a) A is an r-vector iff A = {A} . 

(b) (A + B) r = (A) r + {B) r . 

(c) (XA) r = (AX) r = X (A) r . 

(d) {{A) r ) = (A} r S rs . (Thus the () r are independent projection operators.) 

(e) ^ r (A) = A for any A £ G- (Thus the () are a complete set of projection operators.) 

It will turn out to be convenient to define () even when r is negative, so let me add one final property: 

(f) (A) r = if r < for all A G Q. 

Because we take the scalar part of multivectors so often, I will let ( ) mean ( ) . 

The notation A r will usually mean that A r is an r-vector. The exception is vectors: a\, 02, etc, are 
all vectors in a single enumerated set (not objects of increasing grades). Sometimes A r will represent the 
grade-r component of multivector A, which is more properly denoted (A) r , but that notation is cumbersome 
so sometimes I drop it. You can always teil from the context. A blade is indicated by boldface; for example, 
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A r is an r-blade. The exceptions are scalars (0-blades) and vectors (1-blades). I want a special notation for 
blades because they have geometrie meaning while general r- vectors don't, as Til show in Section [3] 

Axiom [ü] tells us that the relation XA = AX, which I proved above for any homogeneous multivector A, 
is true for any A G Q, homogeneous or not. Another consequence of Axiom [ü] is that Q is the direct sum 
of subspaces Q + and Q- consisting of the even-grade and odd-grade multivectors respectively. Since many 
identities will contain factors of (— l) r , they will only depend on whether the multivectors are even or odd. 
Also, Fll show in Section |4] that the product of even multivectors is also even; this means that the even 
subspace is actually a subalgebra, which will turn out to be important. For these reasons it's good to extend 
some of my notation to cover even and odd cases; the notations A + and A_ will mean that these objects 
have only even-grade or odd-grade terms, respectively, and for any multivector A, (A), (resp. (A)_) is the 
even-grade (resp. odd-grade) part of A. 

By the way, I haven't actually proved that anything satisfying these axioms exists. That's done in [5]. 

3. The contents of a geometrie algebra 

According to Axiom |6l a geometrie algebra consists of r-blades and their sums. However, the axioms and 
my comments at the end of Section 11.11 give us two different pictures of what r-blades are. According to 
the axioms, an r-blade is a product of r orthogonal vectors; according to the end of Section ll. 11 an r-bladc 
is an outer product of r arbitrary vectors. I also said that blades represent subspaces, with weights and 
orientations when the scalars are real. I'll spend this section relating these two pictures; first Fll show the 
two definitions of r-blades are equivalent, and then I'll justify the geometrie interpretation. Then we'11 have 
a good intuitive feel for what a geometrie algebra really is: sums of subspaces, with orientations and weights 
if the algebra is real. To do this, I'll be using some concepts that I haven't fully explained yet. Everything 
left hanging here will be fixed up in the next few sections. 

First I want to show that outer products of vectors really are r-blades in the axiomatic sense of Section 
[2j To do this, I define the outer product of vectors {cij}j = i,...,r to be their fully antisymmetrized product, or 



where a is a permutation of 1 through r, sgn a is the sign of the permutation (1 for even and —1 for odd), 
and the sum is over all r! possible permutations. If r = 2 this reduces to the outer product of two vectors 
dcfincd previously. Here's the result I need. 

Theorem 1. The outer product ofr vectors is an r-blade, and every r-blade is the outer product ofr vectors. 

A corollary is that the outer product of two vectors is a 2-blade, as I said in Section 11.11 This means 
that I could have used this as the definition of an r-blade in the axioms, but the definition I did use is more 
convenient in many situations. Now, however, I'll use either definition as I need to. 

Proof. To begin, I'll show that if the dj all anticommute then the outer product reduces to the 
geometrie product, so the result is an r-blade. Let {ej }j=i,. ,., r anticommute with one another, 
and consider their outer product 




(30) 



(J 




(31) 



In each term, the e$ can be reordered so they're in ascending numerical order, and each interchangc 
of two ei introducés a minus sign. The end result is a factor of the form sgn er, which cancels the 
sgn er that's already there. The result is 



ei A €2 A • • • A e r = — 
r! 




= eie 2 • • • e r 



(32) 
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since there are r! permutations to sum over and all r! terms are the same. So when the vectors 
all anticommute, the wedges can be retained or dropped as desired. The r = 2 version of this 
result, 

ei A e 2 = e\ e 2 , (33) 

was already obvious from Eq. ([6]) since e\ ■ e% = 0, or e\ J e 2 = as I would say it now. 

Turning to the general case, I can show that this is an r-blade by examining the matrix 
M with entries Mij = ai J dj . This is a real symmetrie matrix, so it can be diagonalized by an 
orthogonal transformation, meaning that there exists an orthogonal matrix R and a set of vectors 
{&i}i=i,...,r such that 

a t = ^ Ri o e J and e » J e J = e i <%> ( 34 ) 

3 

so the ei anticommute with each other. In that case 

a\ A a,2 A • • • A a r = RuCi A RijCj A • ■ ■ A R r m e m 

i,j,...,m 

= det(iï)ei Ae 2 A---Ae r . (35) 

Now det(iï)=±l, and if it equals —1 I interchange e± and e 2 and relabel them, with the result 

ai A a 2 A • • • A a r = e\ A e 2 A • • • A e r 

= eie 2 ---e r (36) 

where the final line relies on the result from the previous paragraph. So the outer product of 
r vectors can be re-expressed as the product of r anticommuting vectors, making it an r-blade. 
Further, every r-blade is such an outer product (since for anticommuting vectors the wedges can 
be added or dropped at will), so an object is an r-blade iff it's the outer product of r vectors. □ 

Since every multivector is a unique sum of r- vectors by Axiom[6l and every r- vector is a sum of r-blades, 
I can now say that a multivector is a sum of a scalar, a vector, and a bunch of outer products. Now let's 
take the geometrie point of view. I know what scalars and vectors are geometrically, but what are the 
outer products? To answer that, I need to look at when they vanish. For example, a A a = for any a by 
antisymmetry. The more general case is given by this theorem. 

Theorem 2. The simple r-vector a\ A a 2 A • • • A a r = iff the vectors {ai\i=i t ..., r are linearly dependent. 

Proof. The outer product is clearly antisymmetric under interchange of any pair of factors, 
so it vanishes if any factor repeats. It is also linear in each of its arguments, so if one factor is 
a linear combination of the others, the outer product vanishes. So if the vectors are dependent, 
their product vanishes. The other half of the proof, that the product of independent vectors 
doesn't vanish, is given in Theorem [25j which we don't have the tools to prove yet, so I'll defer 
it until later. □ 

So an r-blade A r is nonzero exactly when its factors span an r-dimensional subspace. Thus I associate A r 
with that subspace (attitude). 

To solidify the connection between subspaces and r-blades, here's a really cool result. It uses af\A r , 
which I haven't defined yet, but for now let's just say that it equals the outer product of a and the factors 
of A r . 

Theorem 3. If A r is a nonzero r-blade with r > 1, then vector a lies in the span of the factors of A r iff 
aAA r = 0. 

Proof. a A A r = iff a and the factors of A r are linearly dependent. Now the factors of A r 
are themselves independent because A r is nonzero, so a A A r vanishes iff a is a linear combination 
of the factors of A r . □ 
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Therefore A r does indeed define a subspace: the set of all vectors a such that a A A r = 0. 

The proof of this theorem actually shows a bit more. If a A A r =t 0, then it's an r + 1-blade, and it 
represents the direct sum of A r and the one-dimensional subspace defined by a. I'll use this fact later when 
I show how to interpret outer products geometrically. 

Theorem [3] implies another useful fact. 

Theorem 4. Two nonzero r-blades A r and B r define the same subspace iff each is a nonzero multiple of 
the other. 

Proof. If A r = \B r for some nonzero A, then clearly a A A r = iff a A B r = 0, so they 
represent the same subspace. Conversely, suppose A r and B r represent the same subspace; then 
A r = ai A a 2 A • • • A a r and B r = b± A &2 A • • • A b r for some linearly independent sets of vectors 
{ai}i=i,..., r and {0j}j=i,...,rj & n d each of the bj is a linear combination of the ai. Substituting 
those linear combinations into the expression for B r , removing the terms where any ai appears 
twice, and reordering the factors in each term, I find that B r equals A r multiplied by some 
scalar. This scalar can't be zero because B r is nonzero, so that completes the proof. □ 

Frankly, something would be wrong if this weren't true. In turn, Theorem [5] gives me another result I'll use 
a lot. 

Theorem 5. If A r represents a proper subspace of A s , then A r can be factored out of A s ; that is, there 
exists a blade A s _ r such that A s = A r A A s — r . 

Proof. Let A r = ai A • • • A a r ; then {aj}j=i,..., r is a linearly independent set lying in A s , so it 
can be extended to a basis {oj}j=i,...,s of A s . That means a± A • • • Aa s defines the same subspace 
as A s , so it differs from A s by a scalar factor; absorb that factor in the newly-added vectors and 
we have A s = A r A A s — r , where A s — r is the outer product of the newly-added vectors. □ 

I'll show later that A r A B s = (—l) rs B s A A r for any r- and s-vectors, so A r can be factored out of A s from 
either side and the other blade doesn't have to change by more than a sign. 

While I'm here, let me also give a necessary and sufficiënt condition for a vector to be orthogonal to a 
subspace. This theorem uses aj A r , the left inner product of a vector and an r-blade, which once again I 
haven't defined yet. For now, think of it as taking the inner product of a with each factor of A r separately, 
as in Eq. (|57|) below. (Now that Fm taking the inner product of objects of different grades, it matters which 
of the two products I use; notice that the "floor" of the product points toward the vector.) 

Theorem 6. If A r is a nonzero r-blade with r>l,a is orthogonal to the subspace A r iff a J A r = 0. 
Proof. To show this, let A r = a\ A a 2 A ■ ■ ■ A a r and let's look at 
a J A r — a J (ai A a 2 A • • • A a r ) 

r 

= ^(-l) j_1 (a J dj) ai A • • • A %• A • • • A a r , (37) 
i=i 

where I used Eq. ([77]) in the second line. (I'll derive it later. The check over aj means it's 
not included in the outer product.) If a is orthogonal to A r then it's orthogonal to all the aj 
and a\ A r = 0. If instcad a is orthogonal to all the aj but one, then a\ A r is the product of 
a nonzero scalar and the outer product of the remaining aj, which is nonzero because they're 
linearly independent. So aj A r ^ 0. The remaining case is a nonorthogonal to multiple aj, in 
which case let {aj}j=i,...,s where 1 < s < r be the vectors for which a J Oj 7^ 0, and for j = 2, . . . , s 
let bj be defined by 

h = <>3 (4^) °i- ( 38 ) 

\a\ai) 
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None of the bj equal because the aj are linearly independent, and a \ bj — 0; better yet, because 
each of the bj is just aj with a multiple of a\ added to it, the outer product is unchanged by 
replacing the üj with the bj: 

A r = ai A ü2 A • • • A a r 

— ai Ab 2 A • • • A b s Aa s+ i A ■ • • A a r . (39) 

Now I'm back to the previous case where only one vector in A r is nonorthogonal to et, so I get 
the same result as before. Therefore if et is not orthogonal to A r then a \ A r ^ 0. □ 

The orthogonal complement of a subspace is the set of all vectors orthogonal to every vector in the 
subspace. Theorem |6] says that the orthogonal complement of A r is the set of all vectors a satisfying 
etj A r = 0. 

Just as with Theorem [3l the proof of Theorem [6] actually shows a bit more. If aj A r ^ 0, then it's an 
r — 1-blade, and it represents the subspace of A r that is orthogonal to a. Hl use this fact later when I 
interpret inner products geometrically. 

So not only do we know that r-blade A r represents an r-dimensional subspace, we have an easy way to 
teil whether vector a is in that subspace (a A A r = 0) or orthogonal to it (a\ A r = 0). Theorems[3]and[ü]are 
also our first examples of a general fact: algebraic relations between multivectors reflect geometrie relations 
between subspaces. We '11 see more advanced examples later. 

Now lct's suppose the scalars are real, in which case blades are also supposed to have orientation and 
weight. To give A r an orientation, I note that it's the product of r vectors in a given order. That order 
dennes an orientation: follow the vectors in their given order, and then follow their negatives in the same 
order until you're back where you started. For example, the orientation of the 2-blade a A b is found by 
moving along a, then 6, then —et, then —b back to the beginning. A little experimentation shows that 
interchanging any two vectors reverses the orientation, and it also changes the sign of the blade. Therefore 
there are two orientations, and changing orientations is associated with changing the sign of the blade or 
equivalently interchanging vectors. 

Now this definition had nothing to do with what the scalars are; the problem with non-real algebras 
arises when you try to decide what scalar multiplication does to the orientation. In real algebras it's easy: 
every nonzero scalar is positive or negative, which either leaves the orientation alone or reverses it. Other 
sets of scalars are not well-ordered like this, so we can't say unambiguously what they do to the orientation 
of blades; this is why I define orientation only for real algebras. 

All that remains now is to define a blade's weight, which should include the notions of length, area, 
and volume, but generalize them to arbitrary dimensions. That's most easily done in Section liOl on scalar 
products, so I will defer it until then. (F11 also explain why I define it only in the real case.) Taking that 
for granted at the moment, I can now conclude that a general multivector is a sum of terms that represent 
different subspaces of all dimensions from on up, with weights and orientations if the scalars are real. This 
of course is what I was after in the very beginning. 

4. The inner, outer, and geometrie products 

Before we go on, Hl repeat what we know about the product of two vectors using all my new notation from 
Section O 

uv = u\v + uAv (40) 

where the two terms are the symmetrie and antisymmetric parts of the product. We also know that the first 
term is a scalar and the second is a bivector. That means the inner and outer products can also be written 

u J v = (uv) 

uAd = (uv) 2 . (41) 

My next job is to extend this knowledge to the inner, outer, and geometrie product of any two multivectors 
at all. Once that's done, I will also geometrically interpret the inner and outer products. Along the way Hl 
build up a set of tools and identities Hl use later to do calculations. 
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4.1. The inner, outer, and geometrie products of a vector with anything 

As a steppingstone, I first define the inner and outer products of a vector with any multivector. Sometimcs 
I will want to explicitly indicate that a particular vector is absent from a product; I do this by including the 
vector anyway with a check'over it. For example, if {tii}i=i,...,r is a collection of vectors, then a\02 ■ ■ ■ cij ■ ■ ■ a r 
is the product of all the vectors except for aj . 

I define the left inner product of vector a and r- vector A r (also called the inner product of a into A r ) to 

be 

a\A r :=-[aA r -(-l) r A r a], (42) 
and I also define the right inner product of A r and a (or the inner product of A r by a) to be 

A r [a := - [A r a- (-l) r aA r ] 

= (-l) r - 1 a\A r . (43) 

(Just as in Section[31 the "floor" of the inner product always points toward the vector. Later Fll show how to 
calculate it the other way.) When r = 1, 1 recover the inner product of vectors defined earlier. An equivalent 
way to write the relation between these products under interchange is 

a\ A + = -A + la while a\ A- = A- [a. (44) 

Here's why I define inner products this way. 

Theorem 7. a \ A r and A r \_a are both r — 1-vectors, so the left or right inner product with a vector is a 
grade lowering operation. 

Proof. To show this, I start by proving this relation: if a, a\, 02, . . . , a r are vectors, then 

1 r 

— [aai<22 • • • a r — (— l) r a\a2 • • • a r a\ = l) 5-1 ^] aj)aiü2 ■ ■ ■ dj ■ ■ ■ a r . (45) 

j'=i 

I proceed by induction. If r = 1 the result is true because it reduces to the definition of the inner 
product. Suppose the result holds for r — 1, so 

1 1 

-(aaia 2 • • ■ a r _i) = -(-l) r 1 {a\a 2 ■ ■ ■ a r _ia) + 1) J ' 1 {a J aj)aiü2 ■ ■ ■ a? • • • a r _i. (46) 

i=i 

Then since 

— {aa\a2 ■ ■ ■ a r ) — — (aaiö2 • • • a r _i)a r 
1 _ 

= — (— l) r 1 (aia2 • • • a r _ia)a r + 1) J ' 1 (aJ aj)aia2 • • • &j ■ • • a r -\a ri (47) 

we find 
1 



[aaiQ2 ■ ■ ■ a r — (— l) r ai02 • • • a r a] 



r-1 



1 

— (— l) r ~ 1 (aia2 ■ ■ • a r _iaa r + a\02 • • • a, — \a r a) + ^ l)-'~ 1 (a j aj)a\02 ■ ■ ■ dj ■ ■ ■ a, — \a r 



i=i 



(— l) r 1 {a J a r )aia2 ■ ■ • a r _i + 1)^ 1 (a J aj)aia2 ■ • • a>j ■ ■ • a. r _ia r 
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= l)-' 1 (a J aj)a\ü2 ■ ■ ■ dj ■ ■ ■ a r , (48) 

i=i 

which is the desired result. 

Now let's look at the special case that the numbered vectors are an anticommuting set: 

1 r 

-[aeie 2 ■ • • e r - (-l) r eie 2 • • • e r a] = ^(-iy^ 1 (a J e,)eie 2 • • • ëj ■ ■ ■ e r . (49) 

The right hand side is a sum of r — 1-blades, making it an r — 1-vector. Now a generic r-vector is 
a sum of r-blades, and any r-blade can be written eie 2 ■ • • e r , so it follows that for a vector a and 
r-vector A r the quantity | [aA r — (— l) r A r a] , which is the left inner product, is an r — 1-vector. 
Since the right inner product differs only by a sign, it's an r — 1-vcctor too. □ 

This begins to show why we have two inner products: the vector operates on the r-vector to lower its grade, 
not the other way around. Two products allow the vector to do this from either side. Notice that when 
r = (so A r = A) both inner products reduce to 

aj A = A [a = 0. (50) 

In retrospect this makes sense: these products lower grade, and the scalars have the lowest grade in the 
algebra, so there's nothing other than zero for them to be. 

At this point it can't be much of a surprise that the outer product of vector a and r-vector A r is defined 
to be ^ 

aAA r :=-[aA r + (-l) r A r a], (51) 

and the outer product with the order reversed is given by 

A r A a := ^ [A r a + (-l) r aA r ] 

= {-l) r ahA r . (52) 

When r = 1, of course I recover the outer product of vectors. The behavior of the outer product under 
interchange is the opposite of the inner product, and it can also be written 

aAA + = A + Aa whilc sAA_ = -i_Ao. (53) 

This theorem is probably no surprise either. 

Theorem 8. aAA r is an r + 1-vector, so the outer product with a vector is a grade raising operation. 
Proof. To show this, I again need to prove a prcliminary result: 

^ [a(ai A a 2 A • • • A a r ) + (-l) r (ai A a 2 A • • • A a r )a] = a A a\ A a 2 A • • • A a r . (54) 

Again I use induction. If r = 1 the expression reduces to the dchnition of the outer product 
of two vectors, so suppose it's true for r — 1. Let a\ A a 2 A • • • A a r = eie 2 ■ ■ ■ e r where any two ei 
anticommute; thcn 

a A oi A a 2 A • • • A a r = a A e\ A e 2 A • • • A e r (55) 

because the substitution of the for the at yields a factor of det(iï)=±l which can be eliminated 
as before, so the prcliminary result becomes 

i [a(eie 2 • • • e r ) + (— l) r (eie 2 • • • e r )o] = a A e\ A e 2 A • • • A e r . (56) 

To prove it, let's begin by looking at the term on the right hand side, which we know is the 
sum of (r + 1)! permutations. I want to regroup it into r + 1 terms, each of which consists of 
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all permutations that put a particular one of the r + 1 vectors in the first position. Let the 
permutations that put a first be called n, and the ones that put a first be called 71^; thcn 



a A ei A e 2 A • • • A e r 



1 



r + 1 



^ X^( sgn7r ) ae7r ( i ) e ' r ( 2 ) " ■ e * 

r 1 



1 r! 

J = l TT, 



(2)---«7T J (r) 



(57) 



whcre onc of the a 



7r,(i) 



and the others are the e,; other than e,-. Now the e„(j) in the first 



term on the right hand side can be rearranged, canceling the sgn ir factors just as before, so 



r + 1 



~7 X!( sgn7r ) ae7r ( 1 ) e7r ( 2 ) ' ' ' e?r ( r 



r + 1 



aeie 2 ■ • • e r . 



(58) 



As for the other terms, irj as a permutation of a and all the ej that puts e^- in the first spot 
has the same sign as the corresponding permutation of just a and the other ej times a factor 
(— iy , because this is the factor gained by moving ej from its original j + 1 position to the 
front. Thcrcfore with that factor added each Uj may be thought of as a permutation of only the 
remaining r vectors, or 



^ ( s g n ^(i)^ (2) ■ ■ ■ ^ 



(r) 



= (-l) j ej(a Aei A e 2 A • • • A éj A • • • A e r ) 



(— l)- , ej[aeie 2 • • • ëj • • • e r + (— l) r eie 2 ■ ■ ■ ëj ■ ■ ■ e r a] 



= - [(-l)- 7 e :) aeie 2 ■ ■ ■ ëj ■ ■ ■ e r + (-l) r eie 2 • • • e r a\. 



(59) 



In the second line I used the fact that the relation is assumed true for r — 1, and in the second 
term on the third line I moved ej past the first j — 1 ej. Now since 



CjCL — 2a J — ae 



the relation above becomes 



i ^(sgnTr^eja^^!)^.^) • • • a Wj{r) 



(60) 



(61) 



-(-ir 1 a e,e 1 e 2 ...ë r ..e r + (-ir (a je,)e 1 e 2 ...ë r 
-aeie 2 ■■■e r + -(-l) r eie 2 ■■■e r a + (-iy(a\ ej)e Y e 2 ■ 



Putting all this back together, I get 

a A e\ A e 2 A • • • A e r 
1 



e T + -(-l) r eie 2 ■ 



1 



aeie 2 • • • e 



r [1 1 
r + ^2 ( -aeie 2 ■ ■ ■ e r + -(-l) r eie 2 • • • e r a+ 

.7 = 1 



(-l)-'(aj ej)eie 2 ■ ■ ■ ëj ■ ■ ■ e r 



1 r 

r+l + 2r + 2 



aeie 2 • ■ ■ e r + 



2r + 2 



(-l) r eie 2 • • • e r a- 



■ e r a 



(62) 
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\ae\e-2. • • • e r — (— l) r eie2 • • • e r a\ 



2r + 2 

= ^ [oeie 2 • • • e r + (-l) r eie 2 • • • e r a] , (63) 

which proves the preliminary result. 

The right hand side of Eq. (I5H) is an r + 1-blade. Now a generic r-vector is a sum of r-blades, 
and any r-blade is of the form a\ A a 2 A • • ■ A a r , so it follows that for a vector a and r-vector A r 
the quantity | [aj4 r + (— l) r A r o] , which is just the outer product, is an r + 1-vector. And that's 
that. " □ 



You may wonder why I didn't define two outer products. In this case, it makes equal sense to think of the 
vector raising the r-vector's grade by 1 or the r-vector raising the vector's grade by r. 
So now we know that 

aA r = a\ A r + aA A r , (64) 

and further 

a\ A r = (aA r ) r _ 1 

aAA r = (aA r ) r+1 , (65) 

and similar results hold for A r a: 

A r a = A r [a + A r Aa, (66) 

where 

A r la = (A r a) r _ 1 

A r Aa=(A r a) r+1 . (67) 

These expressions generalize Eqs. (|4Ü| and (|4T]) and reduce to them when r = 1. In fact, summing over 
grades r in Eqs. (|64[) and f|66[) shows that they're true for any multivector A. So I've achieved my goal from 
the beginning of this section, at least for the special case of multiplying by a vector. The expressions for J 
and |_ hi Eqs. (jïï5"|) and (jBTj) work even when r = 0, because back in Section [5] I made a point of defining all 
negative-grade components of a multivector to vanish. This is why. 



4.2. The general inner product, outer product, and geometrie product 

So far I have shown that the product of a vector and a multivector is the sum of two terms; if the multivector 
is grade r, the two terms have grades r — 1 and r + 1. I've also shown how to calculate each term separately. 
Now I'll use this information to characterize the product of any two multivectors, and aft er that I'll introducé 
the most general forms for the inner and outer products. 

Let A = J2 r Ar and B = Bs> then AB = s A r B s , so I'll consider each term separately. 

Theorem 9. A r B s consists of min{r, s} + l terms of grades \r — s\, \r — s\ + 2, \r — s\ + 4, . . . , r + s, or 

min{r,s} 

A r B s = ( A r B s)\r-s\+2 3 - (68) 



Proof. If r = or s — this expression is obviously true, so next I'll consider the case 
< r < s, assume A r is an r-blade A r (if it's true for a blade it's true for sums of blades), and 
proceed by induction on r. If r = 1 the expression becomes Eq. (|64|) . which I've proved already; 
so assume it's true for r — 1. A r can be written aA, — i where a is a vector and A r —\ is an 
r — 1-blade, so 

A r B s = aA r —±B s 
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min{r — 

E a ( A r-l B s)\ r -l-s\+2j 
3=0 

r-1 

= E a ( A r-l B s) s - r +2j + l 
3=0 

where the second line uses the fact that the relation is assumed true for r 
follows from the inequality r < s. Now applying Eq. (|64p. 



(69) 

1 and the last line 



r- 1 

A r B s = ^2 [ a \ ( A r 



-xB 



s l s-r+2j + l 



■ah (Ar-iBg 



s-r+2j+l 



3=0 

-- a\ (A T 

r-X 

E 



L-B,. 



[aj 



s-r+2j+l 



aA(A P _iB s ) 8 _ r+2j _ 



(70) 



s-r+2j-l 



(71) 



Noting the grades of the various terms in the sum, I identify 

(A r B s ) s _ r = a\ {A r _xB s ) s _ r+l 

{A r B s ) s _ r+2j = a\ (A r _ 1 B s ) s _ r+2j+1 +ah(A r _ 1 B s 

(A r B s ) r+s =ah(A r -iB 3 ) s+r _ 1 . 

Since s — r = \r — s\, A r B s is now expressed as a sum of terms of grade \r — s\, \r — s\ + 2, 
\r — s\ + 4, . . . , r + s, which proves the result for r. The remaining case is < s < r, which is 
proved by induction on s. □ 

This proof actually gives somewhat explicit formulas for the terms in the product. To illustrate this Fll 
consider Eqs. ([7Tt for the special case r — 2, so A? — e\e2' 

A 2 B S = ei J (e 2 B s ) s _ 1 + ei J (e 2 B s ) s+1 + ei A (e 2 B s ) s _ 1 + e 1 h (e 2 B s ) s+1 

= ex J (e 2 J B s ) + e x J (e 2 A B s ) + ex h(e 2 J B s ) + ex h(e 2 h B s ). (72) 

I could have arrived at the same result by writing A 2 B S = exe 2 B s , using Eq. (|64[) to expand e 2 B s , and using 
Eq. (|64[) again to expand the product of ex with each term. The first term is of grade s — 2, the middle two 
terms are of grade s (one grade lowering and one grade raising operation applied to B s ), and the final term 
is of grade s + 2. 

Theorem|ÏÏ]tclls us something important: while A r B s is not an r + s- vector, every term has grade r+s — 2j 
for some j, so A r B s is even if r + s is even and odd if r + s is odd. That means that the product of two even 
grade elements is itself an even grade element, so the even grade subspace of any geometrie algebra (defincd 
at the end of Section [2]) is not just a subspace but a subalgebra. (This is not true of the odd subspace, 
because the product of two odd elements is also even.) Also, since (— l) r + s ~ 2 J = (— l) r + s , it follows that 

a\(A r B s ) = -(aA r B s - {-l) r+s A r B s a) 



ah(A r B s 



(aA r B s + (-l) r+s A r B s a 



(73) 



for any vector a, which is kinda nice. Four identities follow from this. 
Theorem 10. 

a\{A r B s ) = (aj A r )B s + (-l) r A r {a\ B s ) 
= (ahA r )B s - (-l) r A r (ahB s ) 
ah(A r B s ) = (ahA r )B s - (-l) r A r (a\ B s ) 
= (a\A r )B s + (-l) r A r (ahB s ). 



(74) 
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Proof. These are all proved the same way, so I'll show only the first one. Starting with the 
first of Eqs. (1751) and then adding and subtracting |(— l) r A r aB s , 

a\(A r B s ) = ^(aA r B s - {-l) r+s A r B s a) 

= ^(aA r B s - (-l) r A r aB s ) + h-l) r (A r aB s - (-l) s A r B s a) 

= (a]A r )B s + (-l) r A r {a\B s ) (75) 

where I reassembled the terms into inner products using the first of Eqs. (|73[) again. □ 

By summing over grades s, you can see that these identities are valid even if B is a general multivector. 
(Thcy're also valid for general A with a little twcaking; see Section lOl ) 
An obvious generalization of Eqs. (|73]) is 

a J (ai<i2 • • • a r ) — —(aai 02 ' * * a r ~ ( — l) r aiö2 ■ • ■ a r a) 

a A(aia2 • • • Of) = —(aai 02 ■ ■ ■ a r + (— l) r aia2 • • • a r a) (76) 

and the first of these equations can be used to prove another handy result. (I used this result to prove 
Theorem [üj you may recall.) 

Theorem 11. 

r 

a\(a 1 Aa 2 A---Aa r ) = V^(— l) J ' _1 (aJ dj) a± A a2 A • • • A üj A • • • A a r . (77) 

j'=i 

Proof. Using the first of Eqs. (f76]) and Eq. (|45|) . I can write 

a\ (aiü2 ■ ■ ■ a r ) — — {aa\a2 ■ ■ ■ a r — (—l) r a±a2 ■ ■ ■ a r a) 

r 

= ^^(— l)- ,_1 (a J üj) a\ü2 ■ ■ ■ dj ■ ■ ■ a r . (78) 
j=i 

I'll prove just below that the grade-s term in the product of s vectors is their outer product (see 
Eq. (|85p). so by taking the r — 1-grade term of both sides and using that result, the identity 
follows. □ 

Here's a nice mnemonic for remembering the coefficients in this sum. In the jth term, a acts on aj , so imagine 
that aj first has to be moved to the far left side of the outer product, which requires j — 1 interchanges of 
adjacent vectors and introducés a factor of (— l)- 1 ' -1 . 

I can use the method of proof of this theorem to prove a fact about versors. 

Theorem 12. An r-versor a\a,2 ■ ■ ■ a r is a linear combination of terms, each of which is an outer product 
of some subset of {oj}j=i,. ..,»•• The number of factors in each term is even or odd as r is even or odd. 

Proof. As usual, the proof is by induction. The result is truc if r — 0, 1, or 2, so assume it's 
true for r — 1; then 

a\a2 • • • a r = fli J (02 ■ ■ ■ a r ) + a\ A(a2 • • • a r ) 

r 

= l) J_2 ( a i J a 3 )a2 ■ ■ ■ dj ■ ■ ■ a r + a\ A(a 2 • • • a r ), (79) 

where I used Eq. (|78[) to go from the first to the second line. The first term is a linear combination 
of products of r — 2 of the , so by the 1 — 2 result the first term is a linear combination of outer 
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products. The number of factors in each term is even or odd as r — 2 is even or odd, or as r is 
even or odd. The second term is the outer product of a\ with the product of the remaining r — 1 
vectors. By the r — 1 result, that product is a linear combination of outer products, and each 
term is even or odd as r — 1 is even or odd. When you take its outer product with a±, it's still a 
linear combination of outer products, and each term is even or odd as r is even or odd. □ 



Now I define for any A — A r and B — J2 S 



A\ B 



r,s 



A[B 



AAB 



(80) 



AU previous expressions for the inner and outer products of two objects are special cases of these definitions. 
Further, this dehntion for the outer product of two objects and the definition for the outer product of 
arbitrarily many vectors in Eq. (|30[) are consistent with each other thanks to Eq. (|54p. which shows that 



Some other facts are worth mentioning. 

1 . A.f J Bf = A.y [ B>f = {Aj°By°) . 

2. If r > s then A r J B s — B s \_A r = because all negative-grade multivectors vanish. 

3. The lowest grade term in A r B s is A r J B s if r < s and A r [ B s if r > s. 

4. The highest grade term in A r B s is A r A B s . 

5. For any A, XA = \\A = \f\A and AX = A[X = AAX. The product of a vector and any multivector 
is a sum of inner and outer products, as shown by Eqs. (|64j) and (|66j) . In all other cases there are 
additional terms of intermediate grades. 

These definitions make the inner and outer products much easier to work with, because in general the 
geometrie product has nicer algebraic properties. The main advantage of the inner and outer products over 
the full product is nice behavior under interchange of the factors; as you'11 see in Section RT21 A r B s and 
B s A r are related, but not in a way that's easy to use, while going from A r J B s to B s [ A r or from A r A B s 
to B s A A r is just a matter of a sign change (Eqs. (|141[) and (|143[1 ). 

The definitions also allow me to deduce three more identities from Eqs. ([74")) ; I take the r + s — 1-grade 
term of the first of Eqs. (|74p . the r — s + 1-grade term of the third, and the s — r + 1-grade term of the last, 
with the results 



(Taking an appropriate-grade term of the second identity only yields a special case of the third of Eqs. ([55)1 
below.) Again, these expressions are actually valid for general B, and also for general A when I use some 
results from Section [5Tl 

Like the geometrie product, the inner and outer products are distributive, and they obey these identities. 
Theorem 13. 



a A(oi A ü2 A • • • A a r ) — a A ai A ü2 A • • • A a r . 



(81) 



a J (A,. A B s ) 
aA(A r [B s ) 
aA(A r J B s ) 



(oj A r )AB s + {-l) r A r A(a\B s ) 
(aAA r ) [B s - (-l) r A r L(aJ B s ) 
(a\A r )\Bs + (-l) r A r \(aAB s ). 



(82) 



AA(BAC) = {AAB) AC 
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A](B[C) = (A\B)[C 
A\(B\ C) = {A/\B)\ C 
A[(B AC) = (A[B)[C 



(83) 



So the outer product and certain combinations of left and right inner products are associative, but neither 
left nor right inner products are associative by themselves. In the homogcncous case, the first rclation abovc 
becomes 

A r A{B S A C t ) = (A r B s C t ) r+s+t . (84) 

An important special case of this is 

a± A ö2 A • • • A a r = (ai<X2 • • • a r ) r ■ (85) 
Associativity of the outer product plus its properties under interchange of factors leads to the result 

aAAAb=-bAAAa (86) 
where a and b are any vectors and A is any multivector. 
Proof. For the first relation, I note that 

A r A(B s ACt) = A r A{B s Ct) s+t 

= (A r (B s C t ) s+t ) r+s+t 

= (A r B s C t ) r+s+t 

= ((A r B s ) r+s C t ) r+s+t 

= (A r AB s )AC t . (87) 

The crucial step is taken on the third line, where the ( ) s+t is dropped. This can be done because 
the only term in A r B s Ct that has grade r + s + 1 is the term that comes from multiplying A r by 
the highest grade term in B s C t . From this it follows that AA(B AC) = (AAB) AC for any A, 
B, and C. The third line also gives me Eq. For the second relation, consider 

A r \(B s [C t ) = A r \ (B s C t ) s _ t 

= (A r (B s C t ) s _ t ) s _ {r+ty (88) 

Now this vanishes automatically unless r < s — t, in which case 

(A (B s C t ) s _ t ) s _ (r+t) = (A r B s C t ) s _ (r+t) (89) 

because when the inequality is satisfied, the only term in A r B s Ct that has grade s — (r + t) is 
the term that comes from multiplying A r by the lowest grade term in B s Ct- Therefore 

A r \(B s [C t ) = (A r B s C t ) s _ {r+t} 
= {A r B s C t ) {s _ r) _ t 
= {{A r B s ) s _ r C t ) {s _ r) _ t 
= (A r B s ) s _ r [C t 

= (A r \B s )[C t . (90) 

This expression also vanishes unless the inequality is satisfied, so A\(B |_C) = (A\ B) [C in 
general. Finally, 

A r \{B s \C t ) = A r \(B s C t ) t _ s 
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= (A r (B s C t ) t _ s ) t _ {r+s) . (91) 

Now this vanishes automatically unless r + s < t, in which case logic similar to that used above 
yields 

A r \{B s \C t ) = (A r B s C t ) t _ [r+s) 

= ((A r B s ) r+s C t ) t _ (r+s) 
= (A r B s ) r+s \C t 

= {A r AB s )\C t . (92) 

This final expression also vanishes unless the inequality is satisfied; therefore AJ(BJC) = 
(AAB)\ C for all A, B, and C. A similar proof shows that (A [ B) [ C = A [(B A C). " □ 

It's useful to introducé an order of operations of these products to cut down on parentheses. The order 
is outer products, followed by inner products, followed by geometrie products. Thus, for example, 

A\ BACD = {A\(BAC)}D. (93) 

Despite this convention, I'll occasionally put the parentheses back in for clarity. However, I will use it a lot 
in Section[7]on projections, rotations, and reflections. (This is not the only convention in use; the other one 
reverses the order of inner and outer products. I picked this one.) 

Now that I've defined two separate inner products, I should probably say why Generally, we associate 
inner products with projections, but that turns out not to be quite right. Look back at Eq. ([5]) for the 
projection of v along u; noting that the inner product of a scalar into a vector is actually their product, 
you'11 see Eq. ([5]) can also be written 

P u (v) = (v\u)\u- 1 . (94) 

So orthogonal projection is actually a doublé inner product. The geometrie meaning of a single inner 
product, as I'll show in Section FOl is this: for blades A r and B s , A r J B s is also a blade, and it represents 
the subspace of vectors orthogonal to A r and contained in B s . So the roles played by the factors in the 
inner product are not the same; that's why the product is asymmetrie and there are two of them, so either 
factor can play either role. 

Finally, I can use the ncw dcfinition of the left inner product to generalize Thcorcm llll 

Theorem 14. If r < s then 

B r J(ai Aa 2 A • • • Aa s ) = y"](-l)£J=ife~j) (ff r j ^ a a i2 A ■ ■ -Aa. lr )a ir+1 A ■ ■ ■ A a is , (95) 

where the sum is performed over all possible choices of {ffli<}j=i,...,r ou t °f {a>i}i=i,...,s> an d i n each term i\ 
through i r and i r +i through i s separately are in ascending order. 

The coefficients in this sum can be remembered using the same trick used for Theorem II 11 In that case, 
you imagine that you need to permute each vector to the far left in order to act on it with a. For this 
theorem, you imagine that you need to permute each distinct subset of r vectors to the far left, keeping 
them in their original order, in order to act on them with B r . In both cases, permuting vectors to the far 
left introducés a power of —1 equal to the required number of interchanges of adjacent vectors. 

Just as in Theorem II 11 each inner product in the sum is a scalar. Since each term picks r elements out 
of a set s elements, the sum has r) terms. 

Proof. If the result is true for an r-blade, it is true for any r-vector, so assume B r is an 
r-blade. Now I'll proceed by induction on r. If r = 1 this becomes Theorem [TT1 so assume it's 
true for r — 1. Now B r = a A B-, — i where a is a vector and -B r _i is an r — 1-blade, so using the 
third of Eqs. (|83|) I can write 

B r J (ai A a 2 A ■ ■ • A a s ) =(aA -Br-i) J («i A a 2 A • • ■ A a s ) 
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= a J [B r -x J (ai A a 2 A • • • A a s )] 

= H(" 1 ) 5:j " ll(i3 ^ )Bl - 1 J( at i Aa 4a A---Aa ir _ 1 ) x 

aj [a ir A • • • AfljJ (96) 

where i\ through i r -\ and i r through i s are in ascending order separately. Now I use Theorem 
[TTIto get 



B r J( ai Aa 2 A---Aa s ) = ^(-l)^= ll(ij ^-B^-i J (ofe A Oj. 



■ A a* 



(-l) fe r a J a lfe a lr A • • • A a lk A ■ ■ • A a. 



(97) 



This expression has (s — r + 1) ( = rfj terms, which is too many by a factor of r, so it's time 
to do some grouping. 

First, notice that each term is a scalar calculated using r of the vectors, multiplicd by the 
outer product of the remaining s — r vectors arranged in ascending order, and that every possible 
choice of r vectors occurs. That means that for some choice of scalars C(üi 1 , . . . , a ir ), 

B r J [a\ A a 2 A • ■ ■ A a s ) = ^ C(a ll , . . . , a ir ) a lr+l A . . . A a is (98) 

where the sum is over all choices of r vectors out of the set of s. All that remains is to figure out 
what the coefficients C{ai 1 , . . . , üi T ) are. Well, for a given choice of through ai T , the coefhcient 
will include terms in which one of the is in the inner product with a while the others are in 
the inner product with B 1 — 15 or 

r 

C{a h , . . . , a ir ) = ^(-l) ej B r _i J (a h A • • • A A • • • A a ir ) a\ a tj (99) 
i=i 

for some value of ej for each j. This sum has r terms, one for each of the a ij , which is exactly the 
number I need, so now I need to figure out the exponents ej. Remember the mnemonic Fve been 
using: each vector's contribution to ej equals the difference between its position in the original 
outer product and the position to which it is moved to compute the inner product. (You can 
verify that this is true for every vector in Eq. ([97jl by inspection.) So all we have to do is figure 
out those positions. First considcr every at k in the inner product with B r —± where k < j: each 
one is moved from position ik to position k, so it contributcs ik — k to ej. Now consider the a,i k 
where k > j: each one is moved from position ik to position k — 1 (because position j is empty), 
so it contributes ik — k + 1. 

Finally, let's take a look at ct^. . If its inner product is taken with o, then Eq. ((97)) tells me that 
it is part of the second group of vectors. Therefore vectors a ij+1 through a ir had to be moved to 
its lcft, moving it from its original position ij ahead to ij — j + r. It is then moved to position r 
for the inner product with a, so its contribution to ej is ij — j. Therefore 

j— 1 r 

ej = y](ik - k) + ^ {h-k + l) + ij - j 

k=l k=j + l 

r 

= J2( i k~k) + r-j. (100) 
fe=i 

Putting this in Eq. (f9"9"| gets me 

r 

C{a n a ir ) = 5^(-l)SUi(**-*)+'--i J B r _ 1 j ( üil A • • • A ij A • • • A a ir ) a\ a tj 
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= (-l)EJ= 1 fe-fe)(_i)r-i Br ,_ 1 j ^(_i)i-i( a j ai .) aii A . • • A a % A- • • A m 

= (_l)E^ 1 fe-fc)(_l)r-l Br ,_ 1 j [a A • • • A OjJ] 

= (-l)S= 1 fe-^(-l)'- 1 (B r ._ 1 A o) J K A • • • A a ir ) 
= (-l)S=ife-')( a AB P _0 J (a n A • • • A a ir ) 

= (_i)EU(h-*)b p J A • • • A a ir ). (101) 



Comparing this expression for C(at 1 , 
proved the result. 



, a,i r ) with the statement of the theorem, I see that I've 

□ 



It's easy to verify that 



E&-i) = 



(102) 



so this is another way to write the exponent of —1 in the statement of the theorem. Pil use this later. 

Now I can answer some questions left hanging in Section [F3l I asked how you could teil what nvnT 1 



and v A \ 



are without knowing how they were derived. (Notice Fm using order of operations to drop 



parentheses.) First let's do nvn~ x \ since this is proportional to nvn, Fll look at that instead. We know now 
that the product of three vectors will in general by the sum of a vector and a trivector, and the trivector 
is the outer product of the factors. In this case the outer product is nAuAn, which vanishes because n 
appears twice, so the product must be pure vector. 

Next let's look at «Ann -1 ; again Fll consider vAnn because the answer will be the same. The easiest 
thing to do here is expand the product with the final n into inner and outer products: 



v A nn = v A n [n + v A n A n 
= v A n |_ n 



(103) 



because the v AnAn term vanishes. Even if you didn't know the remaining term had to be a vector, you 
could figure it out because it starts as a vector and has one grade raising and one grade lowering operation 
applied to it. 

I also asked how to calculate {y Ann~ 1 ^ J n to verify that vAnn" 1 really is perpendicular to n. (Re- 
member that • has changed to J since we got through Section [2]) Since the inner product of vectors is just 
the scalar part of their geometrie product, 



(v Arm l )\n — (vAnn 
= (u Ara) 
= 



(104) 



since uAn has no scalar part. Easy, huh? 



4.3. The geometrie meaning of the inner and outer products 

Fve now accomplished all I set out to do in this section except for geometrie interpretation. First Fll handle 
the outer product. 

Theorem 15. Let A r and B s be nonzero Mades where r, s > 1. 

(a) A r A B s = iff A r and B s share nonzero vectors. 

(b) A r AB s , if nonzero, represents the direct sum of the corresponding subspaces. 
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Proof. This is true because, by the same reasoning used in Theorem [31 A r A B s = iff the 
factors of A r and B s form a linearly dependent set, which is true iff the subspaces share a 
nonzero vector. The fact that nonzero A r A B s represents the direct sum of A r and B s follows 
immediately, since a A A r A B s = iff a is a linear combination of the factors of A r and B s . □ 

So forming the outer product is equivalent to taking the direct sum, and it's nonzero iff the direct sum can 
be taken. 

Next let's move on to the inner product. I already said at the end of Section [4751 that the inner product 
combines inclusion in one subspace and orthogonality to the othcr, and I need some terminology to more 
conveniently describe this. In Section [3] I defined the orthogonal complement of a subspace: it's the set of 
all vectors orthogonal to every vector in the subspace, or equivalently the orthogonal complement of A r is 
the set of all vectors a satisfying ctj A r — 0. Here's another definition: the orthogonal complement of A r 
in B s is the intersection of B s and the orthogonal complement of A r . Algebraically, a is in the orthogonal 
complement of A r in B s iff a \ A r = and a A B s = 0. Now I'm ready to prove the result. 

Theorem 16. Let A r and B s be nonzero blades where r, s > 1. 

(a) A r J B s — iff A r contains a nonzero vector orthogonal to B s . 

(b) If r < s then A r J B s , if nonzero, is an s — r-blade representing the orthogonal complement of A r in 



Again, this is why two different inner products are defined; the geometrie roles played by the two factors 
in the product aren't the same. In contrast, Theorem 1 1 51 shows that the roles played by the two factors in the 
outer product are the same, which is why there's only one outer product. This also explains geometrically 
why A r J B s = when r > s; as I show in the proof, if one subspace is higher-dimensional than other, the 
larger subspace always contains a nonzero vector that is orthogonal to the smaller one. 

Proof. First I'll consider the case r < s; try to be surprised that the proof is by induction 
on r. The r = 1 result is taken care of by the proof of Theorem [6] (see both the theorem 
and the discussion right after the proof), so assume the results have been proved for r — 1 and 
consider A r J B s . To prove part (a), let a\ be any vector in A r ; then for some A r —! I can write 
A r = A r —i A ai, which means 



Suppose A r contains a vector orthogonal to B s ; then let that vector be a%, so a\ J B s = 0, so 
A r J B s = 0. For the converse, assume A r J B s = 0; then it follows that A r -i J(ffliJ B s ) = 0. 
There are now three possibilities. The first is a\ J B s = 0, in which case a\ is orthogonal to B s 
and I'm done. If not, then by the r — 1 result A, — i contains a vector ai orthogonal to a\ J B s ; 
if that vector happens to be orthogonal to all of B s then Fm also done. Now for the third case: 
a% is orthogonal to ai J B s but not B s . The proof of Theorem [S] showed that I can factor B s as 
b\ A • • • A b s where a\ J i>i 7^ while the other bj are orthogonal to a\, so 



A r J B s = A r -! A fli J B s = A r -! J (o X J B s ). 



(105) 



a 1 \B s = (01 J 6i)& 2 A- • • Ab s . 



(106) 



The only way 02 can be orthogonal to a\ J B, 
to all the other b j . In that case consider 



but not B s is if a-i J b\ ^ while 02 is orthogonal 




(107) 



This vector lies in A r ', it's nonzero because ai and are linearly independent; and it's orthogonal 
to all the bj, including 6 1; by construction. Thus it's orthogonal to B s . So in all three cases 
A r contains a vector orthogonal to B s . To prove part (b) , assume r < s and A r J B s = 
A, — 1 J (ai J B a ) 0. Then by the r — 1 result A r J B s is the space of all vectors in a% J B s that 
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are orthogonal to A r —±. However, by the r = 1 result a± \ B s is the orthogonal complement of 
a± in B s , so a vector lies in A r J B s iff it lies in B s and is orthogonal both to a\ and to A 1 — i, 
and thus to all of A r . This proves both parts for r < s. 

Now for the case r > s; A r \ B s = automatically, so I can forget about part (b) and 
I only need to show that A r always contains a vector orthogonal to B s . Consider B s J A r ; 
either it vanishes or it doesn't. If it doesn't, then it contains vectors in A r orthogonal to B s 
and I'm done. If it does vanish, then B s contains vectors orthogonal to A r ; let B p represent 
the subspace of all such vectors. If p = s, then B s is orthogonal to A r and Fm also done. If 
p < s, then B s = B p A B s — p where B s — p contains no vectors orthogonal to A r , which implics 
B s — p \ A r =/= 0. Then any vector in _B s _ p J A r lies in A r and is orthogonal to _B s _ p ; but by 
lying in A r the vector is already orthogonal to B p , so it's orthogonal to all of B s and the result 
is proved. □ 

Incidentally, Hl show in Section 15721 that A r IB S = (— l) r ^ s_1 ^B a J A r for any r- and s-vectors, so the 
geometrie interpretation of the right inner product is the same as the left product, but with the factors 
reversed, which certainly seems reasonable. 

As a fun exercise, at this point you might look back at the identities in Eqs. ([55]) in the special case that 
A, B, and C are blades and try to hgure out the geometrical meaning of each one. 

To explain the next few theorems, I need some facts about blade A to be proved in Section EU 

1. A 2 is a scalar. 

2. A is invertible iff A 2 ^ 0, and A^ 1 = A/A 2 . Therefore A and A~ x represent the same subspace. 

3. A is invertible iff the inner product is nondegenerate on A. 

The inner and outer products of blades can shed light on how their subspaces are related. For example, 
if one subspace lies inside another, their blades are related as follows. 

Theorem 17. Let A r and B s be nonzero blades where 1 < r < s. 

(a) If A r is a subspace of B s , then A r B s = A r \ B s . 

(b) The converse is true if either (1) r = 1 or s = 1 or (2) A r or B s is invertible. 



Proof. First let r = 1; then aB s = a\ B s iff aAB s = 0, which is true iff a belongs to B s . 
Now assume the result is true for r — 1 and let A r be a subspace of B s ; I can write A r = A r —\a 
for some vector a such that a and A, — i are orthogonal, so the r = 1 result Iets me write 

A r B s = Ar-^aBg = Ar-xa J B s . (108) 

Now a J B s is the orthogonal complement of a in B s , which means a\ B s contains A, — i, so by 
the r — 1 result 

A r B s = A r _ ± a\ B s = A r _ 1 J (a J B s ) 
= (Ar- X Aa)J B s 

= A r \B Sl (109) 

which is the desired result. (If aj B s — then A r \ B s = also by Theorem [T6l so if one side 
vanishes then so does the other.) Fve already shown the converse is true when r = 1, and s = 1 
implies r = 1, so assume A r is invertible; then A r B s = A r J B s implies B s = A~ 1 A r \ B s . By 
assumption B s ^ 0, so A r J B s ^ also. Therefore if r = s, then B s is just a nonzero multiple 
of A~ x . Since A r and A~ l represent the same subspace, A r and B s also represent the same 
subspace. If r < s, then since B s is an s-vector, A~ 1 A r J B s must be too; but only its highest 
grade term, its outer product, has grade s, so for this relation to hold the product must equal the 
outer product, so B s = A~ Y A(A r J B s ). Therefore B s is the direct sum of A r J B s and A^ 1 ; 



29 



but A^ 1 represents the same subspace as A r , so A r is obviously a subspace of B s . The proof 
when B s is invertible is similar; A r B s = A r J B s implies A r = A r J B s Bj x , so if r = s, A r is 
a nonzero multiple of _B S . If r < s, then since Ap is an r- vector, A r J BgB^ 1 must be too; but 
only its lowest grade term, the inner product, has grade r, so the product must equal the inner 
product, so A r — (A r J B s ) J B^ 1 . Therefore, A r is a subspace of -B" 1 , and thus of B s . □ 

Anothcr possible relationship is orthogonality: two subspaces are orthogonal if every vector in one is 
orthogonal to every vector in the other. In that case, their blades are related as follows. 

Theorem 18. Let A r and B s be nonzero blades where r,s > 1. 

(a) If A r and B s are orthogonal, then A r B s = A r AB S . 

(b) The converse is true if either (1) r = 1 or s = l or (2) A r or B s is invertible. 



Proof. To begin, I note that A r can be written a\ü2 ■ ■ ■ a r where the are orthogonal to 
each other; similarly, B s can be expressed &1&2 ■ • • b s where the bj are also orthogonal to each 
other. Now suppose that A r and B s are orthogonal; then all of the and bj are orthogonal to 
each other as well, so now I can use the rule that the product of orthogonal vectors equals their 
outer product to get 

A r B s = a\ai ■ ■ ■ a r b\b2 ■ ■ ■ b s 

= a\ A 02 A • • • A a r A b\ A b% A • • • A b s 

= A r AB s . (110) 

To prove the converse, first let r = 1; then aB s = a/\B s iff a\ B s = 0, which is true iff 
a is orthogonal to B s . Now assume A r is invertible and let A r B s = A r AB s ; then B s = 
A~ x A r A B s . Since A~ x represents the same subspace as A r , which is a subspace of A r A B s , 
it follows that B s = A~ x J (A r A B s ), so B s is orthogonal to A r . The proof when s = 1 or B s 
is invertible proceeds similarly. □ 

You may be surprised that the converse parts of these theorems aren't generally true; let me give an 
examplc to show why not. Consider an algebra with orthogonal vectors ei, e2, 03, and such that is null 
but the others aren't. (These four vectors define only a subspace of the full space of vectors, or Axiom[5]would 
be violated.) Let A3 = eie2e3 and Bs — 626364; then e\ = implies A3B3 = 0, so A3 J B3 — A3 A B3 = 
also. However, neither A3 nor B3 is a subspace of the other, so the converse part of Theorem [TT] doesn't 
hold, and the subspaces are not orthogonal to each other (they both contain non- null 03), so the converse 
part of Theorem 1181 doesn't hold either. Null vectors make life hard sometimes. 

Incidentally, A r B s = A r \B s implies B s A r — B s \_A r for any r- and s- vectors, and the same is true if 
the inner product is replaced with the outer product (this follows from Eq. (|140|) ). so the last two theorems 
don't depend on the order of A r and B s . Really, it would be weird if they did. 

Let me end with a result that combines the previous few theorems in an interesting way Suppose A r is a 
subspace of B s ; then it seems plausiblc that B s should be the direct sum of A r and its orthogonal complement 
in B s . (For example, three-dimensional Euclidean space is the direct sum of the z axis and its orthogonal 
complement, the xy plane.) Using our theorems, that suggests something like B s = A r A(A r \ B s ). Now 
that can't be right as it stands because the result shouldn't depend on the weight of A r , just its attitude. 
That's easy to fix, though: maybe B s = A r A(A T T 1 J B s ) instead. That turns out to be right, but with some 
caveats. Let me prove a theorem I need first. 

Theorem 19. If 1 < r < s and nonzero blades A r and B s satisfy A r B s = A r J B s , then A^.B S — 
A r A(A r \ B a ). 

Proof. Believe it or not, this result can be proved without induction. Suppose the condition 
is true; then 

A 2 r B s = A r (A r B s ) = A r (A r J B s ). (111) 
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Since A 2 r is a number, the left hand side is an s- vector; so the right hand side must be also. Using 
the same logic as in the proof of Theorem II 71 the product on the right hand side must equal the 
outer product, and that proves the result. □ 

To interpret this theorem, suppose A r is a subspace of B s , so A r B s = A r J B s by Thcorem [TT] Now 
either A r is invertible or it's not; first suppose it is. Then A 2 , ^ 0, so taking the result of Theorem [TÖl and 
dividing by A 2 yields B 3 = A r A(A r T 1 J B s ), which is the result I was after. Notice I had to assume A r was 
invertible to get this, though. What if it isn't? In that case A r contains a nonzero vector orthogonal to all 
of A r , so one of two things can happen: either that vector is also orthogonal to all of B s , so A r J B s = 0, 
or it isn't, in which case that vector also lies in A r J B s , so A r A(A r J B s ) = 0. Sincc A 2 = 0, the theorem 
nicely includes both of those cases. 

5. Other operations 

Our computational powers have grown by leaps and bounds, but they're not yet complete. Some more 
operations will be useful to us later, so I'U describe them all here. Please be aware that different authors use 
different symbols for some of these operations; Fve listed all my symbol choices in Appendix [Al 

5.1. Grade involution 

The first has the formidable name grade involution. It is represented by an * and defined as follows. 



A* 
a* 
(AB)* 
(A + B)* 



= X 
= —a 
= A*B* 

= A*+B*. (112) 



This operation takes reflection of vectors through the origin (the second linc in the defmition), sometimes 
called the parity operation, and extends it to the whole algebra. From these rules it follows that 

(0102 • • • a r )* = (-l) r aia2 ■ • • a r , (113) 

which implies 

A* r = (-l) r A r , (114) 
so grade involution leaves even grades alone while changing the sign of odd-grade multivectors, or 

A* = (A) + -(A)_. (115) 

This is equivalent to the occasionally handy result 

{A) ± = l -{A±A*). (116) 

Eq. ([TT!) tells me that 

(A)* r = (A*) r , (117) 
so grade involution commutes with taking the grade-r part, and 

A**=A (118) 

for any multivector A. (That's what makes it an involution.) Suppose A is invertible; then since (A~ 1 )*A* = 
(A~ 1 A)* = 1* = 1, 

(A- 1 )* = (A*)- 1 . (119) 
By projecting onto terms of appropriate grade, the third rule in the defmition becomes 

{A\ B)* =A*\B* 
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(A[B)* = A*[B* 

(AAB)* = A* AB* . (120) 

Formulas with factors of (— l)' r are usually simplified by grade involution. For example, Eqs. (H2l and (jÏÏTj) 
for the inner and outer products of vector a with multivector A become 

a\ A = l;(aA - A*a) 

aAA = ^(aA + A*a), (121) 
the definitions of A [a and AA a from Eqs. (|4"3")l and (|52l) become 



Ala = ~a\ A* 

AAa = aAA*, (122) 
and finally the identities from Eqs. (IT4")) and (|82p can now be written 

a\{AB) = (oj A)B + A*(aJ fl) 

= (aAA)B — A*(aAB) 
aA(AB) = {aAA)B - A*{a\ B) 

= {a\A)B + A*(aAB). (123) 

a\{AAB) = (aj i)AB + A*A(aJ B) 
aA(A L-B) = {aAA) [B- A* [(a \ B) 

aA{A\B) = {a\A)\B + A* \{aAB). (124) 

Every once in a while I'll use * r to indicate grade involution taken r times. A* r equals A if r is even and A* 
if r is odd. 

This operation is called "inversion" by some authors, but that can be confused with the multiplicative 
inverse, which would be bad because both are important and are used frequently (sometimes at the same 
time). 

A blade and its grade involution represent the same subspace since each is a multiple of the other by 

Eq. 

5.2. Reversion 

The second operation is called reversion or taking the reverse and is represented by a t. It's a little more 
complicated, and it's defined as follows. 

At :=A 
a' := a 
{AB)ï := BU f 

{A + B)ï := A' +B ] . (125) 

From this it follows that, for example, 

(a,\(i2 ■ ■ ■ a r )t — a r ■ ■ ■ a^ai, (126) 

which shows that the reverse of any multivector is found by writing it as a sum of blades and reversing the 
order of the vectors in each blade. Hence the name. This also implies 



A» = A (127) 
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for any multivector A. So reversion is also an involution. 
Let {&i\i=i,...,r be an anticommuting set; then 

(eie 2 • • • e r ) t = e r ■ ■ ■ e 2 ei 

= (_l)^-i)/2 ei62 ... er (128 ) 

because r(r — l)/2 interchanges are needed to return the vectors to their original order; therefore 

At = {-l) r{ - r -^' 2 Ar. (129) 

If you evaluate this expression for different r, you quickly find (a) multivectors two grades apart behave 
oppositely under reversion, and (b) two adjacent grades behave the same under reversion ifT the lower grade 
is even (scalars and vectors, for example). Since the effect of either reversion or grade involution is to change 
signs of a multivector grade by grade, these operations commute: 

A't = A f *. (130) 

Eq. (|129p also shows that 

(A)t = (A%, (131) 

so taking the reverse commutes with taking the grade-r part, just as grade involution does. Suppose A is 
invertible; then since (A~ 1 ) t A t = (AA -1 )' = 1+ = 1, 

(A- 1 )! = (At)" 1 . (132) 

By projecting onto terms of appropriate grade, the third rule in the definition bccomes 

{A\ B) f =B+ L^ f 

(A A B)^ = A A*. (133) 

Among other things, this shows that the right inner product can be defined in terms of the left inner product 
and reversion and is thus technically redundant. Oh well. 
Here's an easy application of Eq. (|129[) : 

(Aö) J , = (-l)^- 1 )/ 2 ((Aö)t) r 

_ (_x)Kr-i)/2 (sUt> r (134) 

with the nice special case 

(AB) = (StAt) . (135) 

The even more special case where A and B are homogeneous is also useful. The product of two homogeneous 
multivectors of different grades doesn't have a scalar part, so (A r B s ) = (B s A r ) — if r ^ s, and when r = s 
we get 

(A r B r ) = (BtAl) 

= (B r A r ) . (136) 
Therefore (A r B s ) — (B s A r ) in general. Since the geometrie product is distributive, I can go all the way to 

(AB) = (BA) (137) 

for any A and B, or better yet 

(AB ■ ■ ■ CD) = (DAB • • • C) , (138) 
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so the scalar part of a product is cyclic in its factors. This is very useful. Eqs. (I135[) and (|137|) together also 
imply 

(AB) = (A^B^) . (139) 

Almost every identity involving reverses is proved by successively applying Eq. (1129[) and using the rules 
in Eqs. (|125[) . For example, let's examine a generic term in A r B s : 



= (- 


.1)(H 


-s-2j)(r+s-2j- 


- 1)/2 ((^) f ), +s - 2j 


= (- 


.1)(H 


-s-2j)(r+s-2j- 


1)/2 ( S ^4) r + 5 -2j 


= (- 


.1)(H 


-s-2j)(r+«-2j- 


" 1 )/2(_l) r ('- 1 )/2(_ 1 ^(s-l)/2 


= (- 






-2j ' 



(140) 

So multiplication may not commute, but A r B s and B s A r aren't totally unrelated; term by term, they're 
actually equal up to signs. (Also, successive terms of A r B s , whose grades differ by 2, have opposite behavior 
undcr reversion, which I expected given what I said after Eq. (|129|> .) This result also has two important 
special cases. First, suppose r < s; then Eq. (|140|) with j — r refers to the lowest grade term, so 

A r \B s = (-l) r ^B s [A r . (141) 

Notice that this expression also holds when r > s because both sides vanish. So the inner product of an 
odd-grade multivector into an even-grade multivector anticommutes (with left changing to right and vice 
versa), as in 

A-\B+ = -B+[A-, (142) 

but in all other cases it commutes. Without any restrictions on r and s, Eq. (|140p when j — gives for the 
highest grade term 

A r A B s = {-l) rs B s A Ar, (143) 
so the outer product of two odd-grade multivectors anticommutes like so, 

A_ AB_ = -B- A A-, (144) 

with all other cases commuting. (These last few results are equivalent to Eq. (j 133[) . by the way.) 

The properties of objects under reversion are sometimes helpful in sorting out their grades. As an 
example, let me reconsider the product nvn of three vectors from Section ll.3l Notice that (nvn)^ = nvn. 
Now vectors don't change sign under reversion but trivectors do. Therefore nvn has no trivector component 
and is pure vector. 

A blade and its reverse represent the same subspace since each is a multiple of the other by Eq. (|129p . 
5.3. Clifford conjugation 

The third involution in a geometrie algebra is called Clifford conjugation or taking the Clifford conjugate. 
It's represented by a * and defined as follows: 

A*:= A 
a} := — a 
{ABf := B t A t 

{A + Bf := A t +B X . (145) 
If this looks like a mixture of grade involution and reversion, that's because it is; in fact, 

A t = A*K (146) 
This immcdiately tells us that Clifford conjugation really is an involution, 

A n = A, (147) 
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that it commutes with taking the grade-r part, 

(A)* = (A*) rt (148) 

that when A is invertible 

(A-y = (A*)-\ (149) 

and finally 

(A\ Bf = B x [A t 
(A[B)i = B i \A t 

{AABf = B x AA*. (150) 
The Clifïord conjugate of an r-vector is given by 

At — 4*t 

= (-l) r (-l) r < r - 1 >/ 2 j4 r 

= (-l) r{r+1 ^ 2 A r . (151) 

This looks a lot like reversion. If you evaluate this for different r you find that (a) multivectors two grades 
apart behave oppositely under Clifford conjugation, just as with reversion, but (b) two adjacent grades 
behave the same under Clifford conjugation iff the lower grade is odd, not even (vectors and bivectors, for 
example). So Clifford conjugation resembles reversion with grades shifted by 1, so to speak. 

A blade and its Clifford conjugate represent the same subspace since each is a multiple of the other by 
Eq. (HST). 

5.4. The scalar product 

Next is the scalar product, defined by 

A* B := (A*B) . (152) 

(Some authors define A * B = (AB). I'll teil you why I don't shortly.) First consider the scalar product of 
homogeneous multivectors. Only the lowest-grade term in the product (the inner product) has any chance 
of being a scalar, so it's certainly true that 

A r * B s = (At \B S )= (At [ B.) . (153) 

Since the scalar product and inner products are distributive by construction, it follows that 

A*B=(AÏ\B) = (AÏ[B) (154) 

for any multivectors. Now of course we actually know a little more than Eq. (|1 53|) Iets on. Only the product 
of two equal- grade homogeneous multivectors has a scalar part, so 

A r * B s = (At J B s ) S rs - {Al [ B.) S rs - (155) 

Therefore, homogeneous multivectors of different grades are orthogonal under the scalar product. That 
means the scalar product of two general multivectors may also be written 

* B r 

r 

= J2 A i\B r = J24-[B r . (156) 

r r 

This also makes it clear that 

A*B = A**B*. (157) 
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My results for reversion and Clifïord conjugation also establish some properties of this product; for example, 

A*B = B*A = AÏ*BÏ=A t *B t (158) 

The next to last equality shows that an equivalent definition of the scalar product is (AB>\. The scalar 
product interacts with the other products we know this way. 



Theorem 20. 



A * (BC) = (B^A) * C 
A*(B[C) = (B 1 " [A)*C 
A*(B\ C) = (Bt AA) * C 

A*(BaC) = (B^ \A)*C (159) 



Proof. The hrst idcntity is proved as follows: 

A * (BC) = (A^ BC) 

= (A^B)UC 

= (B t A) * C. (160) 

The remaining three are proved roughly the same way, so I'll prove only the fust one. Using 
Eqs. (TÏ56| and (fT33"l) and the second of Eqs. (|83| . 

A*(B[C) = (AU(B[C)) 
= ((A^\B)[C) 
= (Aï \B)UC 

= (B^[A)*C. (161) 

□ 

The last of Eqs. (j 159[) is the basis for a different approach to geometrie algebra, followed for example 
in [5]. You start by defming the outer product and the scalar product; then you decide you'd like to be 
able to factor the term B out of expressions like A * (B AC). You do this by defining an inner product 
that obeys the last of Eqs. (| 1 59|) . Then you define the geometrie product of two vectors to be the sum of 
their inner and outer products, and you're off and running. This has the advantage that it starts with two 
products that have clearly separated geometrie functions: the outer product builds subspaces out of vectors, 
and the scalar product carries all the metric information. It's thus more congenial to the point of view 
inherent in diffcrential forms, which are built using only an outer product and which clearly separate metric 
and non- metric properties. Personally, I think it's cleaner to start with the fundamental product and define 
every other product directly in terms of it, which is why I follow the approach given here. 

I use the scalar product to define the magnitude or norm of a multivector by 

\A\ 2 :=A*A. (162) 

A is said to be null if \A\ 2 = and a unit multivector if \A\ 2 — ±1. (Despite the notation, \A\ 2 can be 
negative. In fact, \A\ 2 can be all sorts of things, since the scalars aren't necessarily real numbers.) Eqs. (I157[) 
and (|158p imply 

\A\ 2 = \A*\ 2 = l^l 2 = |A*| 2 . (163) 

I define \A\ n for other powers n in the obvious way as a power of \A\ 2 , but due care should be taken that 
the power in question is well-dcfincd (for example, be careful if \A\ 2 is negative). 
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The squared magnitude of a scalar or vector is just its square, and that result can be generalized a bit. 
Suppose A is an r-versor, so it is a product a\ai • ■ ■ a r ; then 

|yl| 2 = A* A = (A*A) 

= ((ai<22 • • • a r )^a\a2 • ■ ■ a r ) 
— (a r ■ ■ ■ a2aiOia2 • • • a r ) 

= a 2 a 2 2 ---a 2 r . (164) 

Thcrcfore |A| 2 is the product of the squares of its factors. (This is why I included the reverse in the definition.) 
Notice also that if A is a versor then A^A also equals \A\ 2 . This gives me a couple of useful results. 
First, versors can be factored out of scalar products in an interesting way. 

Theorem 21. Versor A and gentral multivectors B and C obey 

{AB) * (AC) = (BA) * (CA) = \ A\ 2 B*C (165) 
Therefore if either A or B is a versor, 

\AB\ 2 = \A\ 2 \B\ 2 . (166) 



Proof. 

(AB) * (AC) = ((AB)ÏAC) 
= (B^A^AC) 
= (\A\ 2 B^C) 

= \A\ 2 B*C (167) 

and a similar argument using Eq. (|138[) proves the other part of the equation. The second part 
follows by setting B = C. □ 

Second, versors are easy to invert. 
Theorem 22. Versor A is invertible iff it's non-null, its inverse is given by 

A ~ l = 1^2' ( 168 ) 

and the squared norm of the inverse is given by 

\A- X \ 2 = \A\- 2 . (169) 



Proof. If \A\ ^ 0, then clearly Eq. (|168p gives an inverse of A, so A must be invertible. 
Conversely, suppose A is invertible; then there exists a B such that AB = 1 . Then it follows that 
ö t A t = 1 also, so 

1 = B^A ] AB 

= \A\ 2 B^B, (170) 

so \A\ ^ 0; thus a product of vectors is invertible iff it's non-null, and its inverse is given by the 
above expression. For the squared norm, just calculate |^4 _1 | 2 using Eq. (|168[) . □ 
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An r-blade A r is a special type of r-versor, so these theorems apply to blades too. But for blades, a 
few more results are also true. Since A\ = (~l) r ^ r ^/ 2 A r , \A\ 2 = A^A becomes |A r | 2 = (-l)K'- 1 )/ 2 A 2 r . 
Therefore the norm of an r-blade differs from its square at most by a sign. That means unit r-blades also 
satisfy A 2 = ±1, although that ±1 may not be the blade's squared norm. It also follows that the inverse of 
A r equals the additional expressions 

= (-1)^-1)/^ = (m) 

So the inverse of an r-blade is a multiple of the original r-blade, just as with vectors. Therefore they represent 
the same subspace. 

In Section [1731 I asked how you would calculate the inverse of 2-blade a A b. Well, now we know: divide 
the original blade by its square. I actually calculated (a/\b) 2 in Section [TTT1 and the result was —a 2 b 2 sin 2 9. 
Therefore 

(«Ar' = ^. (172) 
arcr sin o 

By the way, this is also the reverse of a A b divided by its norm squared, as it should be. 
Next, I give a geometrie property of null blades. 

Theorem 23. A nonzero blade is null (and thus noninvertible) iff the inner product is degenerate on the 
subspace it represents. 

Proof. A r — e\e2 ■ ■ ■ e r is null iff e 2 — for at least one i, in which case is orthogonal to 
every vector in the span of {e-,}^^...^, which is just A r . That means that either (a) e, — or 
(b) ei ^ but the inner product is degenerate on A r . Since A r =/= 0, none of the vanish, so 
that leaves case (b): the inner product must be degenerate. Therefore nonzero A r is null iff the 
inner product is degenerate on A r . □ 

So every nonzero blade is invcrtiblc in a Euclidean space, while in non-Euclidean spaces things aren't as 
simple. 

And here's an interesting property of products of versors. 

Theorem 24. If r, s > 1 and nonzero versors A r and B s satisfy A r B s — ; then both versors are null. 

Proof. This one is easy: if A r B s = 0, then |A r | 2 B s = A* r A r B s — also. Now since B s 
is assumed nonzero, it follows that \A r \ 2 = 0, or A r is null. Going back to A r B s — and 
multiplying from the right by B\ establishes that B s is also null. □ 

This also means that the product of a non-null versor and any versor is nonzero. 

A special case of this theorem arises if a vector a both lies in blade A r (a A A r = 0) and is orthogonal 
to it (aj A r — 0). In that case aA r = 0, and we say a annihilates A r . The theorem tells us that a must 
be a null vector and A r must be a null blade, which is clear from Theorem [22] since the existence of such a 
vector makes the inner product degenerate on A r . 

To conclude this section, I assume the scalars are real so I can define the weight of an r-blade, as I said 
I would back in Section [3J The weight is supposed to be a higher-dimensional generalization of volume, 
and one way to get that is the following: express A r as a product of r orthogonal vectors, and define the 
weight to be the product of the lengths of those vectors. Then the weight is the volume of an r-dimensional 
parallelepiped that spans the correct subspace. That's what the norm gives us, as Eq. (|164|1 shows, so I 

define 

weight(A r ) := ^\\A r \ 2 \. (173) 

The extra | | is under the square root because, as I've repeatedly mentioned, the squared norm can be 
negative. By this definition when r = 0, the weight of a scalar is its absolute value. This definition only 
works on scalars for which an absolute value and square root are defined, which is why Fm defining it only 
for real algebras. 

When I get to integral calculus on geometrie algebras, I'll be using the weights of blades not to define 
the theory but to interpret parts of it. Thus integration will be defined on any geometrie algebra, but some 
of its meaning will apply only to real algebras. Since all of our applications will be on real algebras, I think 
we'11 be fine. 
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5.5. The dual 



The next operation is called a duality transformation or taking the dual. Let A r be an invcrtiblc r-blade; 
then the dual of any multivector B by A r is B \ A~ x . (Duality gets its own symbol only in a special case, 
which I'll describe below.) To understand what taking the dual does, let B be a s-blade B s . 

1. If s > r, the dual of B s vanishes. 

2. If s = r, the dual of B 3 is a scalar which is zero iff B s contains a vector orthogonal to A r (Theorem 

3. If s < r, the dual of B s is either zero or an r — s-blade representing the orthogonal complement of B s 
in A r (Theorem [TÏÏl again) . If B s was inside A r to begin with, the dual of B s is just BgA^ 1 (Theorem 
UZI). 

The dual of an arbitrary B is a sum of these results. Duality transformations are useful both for taking or- 
thogonal complements of blades (based on the observations above) and for performing orthogonal projections 
into subspaces, as I'll show in Section 17.11 

Although one can take the dual by any invertible blade, one class of blades is by far the most important: 
those that represent the entire vector space. The dual by these blades is very useful and also has simpler 
properties than the dual in general. 

Let the dimension of the vector space be n; then all n-blades either vanish identically (if the factors 
are dependent) or represent the same subspace (namely the whole space); therefore Theorem @] says that 
all n-blades are multiples of one another. Since the inner product on the whole space is nondegenerate by 
AxiomO Theorem[22]says that all nonzero n-blades are also invertible and thus non-null, so I define a volume 
element I to be a unit n-blade. This determines I to within a sign. (Some people call a volume element a 
pseudoscalar, but I won't.) In fact, I can calculate \I\ 2 explicitly. Let {ei}j = i r be an orthonormal basis, 
and suppose p of the e, square to — 1 while the rest square to 1. Let I = e±e2 • • • e n ; then using Eq. (|164p . 

|J| 2 =e?ei---4 = (-l)*\ (174) 

Therefore \I\ 2 = 1 in any Euclidean space, while \I\ 2 = —1 in Minkowski space. Since |J| 2 = Pi, this 
implies 

j2 = (_ 1 )n(n-l)/2+p < ( 175 ) 

Given Theorem 1231 we can now see that Axiom [5] is just another way to say "volume elements are 
invertible." I could have used instead a weaker axiom that implies only "volume elements are nonzero," and 
that would have been enough to prove some foundational results, like this one that I've been promising for 
some time. 

Theorem 25. The outer product of linearly independent vectors is nonzero. 

Proof. Let A r be the outer product of linearly independent vectors. Since volume elements 
are nonzero, A r must lie in a subspace rcprcscntcd by a nonzero blade A s : then by Theorem [5] 
there exists a blade A s _ r such that A r A A s _ r = A s . Thus A r is a factor of a nonzero blade, 
so A r is nonzero too. □ 

This theorem is actually equivalent to "volume clements are nonzero" because each implies the other. Because 
of this, some authors take this weaker statement as an axiom instead of my Axiom [5] I still like my axiom, 
though, because if I is invertible, taking the dual by I is also invertible. This makes the dual much more 
useful, as you'11 see below. 

Unless otherwise specified, "the dual of A" means "the dual of A by I" and is denoted A . Let's 
reconsider the three ways the dual of a blade can turn out when we're taking the dual by I. 

1. There are no s-blades for s > n, so the first option can't happen. 

2. Any n-blade B n = XI for some A, in which case the dual of B n is just A. 
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3. If s < n, then B s represents a subspace of the full space, so the dual of B s is just B S I 1 . It cannot 
be zero; if it were, then Theorem [Tïï] would say that B s contains a nonzero vector orthogonal to the 
whole space, which Axiom [5] doesn't allow. The theorem also tells me that the dual of B s represents 
the orthogonal complement of B s . 

So the general formula for the dual of multivector A is 

A x := A\ I 1 =AI~ 1 (176) 

and the dual of a blade represents its orthogonal complement. (Hence the choice of symbol.) Taking the 
dual by J" 1 instead of I is the inverse operation; it's denoted by A . Since I 1 = I/I 2 , A and A^ 1 - 
diffcr only by a factor of I 2 , so the duality transformation is its own inverse up to at most a sign. 
Since the product of any multivector with I is an inner product, it's true for any A r that 

A r i = {-ïy^-^iAr. 

= IA* r in - 1 \ (177) 

so for any multivector A, 

AI = IA< n ~ 1 \ (178) 

This has several consequences. 

1. Even multivectors commute with I regardless of the value of n, so their duals can be taken be taken 
from either side with no difference. 

2. In odd-dimensional spaces, the dual of any multivector can be taken from either side. 

3. In cvcn-dimcnsional spaces, the dual of an odd multivector can still be taken from either side, and the 
results differ only by a sign. 

4. The first of Eqs. (|177[) is true even if A r is an r-versor because all terms in A r are even or odd as r is 
even or odd. 

5. In even-dimensional spaces Eqs. (|116|) and (|178|) can be used to separate the pure even and pure odd 
parts of a multivector: 

(A) ± = ^(A± IAI- 1 ) if nis even. (179) 

Duality Iets me prove a surprising result. 
Theorem 26. If the vector space is n-dimensional, every n — l-vector is an n — 1-blade. 

Proof. Let A n _i by an n — l-vector. The dual of is a vector, so A n _i is the dual 

of a vector. But vectors are 1-blades, and the dual of a blade is also a blade, so A n _i is an 
n — 1-blade. □ 

One corollary of this is that in dimensions below four, all r- vectors are actually r-blades. 0- vectors and 
1-vectors are always blades, n-vectors are always blades (which takes care of bivectors in two dimensions 
and trivectors in three), and bivectors in three dimensions are n — 1-vectors and thus blades. 
The dual is distributive over addition, and it's easy to show that 

(AB) 1 - = AB^. (180) 

Taking appropriate-grade terms also shows that 

(AAB) 1 - = A\ B 1 - 

(A\ B) 1 - = AAB^. (181) 
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Thus the dual relates the inner and outer products. (Here's another way to prove these results: start with 
the third of Eqs. (1551) and set C = I^ 1 . That gets you the first equation. Then replace B with B 1 - and take 
the inverse dual of both sides; that gets you the other equation.) A special case of this is (a J A r ) ± = aA , 
which means vector a is orthogonal to subspace A r iff a lies in A^. That's further confirmation that duals 
represent orthogonal complements. It also shows that any subspace has a direct representation (all a such 
that a A A r = 0) and a dual representation (all a such that a J A^: — 0). These two representations are both 
useful in different situations. 

In the discussion around Theorem 1191 I said that if A r is invertible, then any space that contains A r 
is the direct sum of A r and its orthogonal complement. This is certainly true for the whole space, and it's 
nicely expressed in terms of duals. 

Theorem 27. The whole space is the direct sum of A r and A^ iff A r is invertible. 
Proof. 



A r AA^ = (A r \A r ) ± =A 2 r I- 1 = (182) 



r- 

Now Theorem [22] and Eq. (I171[) teil me A r is invertible iff Al ^ 0. So if A r is invertible, 
Eq. (|182l) shows that I is the direct sum of A r and its orthogonal complement; and if A r is not 
invertible, the equation shows that A r and its dual have vectors in common, so they don't even 
have a direct sum. □ 

If A is invertible, so is A x : 

(A ± )- 1 = (AI- 1 )- 1 
i-i 



IA~ l . (183) 



The dual of a grade involution is given by 



= (-ïy'A*^- 1 )* 
= {-\y\Ai- 1 )* 

= (-1HA- 1 )*. (184) 

The dual of a reverse is 

(aY = [(a^] x 

= [P(A^} ± . (185) 
Combining these results, the dual of a Clifford conjugate is 

(A^ 1 - = [jt(A- L )*] X . (186) 

Finally, Theorem [5T] and the second part of Theorem [22] show that the dual almost preserves scalar 
products: 

A^*B^ = (AI- 1 ) * (BI- 1 ) 
= |J _1 | 2 A*B 

= \I\- 2 A*B. (187) 

So taking the dual preserves scalar products up to a scale factor. 

Occasionally it's convenicnt to take the dual by volume elements that aren't normalized. In that case, 
the dual and its inverse may differ by more than a sign, but the difference is still only a scalar multiple. All 
the results in this section are valid regardless of the normalization of I. 
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5.6. The commutator 



The final operation is called the commutator, defined as follows. 

A x B := -{AB - BA). (188) 

Notice the factor of \, which is not present in the usual definition of the commutator, used for example in 
quantum mechanics. The commutator obeys the identity 

A x (BC) — (Ax B)C + B(A x C), (189) 

which is easily verified by expanding out the commutators. This shows that the commutator is a derivation 
on the algebra (it obeys the Leibnitz rule). Use this identity to expand A x (BC) and A x (CB) and take 
half the difference; the result is the Jacobi identity 

A x (B x C) = (A x B) x C + B x (A x C). (190) 

The presence of the second term on the right hand side shows that the commutator is not associative. This 
identity is often given in the cyclic form 

A x (B x C) + B x (C x A) + C x (A x B) = 0. (191) 

From the defining properties of the three involutions it's easy to see that 

(A x B)* =A*xB* 

(A x B)ï = B^ x A^ 

(AxBf=B t xA t . (192) 

The commutator of any multivector with a scalar clearly vanishes, and the commutator with a vector can 
be expressed nicely by decomposing a general multivector as A = (A), + (A)_ and recalling Eqs. ([32"]) and 
l|51p for the inner and outer products. The result is 

axi = aj (A) + + a A (A) _ 

A x a = (A) + [a+ (A)_ Aa. (193) 
This Iets me prove an important result about commuting multivectors. 
Theorem 28. The following statements are equivalent. 

1. A commutes with all multivectors. 

2. A commutes with all vectors. 

3. A = \ + ^l(I)_. 

Item 3 is my sneaky way of saying A equals A in even-dimensional vector spaces and A + fil in odd- 
dimensional spaces. 

Proof. Since scalars commute with everything, I won't mention them again. If A commutes 
with all multivectors then it obviously commutes with all vectors. On the other hand, if A 
commutes with all vectors then it commutes with all blades, since these are products of vectors. 
Therefore A commutes with all sums of blades, and thus all multivectors. 

Now for item 3. The first of Eqs. (|193[) tells me that axI = aAliï the vector space is 
odd-dimensional and a J I if the space is even-dimensional. Now a A I — and a J I ^ for all 
o, because every vector lies in I and no vector is orthogonal to it; therefore all vectors commute 
with I in odd-dimensional spaces and no vectors commute with I in even-dimensional spaces. 
To finish off, let A r be an r-blade where < r < n. If r is even, then a x A r = a\ A r , and if 
this vanished for all a then A r would bc orthogonal to the whole space, in violation of Axiom [5] 
If r is odd, then a x A r — aAA r . Since r < n there certainly exists a vector a outside A r , in 
which case a A A r ^ 0. □ 
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The most interesting of all is the commutator with a bivector. 
Theorem 29. 

A 2 xA r = (A 2 A r } r , (194) 

so the commutator of a bivector and an r-vector is an r-vector; commutation with a bivector is a grade 
preserving operation. 

Proof. To show this, I note that 

A 2 A r = A 2 J A r + (A 2 A r ) r + A 2 AA r 

A r A 2 = A 2 J A r - (A 2 A r ) r + A 2 AA r . (195) 

The first equation is obvious when r > 2; for r < 2, recall that in such cases A 2 J A r = 0. The 
second equation follows from the first because of the properties of the inner and outer products 
under interchange and Eq. (|140j) when j = 1 . Subtracting these equations yields 

A 2 x A r = \{A 2 A r - A r A 2 ) = (A 2 A r ) r . (196) 

□ 

In particular, the set of bivectors is closed under commutation. That means the bivectors form a Lie 
algebra with the commutator serving as the Lie product. That will be important later when I show how to 
use geometrie algebra to describe Lie groups and Lie algebras. 

Since commutation with a bivector is grade preserving, the identity in Eq. (|189[) still holds if A = A 2 
and I replace all geometrie products with either inner or outer products: 

A 2 x (B\ C) = (A 2 x B)\C + B\(A 2 x C) 
A 2 x (B[C) = (A 2 x B) [C + B[{A 2 x C) 

A 2 x (B AC) = (A 2 x B) AC + B A{A 2 x C). (197) 
The last of these relations can be generalized in this way. 
Theorem 30. 

r 

A 2 x (ai A a 2 A ■ ■ ■ A a r ) = V] A a 2 A ■ ■ ■ A(A 2 [dj) A- •• A a r . (198) 

3=1 

Proof. As usual, I use induction. The result is true when r = 1 because the commutator with 
a vector is the same as the right inner product, and the r = 2 result follows from the last of 
Eqs. (|197[) . so let's assume the result is true for r — 1. Then by associativity of the outer product 

A 2 x (ai A a 2 A ■ ■ ■ A a r ) — A 2 x (I? r _i A a r ) (199) 

where £?, — 1 = ai A a 2 A • • • A a r _!. Applying the last of Eqs. (|197p and the r — 1 result yields 

A 2 x (ai A fl2 A ■ • ■ A a r ) — (A 2 x B r _i) Aa r + -B r _i A(A 2 x a r ) 

r-l 

= [ai A a 2 A • • • A(A 2 [ Oj) A • ■ • A a r _i Ao r ] + 

3=1 

ai A a 2 A ■ ■ ■ A a r _i A(A 2 [a r ) 

r 

= ^aiAa 2 A--' A(A 2 [dj) A • • • A a r , (200) 

3=1 

which complctes the proof. □ 

I expand the order of operations to include all of these new operations as follows: perform the involutions, 
then duals, then outer, then inner, then geometrie products, then scalar products, and finally commutators. 
Following this convention, the parentheses in Eqs. (|159j) and the left hand sides of Eqs. (|189[) . (|197|) . and 
(1198|) (but not the right hand sides) may be omitted. 
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6. Geometrie algebra in Euclidean space 

Now let's apply everything I've done so far to some familiar cases. I'll work through the algebras of two- 
and three-dimensional real Euclidean space explicitly, revealing some neat surprises along the way. 

6.1. Two dimensions and complex numbers 

First Fll consider the real plane M. 2 with the Euclidean scalar product; this is often denoted E 2 . It has an 
orthonormal basis {ei, e 2 }, which produces a geometrie algebra spanned by the elements 1, e\, e 2 , and e\e 2 . 
That last element satisfies 

|eie 2 | 2 = (e 1 e 2 ) t eie 2 = e 2 eieie 2 = 1. (201) 

Therefore it qualifies as a volume element I. Since bivectors change sign under reversion, it also satisfies 
I = — 1. It defines a right-handed orientation, and a few examples show that all vectors anticommute with 
I. This is consistent with Eq. ([TTS]) . 

Now for the geometrie products. We know what a scalar times anything and a vector times a vector look 
like; all that remains is the product of a vector and a bivector, or equivalently the product of a vector and 
I. To see what that does, notice that 

lei = — e 2 

Ie 2 = ex, (202) 

so multiplication of a unit vector by I results in an orthogonal unit vector. (Which it should, since multiplying 
by I takes the dual to within a sign.) Eq. (|202[) actually tells us a bit more: left multiplication by I rotates 
a vector clockwise through tt/2. Similarly, right multiplication rotates a vector counterclockwise through the 
same angle. So I 2 = — 1 means that two rotations in the same sense through tt/2 have the same effect as 
multiplying by — 1. Of course, this is true only in two dimensions. 

The even subalgebra of any geometrie algebra is always of interest, so let's take a moment to look at it. 
A generic even multivector can be written Z = x + Iy where x and y are real numbers and I 2 = — 1; in 
othcr words, the even subalgebra of E 2 is isomorphic to the algebra of complex numbers. Now this may be 
a little bit of a surprise, because the even subalgebra represents scalars and areas, while we normally think 
of complex numbers as vectors in the Argand plane. But there's another way to think of complex numbers: 
the polar form z = ré 10 reminds us that z also represents a rotation through angle 9 followed by a dilatation 
by r. How do these two interpretations of complex numbers relate? 

It works out because we're in two dimensions. Then and only then, the even subalgebra is isomorphic to 
the space of vectors; a generic vector in E 2 takes the form z = xe\ + ye 2 where x and y are real numbers, 
and there's a natural isomorphism between the vectors and the even subalgebra of the form 

Z = e\z 

z = eiZ. (203) 

This isomorphism maps a vector in the e± direction onto a pure "real" number, so ei plays the role of the real 
axis. It also maps a vector in the e 2 direction onto a pure "imaginary" number, so ei is the imaginary axis. 
Now think about complex conjugation: it leaves the real part alone while changing the sign of the imaginary 
part. Therefore complex conjugation is a reflection along the e 2 axis, which takes z to z' = —e 2 ze 2 . What 
happens to the corresponding even element? It gets mapped to 

Z' = e x z' 

= -ei(e 2 ze 2 ) 
= -e 1 e 2 (e 1 Z)e 2 
= -Iei(x + Iy)e 2 
= x-Iy 

= Z\ (204) 
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where I used the fact that I anticommutes with all vectors. Therefore complex conjugation corresponds to 
taking the reverse. Now let w and z be vectors with corresponding even elements W and Z; it follows that 

wz = eiWeiZ 

= ei(x + Iy)e 1 Z 
= (x- Iy)Z 

= W f Z. (205) 

Now let's look at this. The right hand side is the product of one complex number with the conjugate of 
another. That has two terms: the real part equals the dot product of the corresponding vectors, while the 
magnitude of the imaginary part equals the magnitude of the cross product of the vectors. The left hand side 
is the geometrie product of the vectors, which is exactly the same thing. One of the goals of geometrie algebra 
was to take the complex product, which combines the two-dimensional dot and cross products naturally, and 
gencralize it to any number of dimensions. (That was a goal of quaternions too. I'll show how that worked 
out in the next section.) 

Now for rotations. An element W of the even subalgebra has a polar form rexp(— 19) for some r and 9. 
Letting r = 1, multiplication by a vector z produces the vector 

z' = Wz 

= exp(— I9)z 

= exp(-I0/2) cxp(-I9/2)z 
= exp(-I9/2)zcxp(I9/2) 

= RzR- 1 (206) 

where I defined R — exp(— 10/2). Clearly R is a rotor, so multiplication by W performs a counterclockwise 
rotation through 9. What is the corresponding transformation of Z? 

Z' = eiz' 
= eiWz 
— wz 
= W^Z 

= exp(J6»)Z. (207) 

Thus vector z is rotated counterclockwise through 9 when the corresponding even element Z is multiplied 
by exp(I9), exactly as you'd expect. 

In conclusion, the complex numbers are the even subalgebra of the geometrie algebra of the Euclidean 
plane; the identification with vectors is just an accident in two dimensions, just as identifying planes with 
normal vectors works only in three dimensions. Now while complex algebra is useful, so is complex analysis; 
we use its techniques to perform many ostensibly real integrals, for example. If geometrie algebra generalizes 
complex algebra to any dimension, then perhaps calculus of geometrie algebras could generalize complex 
analysis too. I'll describe geometrie calculus later, and I'll show how it generalizes the Cauchy integral 
theorem and other useful results. 

6.2. Three dimensions, Pauli matrices, and quaternions 

Much that is truc in two dimensions carries over to three: E 3 has an orthonormal basis {ei,e2,e3}, so its 
geometrie algebra is spanned by 1, e\, e2, e^, e\e<i, e±e3, e2e3, and the volume element e^e^. This volume 
element is also denoted I, dcfincs a right-handed orientation, satisfies \I\ 2 = 1^1 = 1, and squares to — 1. 
Unlike the two-dimensional case, a few examples show that all vectors and bivectors commute with I, as 
required by Theorem 1281 Now I've repeatedly mentioned that only in three dimensions can you identify 
planes with normal vectors, which is why the cross product works there. The map between planes and 
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normal vectors should be a duality transformation, so the cross product should be the dual of something. 
Well, it is. If a and b are vectors, then 

axb={aAb) ± . (208) 

So cross products are easily converted into outer products and vice versa. Yay. This shows why the cross 
product is not associative even though the outer product is; the dual gets in the way. Since duality is just 
multiplication by — I and vectors and bivectors commute with I, I can use Eq. (|208[) to writc 

ab = a\b + Iaxb, (209) 

which is the three-dimensional analog of the product W' Z of complex numbers W and Z . Another popular 
product in traditional vector algebra is the triple product a-b x c, which relates to geometrie algebra by 

a-b x c= (aA&Ac)^. (210) 

This form makes the cyclic property of the triple product obvious. The triple cross product a x (b x c) is 
also pretty common, and it can be expressed in geometrie algebra as 

a x (b x c) = -oJ(6Ac). (211) 

From here, it's easy to see that the BAC-CAB rule for expanding this product is really just a special case of 
Theorem[III 

Fd like to say a little more about cross products. Since (a A b) 1 ^ = a\b^ and in three dimensions b is a 
bivector, it follows that a x b is also the inner product of a and a bivector. You may recall that in classical 
mechanics, linear and angular velocity are related by a cross product: v = uj x r. The geometrie algebra 
equivalent is v — f2 \_r, where f2 = ui is an angular velocity bivector in the instantaneous plane of rotation. 
(Bivectors figure promincntly in rotational dynamics, as Fll show in Section [10.11 ') You may also recall that 
the magnetic part of the Forentz force on a point charge is F = qv x B, where B is the magnetic field 
vector. In geometrie algebra this becomes F = qv \ B, where B — B^ is the magnetic field bivector. Fll 
show later on that the bivector representation of B is more physically motivated than the vector version. 

A consequence of Thcorem 1261 which I mentioned at the time, is that in dimensions under four, every 
r-vector is actually an r-blade. In two dimensions that was obviously true; we had only scalars, vectors, and 
multiples of I. In three dimensions the scalars, vectors, and trivectors are obviously blades (the trivectors 
are multiples of J), and I can show using geometry that all bivectors are 2-blades. Consider two 2-blades 
and B2; each represents a plane passing through the origin, and any two such planes in three dimensions 
share a common line. Therefore A2 — a A b and B2 — a Ac where a is a vector along the line shared by the 
planes. This means that 

A 2 +-B 2 = aA& + aAc=flA(6 + c) (212) 

is also a 2-blade. Thus any bivector in three dimensions is a 2-blade, as Theorem l26l demands. 

Now for the products. As before, we know what a scalar times anything or a vector times a vector looks 
like; next Fll do a vector times a bivector. Fet a be a vector and B be a bivector; then a = au + a± where 
a, II lies in the plane determined by B and a± is perpendicular to it. In that case there exists a vector b 
perpendicular to au such that B = aub, so 

aB — (au + a±)a\\b 
= a|6 + a±a\\b 

= ap + a±_ A 0|| A6. (213) 

So the product of a and B is the sum of two terms: a vector in the plane of B perpendicular to a, and 
the trivector defined by B and the component of a perpendicular to it. Clearly the vector is aj B and the 
trivector is aAB. The trivector can also be written 

±|aj.||o||||6|J = ±|a||||B|/ (214) 

where | | is the magnitude of a multivector defined in Section [5^41 The ± is there because we don't know the 
orientation of the system defined by the three vectors. 
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Vector times trivector is even easier. If a is a vector and T is a trivector, then T = abc where b and c 
are perpendicular to each other and to a, so 

aT = a 2 bc = a 2 bAc. (215) 

So aT is a bivector representing the plane to which a is perpendicular. This is clearly a\T. 

The product of two bivectors looks like this: since all bivectors are 2-blades representing planes, let 
vector o lie along the direction shared by bivectors A2 and B2, so A2 = ba and B2 = ac where b and c are 
perpendicular to a but not necessarily to each other; then 

A 2 B 2 = baac 
= a 2 bc 

= a 2 b\c + a 2 bAc. (216) 

So the product of two bivectors is a scalar plus a bivector representing the plane normal to their intersection 
lino. The first term is A 2 J B 2 = A 2 [B 2 and the second term is A 2 x B 2 - 

Next, a bivector times a trivector: if bivector B = ab where o and b are perpendicular, then there exists 
vector c perpendicular to a and b such that trivector T — bac, in which case 

BT = baabc = a 2 b 2 c, (217) 

so the product of a bivector and a trivector is a vector perpendicular to the plane of the bivector. This is 
also B \ T. 

The product of two trivectors is just a number. In fact, it's the product of the volumes defined by the 
two trivectors, with the sign determined by their relative orientations. 
The general multiplication rule for the basis vectors can be written as 

e i e 3 = % + X! Ie ijk e k, (218) 

fe 

which is exactly the multiplication rule for the Pauli matrices. Therefore the Pauli matrices are just a matrix 
representation of the basis vectors of three dimensional space. It is well known that the Pauli matrices form 
a Euclidcan Clifford algebra, but the idea that they are literally matrix representations of x, y, and z is not 
so familiar. 

Finally, the even subalgebra of the geometrie algebra on E 3 has some surprises for us too. Let the unit 
bivectors be labeled 

-Bi = e 2 e 3 , £? 2 = eie 3 , and B 3 = eie 2 . (219) 
Notice that this definition is not consistently right-handed because of £?2- These objects satisfy the relations 

B\=B\ = B\ = -1 (220) 

and 

B!B 2 B 3 = -1, (221) 

so the even subalgebra of the algebra on E 3 , which is spanned by 1 and the Bi, is isomorphic to the 
quaternions. The quaternions were created to generalize the complex numbers to three dimensions, of 
course, so something like this was expected; but the quaternions as Hamilton conceived them were intended 
to correspond to the three unit directions, not three planes. The map between them works diffcrcntly in 
two dimensions and three, so while complex numbers can be thought of consistently as either vectors or 
bivectors, quaternions can be mapped from one to the other only by introducing an inconsistency in the 
handedness, as I've done here. 
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7. More on projections, reflections, and rotations 



In Sections 11.11 and 11.21 1 introduced projections along vectors, reflections along vectors, and rotations in 
planes. My purpose was to get you interested in geometrie algebra by showing how well it handled all three 
operations compared to traditional vector algebra. Well, there's more. It turns out that these operations 
can be defined on subspaces just as well as vectors; for example, rotating a subspace means rotating all the 
vectors in it. As I'll show, in geometrie algebra this is very easy, and the resulting formulas are almost the 
same as the formulas for vectors. 

7.1. Orthogonal projections and rejections 

Let's restate what I did in Sections 11.11 and 1X721 a littlc diffcrently. Let u and v be vectors; then the orthogonal 
projection of v along u is given by 

P u (v) = v\ uu' 1 = OJ u)\ u- 1 (222) 

and the orthogonal rejection of v from u is given by 

R u ( v ) = vAuu^ 1 = vAu [u^ 1 . (223) 

(The second parts of each equation are easy to verify.) P u (v) is parallel to u, R u {v) is orthogonal to u, and 
Pu(v) + R u (v) = v. These operations require u to be invertible, so it can't be a null vector. I promised in 
Section [l.2l that this would have geometrical meaning, and now we're about to see what it is. 

7.1.1. Projecting a vector into a subspace 

Let's take a moment to consider the general notion of projection into a subspace. Let S be a subspace (S is 
not a blade this time; it really is the subspace itself) and let a be a vector not in S. Then for any v € S I can 
write a = v + (a — v), which is the sum of a vector in S and a vector not in S. So which v is the "projection" 
of a into S? We can't say without further information. For example, consider two subspaces S± and 5*2 that 
share only the zero vector; then if a vector lies in their direct sum, it can be expressed only one way as a 
vector from S\ plus a vector from S2, and thus has a unique projection into either subspace. Projection into 
a subspace is specified not only by the subspace itself but also by the subspace the rest of the vector will 
belong to, and the operation is well-defined only if the two subspaces share only the zero vector. 

Now consider orthogonal projection as an example of this. The idea is to express a vector as a sum of 
two terms, one in subspace S and one in S- 1 , the orthogonal complement of S. This works only if S and 
S have no nonzero vectors in common, which is true iff the inner product is nondegenerate on 5*. Thus 
orthogonal projection is well-defined only for a subspace with an invertible blade. In that case, I get this 
result. 

Theorem 31. If a is a vector and A r is an invertible blade, then the orthogonal projection of a into and 
the orthogonal rejection of a from subspace A r are given by 

P Ar {a) = a\ A r A- r x = (aj A r )\ A; 1 

RA r ( a ) = fl A A r A- 1 = aA A r [ A' 1 . (224) 



Proof. First, it's clear that PA r (a) + i?A T ,(a) = o. Now, a J A r is the dual of a by A~ x ; it is 
zero if a is orthogonal to A ri and otherwise it is an r — 1-blade representing the subspace of A r 
orthogonal to a. In that case its product with A~ x equals its inner product with A~ x , which is 
just the dual by A r ; the result is a vector that lies in A r . On the other hand, aAA r is zero 
if a lies in A r , and otherwise it is an r + 1-blade that contains A r . In that case the product 
with A~ x equals the right inner product, and is the dual of A^ 1 by (aA A r ) _1 , so the result is 
a vector orthogonal to A r . Both formulas give vectors, they sum to a, the first lies in A r and 
vanishes iff a is orthogonal to A r , and the second is orthogonal to A r and vanishes iff a lies 
in A r . Therefore the two expressions are obviously the orthogonal projection of a into and the 
orthogonal rejection of a from A r . □ 
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So projecting into a subspace is the same as projecting onto a vector; you just replace the vector with 
the blade representing the subspace. We'11 see several more examples of this idea below. 
I can demonstrate directly that Pa, (a) hes in A r : 

P A » A A r = (oj A r A~ l ) f\A T = (oj A r A~ l A r ) r+1 = (aj A r ) r+1 = 0. (225) 

Similarly, I can show that (a) is orthogonal to A r as follows: 

R Ar ( a )\ A r = {aAA r A- l )\ A r = (a A A r A~ 1 A r ) r _ 1 = (aAA r ) r _ 1 = 0. (226) 

This result applics to the Gram-Schmidt process for producing an orthogonal set of vectors from a linearly 
independent set with the same span. Let {o,j}j=i,...,r be linearly independent; then we build the orthogonal 
set {bj}j=t,..., r as follows. Let b\ = ai to start with. Then 62 equals 02 minus its projection onto bi, or 
cquivalently the orthogonal rejection of 02 from b\. Next, 63 equals the orthogonal rejection of 03 from the 
span of 61 and 62, and so on through all of the aj. Therefore we proceed as follows. 

1. Let b\ = a%. 

2. For each j starting with 1, let Bj = b\ A • • • Abj. 

3. Then let b j+1 = a j+ i A BjBj 1 . 

This procedure will work only if each Bj is invertible, which is why it is normally used only in Euclidean 
spaces. 

If A r is a blade, then A^: represents the orthogonal complement of A r . That means that orthogonal 
projection into A^ should equal orthogonal rejection from A r . Using Eqs. (| 18 1[) and (|183[) . this is easy to 
show directly. 

a\ A^{A^)- 1 = (ttAA,) 1 ^ 1 )- 1 
= aAA r I- 1 IA- 1 

= aAA r A~ 1 . (227) 

If A r and B s are orthogonal, then the projection of a vector into their direct sum should be the sum of 
the projections into the subspaces individually. (For example, the projection of a vector into the Euclidean 
xy plane should be the sum of the projections onto the x and y axes separately.) This can also be shown 
directly. By TheoremöH A r A B s = A r B s , so using the first of Eqs. (|74p I find 

-Pa, a b b (a) = -Pa,b b (a) 

= a\(A r B s )(A r B s )- 1 

= [(aj A r )B s + (-l) r A r (a\ B S )]B~ S 1 A~ 1 

= a\ ArA- 1 + (-l) r A r {a] B S )B; 1 A; 1 . (228) 

Now let's work on that last term. If A r and B s are orthogonal, A r and a\ B s are too, so their product is 
an outer product, so I can interchange them and piek up a factor of (— l) r ( s_1 ). And since A" 1 and B^ 1 
are multiples of A r and B s , their product is also an outer product, so I can interchange them and piek up 
a factor of (-l) rs . Putting all this in Eq. (|228]) , 

Pa,ab» = a\ ArA, 1 + (-l) r + r ^+ rs (a\ B s )A r A^ 1 B~ 1 
= a\ ArA' 1 +a\ B a B^ 

= PaM + Pb s {o). (229) 
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7.1.2. Projecting a multivector into a subspace 

Now that I can project a vector into a subspace, how about projecting one subspace into another? As I 
suggested above, this seems simple enough: project subspace B s into subspace A r by taking every vector 
in B Sl projecting it into A r , and seeing what subspace you get. The rejection should be similar: just reject 
all the vectors individually. However, if I am a bit more precise, I discover a wrinkle. I dehne 

Pa„ (h A • • • A b s ) := P Ar (h) A • • • A P Ar (6.) 

R At (h A • • • A b s ) := R Ar (h) A • • • A R Ar (b s ). (230) 

Now suppose the set {bj}j=i,.~, s is linearly independent but their projections are not. That would happen 
necessarily if, for example, I projected a plane into a line. In that case, the projection defined this way 
vanishes. Instead of objecting to this wrinkle, I decide that it provides useful extra information. If by 
chance the projections of the members of B s do not form an s-dimensional space, so be it; I accept that the 
projection is zero. 

These formulas make geometrie sense, but they aren't very easy to use. However, they can be made 
simpler, and extended to all multivectors to boot. Here's how. 

Theorem 32. For any invertible blade A r and vectors {bj}j = i s , 

P A MA---AP AT (b s ) - (6iA-.-A6.)J A r A~ l 

RaMA'-ARaM = (hA-'-Ab^AArA- 1 . (231) 



Proof. I start with the first equation. Since each bj = P Ar (bj) + R A ^.(bj), the outer product 
b\ A' ■ ■ Ab s can be written as a sum of terms, one of which equals P Ar .(bi) A • • • A P Ar (b s ) while 
each of the others contains at least one R Ar (bj). Consider what happens to each term when 
you takc the inner product with A r and multiply by A^ 1 . The term P Ar (bi) A • • • A P Ar (b s ) lies 
inside A r , so by Theorem 1 1 71 the inner product becomes a product, so the A r and A~ l cancel out 
and you're left with P Ar (bi) A ■ • • A P Ar (b s ). On the other hand, each of the other terms contains 
a factor orthogonal to A r , so the inner product with A r vanishes. Thus the first equation is 
valid. 

For the second equation, I again write b\ A ■ • • A b s as a sum of terms, but this time I note that 
one of them equals R Ar . (b\) A • • • A R Ar (b s ) while each of the others contains at least one P Ar (bj). 
Consider what happens to each term when you take the outer product with A r and multiply by 
A^ 1 . The term R Ar (pi) A • • • A R Ar (b s ) is orthogonal to A r , so by ThcorcmfTHlthe outer product 
becomes a product, so the A r and A^ 1 cancel out and you're left with R Ar (bi) A • • • A R Ar .(b s ). 
On the other hand, each of the other terms contains a factor that lies in A r , so the outer product 
with A r vanishes. Thus the second equation is valid too. □ 

Therefore P At .(B s ) = B s \ A r A~ x and R Ar (B s ) = B s A A r A^ for any blade B s . Taking the obvious 
step, I define the orthogonal projection and rejection of any multivector to be 

P A JB) := B\ ArA- 1 = (B\ A r ) J A; 1 

R Ar (B) ■= BAArA- 1 = BA A r [A; 1 . (232) 

You might be surprised that both projection and rejection leave scalars untouched: 

PaAX) = R Ar (X) = X- (233) 

This had to happen for reasons I'll explain in Section [9.6l Projecting into and rejecting from I do what you 
think they should (except for that odd bit with scalars): 

Pi(B) = B 

Ri(B) = (B) . (234) 
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Running it the other way around, here's what happens when you project and reject I: 

P Ar (I)=IS rn 

R Ar (I)=I5 r0 . (235) 

Again, this makes sense; only the whole space is big enough to project I into, and only zero-dimensional 
spaces are small enough to reject I from. 

With a little relabeling and rearranging, the first parts of Eqs. (j232j) become 

A\B s =P Bb (A)B s 

AAB S = R Bb {A)B s . (236) 

This shows that the inner or outer product of a multivector and a blade can also be expressed as a geometrie 
product, as long as the blade is invertible so projection is defined. Using Theorem[2T] this also shows that 
the norm squared of A J B s equals the norm squared of Pb b (A) times the norm squared of B s , and a similar 
rcsult holds for the outer product. 

Comparing Eqs. (|222l) and (|223p with Eqs. (|232[) . the level of generality achieved is astounding. Starting 
with the projection of one vector along another, I've shown that the projection of any multivector into a 
subspace is meaningful and is given by the same expression, with the multivector and blade put in place of 
the two vectors. The rejection of one vector from another follows the same pattern. It is true that we've lost 
one property: we no longer have Pa^{B) + RA r {B) — B in general. This makes geometrie sense, however, 
if you look at the proof of Theorem[32j neither a projected blade nor a rejected blade includes all the terms 
that are partly projected and partly rejected, so to speak. 



7.2. Reflections 

To start, I'll review reflections from Section [L2l I defined the refiection of vector v along axis n as follows: 
the projection of v along n gets a minus sign, while the rejection of v from n is unchanged. If the refiection 
is denoted v', then 

v' = -nvn' 1 . (237) 

What I didn't show in Section 11.21 is that reflections preserve inner products, which I'll show now. Using 
Eq. (|237p . the defmition of the inner product, and the cyclic property of the scalar part of a product, 

a'\ b' = (a'b') 

— (nan~ l nbn~ lx ) 
= (abn^ 1 ^ 
= (ab) 

= a\b. (238) 



7.2.1. Reflecting a vector in a subspace 

Just as I used projection and rejection along an axis to define refiection along an axis, I can use projection 
and rejection in a subspace to define refiection in a subspace. The refiection of a in A r is constructed by 
giving the projection of a into A r a minus sign and leaving the rejection of a from A r alone. Using Theorem 
EU Eq. (Tm]), and Eq. (HU, I nnd 

a' ■= -PaM + RaM 

= -a\ ArA^ 1 +ahA r A~ 1 

= -{-lf- l A r [ aA- 1 + {-l) r A r A aA- 1 

= {-l) r A r aA- x 

= A r a* r A-\ (239) 
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In the last line, a* r means a is grade involuted r times; I introduced the notation back in Section 15.11 
(You may wonder why I did this instead of just leaving in the (— l) r . It will make sense in the next section.) 
Another way to arrivé at this formula is to write A r = a±ü2 • • ■ a r and refiect a along each of the aj in 
succession. Once again, an expression in terms of vectors generalizes to subspaces with only minimal change. 
Reflections in subspaces also preserve inner products; the proof is very similar to Eqs. (|238[t . 

7.2.2. Reflecting a multivector in a subspace 

Now that I can refiect vectors, I can refiect subspaces too: the reflection of subspace B s in subspace A r is 
found by taking every vector from B s , reflecting it in A r , and seeing what subspace you get. That would 
mean something like 

(6iA---A6,)' :=6iA---A6^ 

= (A r bf A- 1 ) A • • • A (A-HfA- 1 ) . (240) 

Again, this is geometrically sensible but not easy to use. Fear not; I can fix that. To start with, notice what 
Eq. (|239[) shows: conjugating a vector by an invertible r-blade gives you back a vector. A more general 
version of that is also true. 

Theorem 33. If A is a versor and B s is an s-vector, then 

AB s ti = (AB s rf) s , (241) 

so conjugation by an invertible versor is grade preserving. 

Proof. The theorem is true for versors if it's true for vectors, so I'll look at aB s a. Using 
Eqs. (|6l"|l and (|66|). I can write 

aB s a = (aj B s ) [a+ (aj B s ) A a + (a A B s ) [ a + a A B s A a. (242) 

The first term is grade s — 2, the middle two terms are grade s, and the last term is grade s + 2, 
so Fm done if I can ditch the first and last terms. The last term vanishes because a appears twice 
in the outer product (compare Eq. ([86])), and the first term vanishes because it can be rewritten 
as (-1) S_1 (S S [a) La = (-1) S ~~ 1 B S L(aAa) = 0. Since Theorem l22l tells me that the inverse of 
a versor, if it has one, is its reverse divided by its norm squared, conjugation by an invertible 
versor preserves grade too. □ 

I'll use this to get the result I really want. 
Theorem 34. If A is a versor, then 

(ABA^) A(ACA*) = \A\ 2 A(B A C)Al (243) 

Therefore if A is invertible, (ABA~ 1 ) A (ACA~ 1 ) =A(BAC)A- 1 . 

Proof. The result is true for general B and C if it's true for B s and Ct, and I've already 
shown that versor conjugation preserves grades, so 

(AB.A1) A(AC t Ai) = (AB s AÏAC t Aï) s+t 
= \A\ 2 (AB s C t AÏ) s+t 
= \A\ 2 A(B s C t ) s+t Aï 

= \A\ 2 A(B s AC t )Aï. (244) 
If A is invertible, then \A\ 2 ^ 0, so dividing both sides by \A\ 4 yields the desired result. □ 
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Now for reflections. If B s = bx A- • • Ab s , then 

B' s :=6iA---A6^ 

= (Arb^A- 1 ) A ■ ■ • A (Arb^A- 1 ) 
= A r (b\ r A ■ ■ • A b* & r )A r 1 
= A r (b 1 A---Ab s )* r A- 1 

= A r B* s r A-\ (245) 

Taking the obvious next step, I define the reflection of multivector B in subspace A r to be 

B' := A r B* r A-\ (246) 

So reflection in A r is done by grade involuting r times and then conjugating by A r . This is a little more 
complicated than the reflection of a vector along an axis that we started with, Eq. (|237j) , but not much. And 
of course it reduces to Eq. (|237j) when A r and B are vectors. 
The reflection of I in a subspace is 

ï = Arl^A- 1 

= i-iy^AriA- 1 

= {-l) nr {-l) r{n -^IA r A- 1 

= {-l) r I. (247) 

This makes sense because r directions in the space were reflected. So the orientation changes iff r is odd. 

You may have noticed that I now have two ways to reflect a vector around the origin. The first is grade 
involution, and the second is to reflect the vector in a volume element. Since both operations have been 
extended to the whole algebra in a way that respects products, they ought to be equal not just for vectors 
but for any multivector, or 

A*=IA* n I- 1 . (248) 

To show that this really is true, start with Eq. (| 1 T8[) . grade involutc both sidcs, and use I* = (— l) n J. Then 
multiply both sides by J _1 on the right and voila. 

Finally, I can relate reflection in A r and reflection in A^: . Reflection of vector a in A r gives the 
component in A r a minus sign and leaves the component in A^r alone, while reflection in A^r does the 
opposite. Therefore one reflection should be the negative of the other, or 

A^a*( n - r \A^ )" x = -A r a* r A-\ (249) 

Extending this to general multivectors, I expect 

A ^B< n - r \A^)- x = {A r B* r A^ 1 )* . (250) 

And indeed that's what I find: 

A^B< n - r \A^)- 1 = (-ïy^-^Arl^B^Arl- 1 )- 1 

_ <y -iy^-r)^_ 1 y(n-l) ArBsI -l IA -l 

= (-ïyi-iyArB.A; 1 

= {A r B*JA- l y. (251) 

7.3. Rotations 

After all this work, rotations are fairly anticlimactic. Once again, I start with a review of Section 11.21 I 
showed there that a rotation in a plane is the product of two reflections along vectors in that plane, so 

v' = RvR- 1 (252) 

where R is the product of two invertible vectors, also called a biversor or a rotor. A rotation clearly preserves 
inner products since it's just two reflections in succession, but you can show it directly by an argument very 
much like Eqs. ([2"gg]) , 



53 



7.3.1. Rotating a multivector in a plane 

The rotation of a subspace is as simple to understand as the reflection. In fact, it's the example I started 
this whole section with: you rotate a subspace by rotating all the vectors in it. The argument is identical to 
the argument for reflections: if A r = a\ A • • • Aa r , then A' r = a[ A • • • A a' r . Therefore the rotation by R is 

A' r := a[ A • • • A a' r 

= (RaiR- 1 ) A • • ■ A (RarR- 1 ) 
= R(a 1 A---Aa r )R^ 1 

= RA r R-\ (253) 

The grade inversion of A r is absent because it is performed twice, once for each factor in the rotor. Therefore 
the rule for rotating any multivector is 

A':=RAR-\ (254) 

which is exactly the same as the formula for vectors. 

Since I commutes with even multivectors (Eq. (|178l) ). rotations leave I alone, 

RIR- 1 = I, (255) 

as expected. 

When I first discussed rotations in Section II. 2[ I said that any two axes in the same plane separated by 
the same angle would generate the same rotation. That means that if I take the two vectors in R and rotate 
them the same amount in the plane of iï, the resulting rotor should perform the same rotation. Therefore, 
if R and S are rotors in the same plane, SRS -1 should represent the same rotation as R. You can show 
this directly: R and S are both scalars plus multiples of the same area element, so they commute. Therefore 
SRS- 1 = RSS- 1 = R. 

Every linear transformation of vectors can be extended to the entire geometrie algebra; IT1 describe that 
process later. These three transformations extend in a particularly compact way, but not all transformations 
do. Rotations and reflections behave as well as they do because they are orthogonal transformations, and 
geometrie algebra is particularly well-suited to represent them. In fact, it's a good idea to pause and notice 
just how good a job it does; compare Eq. (|254|) to the increasingly complicated expressions you get when 
you rotate tensors of ever-increasing rank. One of the great strengths of geometrie algebra is its ability to 
extend orthogonal transformations to the whole algebra in such a simple fashion. 

7.3.2. Rotations in three dimensions 

In a real three dimensional space, rotations have an interesting property that is easy to understand using 
geometrie algebra: the product of two rotations is another rotation. If Ri represents the first rotation and 
i? 2 the second, then their product is R — R 2 R\. We lose no generality by demanding that both Ri and 
i?2 are unit rotors; and that means R is a unit even versor. In three dimensions the only even grades are 
zero and two, so R is actually a scalar plus a bivector: R = (R) + (R) 2 - Therefore |iï| 2 = 1 becomes 
(R) 2 + | (R) 2 | 2 = 1. That tells me that R = cos(0/2) - 5sin(0/2) for some and unit bivector B. And in 
three dimensions every bivector is a 2-blade, so B represents some plane, and thus R — exp(— B0/2), which 
is a rotation through in plane B. As soon as I climb the ladder to four dimensions, though, I lose this 
result, because R 2 R\ could have a 4- vector part. 

8. Frames and bases 

Now I'U consider a geometrie algebra Q n in which the space of vectors has finite dimension n. Let {ai}i=i,..., n 
be a basis for the vector space, which I will also call a frame. (The ai are not assumed orthogonal.) Then 
a generic element of the algebra will be the sum of a scalar and terms of the form a^a^ ■ ■ ■ ai r for r < n. 
Theorcm [T2l tells me that any such element is a linear combination of blades made up of the {dj.}; therefore 
the scalar 1 and the blades A a{ 2 A • • • A a ir generate the whole geometrie algebra. I'U now show that they 
actually form a basis, and I'U also show how to calculate the components of an arbitrary multivector in this 
basis. 
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8.1. Reciprocal frames 

Given a frame {ffl»}i=i nj another frame {a ï }i=i,..., n is called a reciprocal frame to the first if it satisfies 

a l \ aj =ö l r (256) 

If such a set of vectors exists, it is a frame because Eq. (I25ö[) guarantees that the a % are linearly independent, 
so they form a basis. To construct such vectors, consider their definition: a? should be orthogonal to all of 
the at except for aj , so an obvious way to make it is to take the outer product of all of the ai except for üj 
and then take its dual, which is what I'll do. 

Let a/v = di A ei® A • • • A d n ; then a^r is a (possibly unnormalized) volume element. (Even though a^r is a 
blade, Fm not denoting it with capital letters or boldface; you'11 see why in the next section.) Then I dcfine 

é := (-iy^ 1 (aiAa 2 /\---Ad l A---Aa n )a^ 1 . (257) 

{a 1 } is a reciprocal frame because, using the first of Eqs. (1181j) . 

dj J a? = (— l) : >~ 1 ai J (ai A ü2 A ■ ■ ■ A aj A ■ ■ ■ A a n a~^) 

= (-l) j ^ 1 (a i Aai Aci2 A • • • A&j A • • • A a„)a^ (258) 

Now if i =/= j then equals one of the other vectors in the outer product, so the whole thing vanishes. If 
i = j, I move et; past the first i — 1 vectors to its original spot, which cancels out the (— l)^ -1 prefactor. 
Therefore 

ai J a? = (ai A • • • A a n )a~^5l 

= 51 (259) 

This definition exactly expresses the geometrical idea I started with; ajv was chosen to perform the duality 
transform because it gets the normalization right. 

Since both {a{\ and {a 1 } are bases for the vectors, any vector v can be written v = ^ i/ct^ or v = ^2 v j aJ ■ 
In fact, it's obvious that v l = v J d 1 and Vj — v J üj , so the components of v on either basis are easily calculated 
using the other basis. Using the definition of a % and the first of Eqs. (|181|) again, I find that 

v l — v J a l 

= (—iy~ 1 v J (ai A ö2 A ■ ■ ■ A&i A - ■ -AanO^ 1 ) 
= (-l)' _1 («Aai A a 2 A • • • A a % A ■ ■ ■ A a^a^ 1 

= (ai A • • • A ai_i A v A a l+1 A ■ ■ ■ Ad„)d^ r 1 . (260) 

Compare this with Eqs. (jT6j) and (fTT)) back in Section [TT2l 

Since v is a vector, the expressions for its components can be written v % — v * a l and v 3 ■ = v * a^-, where 
* is the scalar product. These forms for the components can be generalized a long way, as I'll show in the 
next section. 

The inner product of any two vectors follows easily from their components: 

b\c = Y,b l c j {a l \a J )=Y,b l c\ (261) 

i,j i 

and switching the frames on which I expand b and c gives me an equally valid result in terms of the 
components b % and c,. 

A frame and its reciprocal satisfy a useful identity. 

Theorem 35. _^ _ 

a ' ai = ^2 ai a * = n - ( 262 ) 
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Proof. 

a,i a l = ai J a* + a, A 

i i i 

= 11 + ^8^0'. (263) 
To evaluate the second term, expand a % on the original frame to get a % = X) 7 ( aï J ö J )aj, so 

a, ; A a' = Oj A a* J a J <Zj 
t < \ j 

ij' 

= (264) 

because a % J is symmetrie in i and j while aj A dj is antisymmetric. The proof that a' ai = n 
is the same except for exchanging superscripts and subscripts. □ 



8.2. Multivector bases 

Before I continue, I need some fancy new notation. Let I be a string of indices i\, 12, ■ ■ ■ , i r , and given a 
string I let aj be defined by 

a/ := ttj! A a i2 A ■ ■ • A a^, (265) 

and similarly for a 1 . I will use the symbol N only to refer to the string 1, 2, . . . , n, to be consistent with a^r 
in the previous section. I also allow / to be the "empty" sequence, in which case I define a/ = a 1 = 1. Then 
I immediately know several things: 

1. ai = a 1 = iff the string / contains at least one index twice. 

2. If I and J contain the same elements but in a different order, then aj = (sgna)aj and a 1 — (sgna)a J , 
where a is the permutation that changes I to J. 

3. Theorem IT2l tells me that given a frame {ai} for the vectors, the set {ai} (or {a 1 }) where / ranges 
over all increasing sequences generates Q n . (I is an increasing sequence if i\ < 12 < ■ • • < ir-) 

To show either set forms a basis, Hl use this result. 

Theorem 36. 

a 1 * a,, = 5 1 , (266) 

where 5j vanishes if either I or J repeats indices or if I is not a permutation of J (including having a 
different length), and otherwise it equals the sign of the permutation that takes I to J. 

Proof. If either string repeats indices then both sides vanish, and both sides also vanish 
when the lengths of / and J are different (the right side by defmition, the left side because 
a 1 and aj have different grades); to get the other results, I let I = i\ < %2 < • • • < ir an d 
J = jx < j'2 < " " • < ir and use Eq. (|156[) to find 

a 1 * aj = ((a H A a t2 A • • • A a lr )\a jl A a n A ■ ■ • A a jr )) 
= ((a v A ■ ■ ■ A o' 2 Aa'^aj, Aa h A- •• Aa jr )) 

= (a v A • • • A a 12 A a 11 ) J (a h A a h A • • • A a jr ). (267) 
Now consider the case where the ik equal the jk] using the third of Eqs. (|53^) and Eq. ([77]). 

(a lr A • • • A a 12 A a* 1 ) * (a il A a i2 A • • • A a,: r ) 
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= (a ir A • • • A a l2 ) J [a 11 J (a h A a l2 A • • • A o ir )] 

r 

= (a ir A • • • A a 12 ) J ^(-l) 3- V 1 J a 4j . Oi 1 Ao ij A-"Aa ij A--'Aflj 

= (a ir A'"Aa i2 )J(a i2 A--Aa ir ), (268) 

which can be repeated for the ii term and for each successive term until the final result 

(a lr A • • • Aa 12 Aa h ) * (a n Aa i2 A • ■• Aa ir ) = 1 (269) 

is reached. Now suppose that ik equals none of the ji; then when the evaluation of the scalar 
product as shown above reaches the fcth iteration, all of the a lk J üj l terms will vanish, and so will 
the scalar product. This establishes the result for I and J increasing; the general result follows 
from the properties of a 1 and aj under rearrangement of clements. □ 

From this result it's pretty obvious that for any multivector A, 

A = ^A I a I where A 1 = A * a 1 (270) 
i 

and the sum extends over all increasing sequences I (including the null scquence to piek up the scalar part). 
Therefore, given a frame {ai}, the elements aj form a true basis for the geometrie algebra, and the equation 
above shows how to expand any multivector on this basis. (Incidentally, the roles of the frame {ai} and the 
reciprocal frame {a 3 } can be exchanged in this expansion, just as vectors can be expanded on either set with 
the other used to compute the coefficients.) Since the number of distinct r-blades in a basis for each r is 
(™) , it follows that 

n n ✓ 

dimö" = ^dim^ = ^( n J =2". (271) 
I can also express the scalar product of any two multivectors B and C in terms of their components: 

B*C = ^B I C J {a I *aj) ^J2 B i C ^ ( 272 ) 
ƒ,.ƒ i 

and switching the bases on which I expand B and C gives me an equally valid result in terms of the 
components B 1 and Cj. 

A consequence of all this is the following theorem. 

Theorem 37. Multivector A is uniquely determined by either of the following. 

1. A* B for every multivector B. 

2. (A) and a\ A for every vector a. 



Proof. Part 1 is obvious. In fact, it's overkill; A*B for all B in a basis for the algebra will do. 
By the distributive property, part 2 is equivalent to this statement: if (A) = and all a\ A = 0, 
then A = 0. So that's what FU prove. To do this, I assume A is an r- vector A r ; if the result is 
true for r-vectors then it's true for general multivectors too. If r = or 1 then Fm done, so let 
r > 1. Let {ai} be a frame and {a 1 } its reciprocal frame; then a component of A r on the basis 
defined by {o^} is (a lr A • • • A a 12 A a n ) J A r for some strictly ascending choice of i\ through i r . 
However, 

(a v A • • • A a 12 A a n ) J A r = (a lr A • • • A a 12 ) J (a 11 J A r ) 

= (273) 

since a 11 J A r = 0. So all the components of A r vanish, so A r = 0. 
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This proves the theorem in every finite-dimensional algebra, but it's usually true in infinite- 
dimensional spaces too. In fact, extra structures are usually imposed on infinite-dimensional 
spaces for exactly this purpose, and I will happily assume henceforth that this has always been 
done. □ 

This result can be extended ad nauseum: A is uniquely determined by (A) , (A)-^, and A2 J A for every 
bivector A 2 , and so on. 

While any r-vector can be expanded using the frame r-vectors, it can also be expanded using only the 
frame vectors; but now the coefhcients aren't necessarily scalars. 

Theorem 38. 

VVüjJ A r = Y^a* A(a, J Ar) = rA r , (274) 

i i 

and the same is true if the frame and its reciprocal are interchanged. 
Proof. For the first equality, note that 

a l a t J A r = ^ a l J (a t J A r ) + ^ a l A(a i J A r ) 

i i i 

= V"(a* Aa^ J A r + VV A(oj J A r ) 

i i 
= ^2 A J A r + ^ A ( Gi J Ar ^ 

= ^VA(aiJA-) ( 275 ) 

i 

because I showed in the proof of Theorem 1551 that J^i a% A a L = 0. 

The second equality is true for all r-vectors if it's true for all members of a basis for r-vectors, 
so I have to prove it only on a basis; and I know just the basis to use. Let A r — a J1 A • • • A a? r 
for some increasing sequence of indices; then 

a l A(o< J Ar) = ^ a i A [04 J (a jl A • • • A a jr )] 

i i 

= ^(-l)* -1 ^ J a jk )a l A a n A • • • A a }k A • • • A a J " 

i , k 

= X)(-l) fc ~V* A a J1 A • • • A a 0k A • • • A a J ". (276) 

k 

In each term of the sum, I move a Jk past k — 1 othcr vectors to return it to its original spot, 
which cancels the (— l)*^ 1 factor, so 



^2a i A(a i }A r )=J2 ajl 



A • • • A a 



3t 



= rA r . (277) 

That completes the first half of the proof, and exchanging superscripts and subscripts provides 
the other half. □ 

Just as the original frame has a volume element ün, the reciprocal frame has a volume element a N defined 
in an analogous way: a N = a 1 A • • • A a n . Now a^r and a N have to be scalar multiples of each other, and 
since Theorem [36] shows that ün * a N = 1,1 conclude that 

a N = (278) 



• r )cS 



A quick calculation shows 

l^l 2 = M~ 2 , (279) 

so for real algebras, the weights of the volume elements of a frame and its reciprocal are themselves reciprocals. 

Given a frame {ai}, I have Eq. (|257l) for the members a l of the reciprocal frame, but I don't have an 
equally direct formula for the reciprocal multivectors a 1 ; all I can do right now is take the outer product 
of Eq. (|257|) several times. However, I can get a nicer formula for a 1 using the same logic that got me the 
reciprocal vectors in the first place. By construction, a l is orthogonal to the outer product of all the frame 
vectors except a,, and similarly for a J ; therefore o' Aa J is orthogonal to the outer product of all the frame 
vectors except a.; and dj. Therefore a 1 A a J is dual to a\ A • • • A öj A • • • A &j A • • • A a n . To make this easier to 
write out, for any string of indices I let me define I c to be its ascending complement, so I c includes exactly 
the indices not in I in ascending order. In these terms, a 1 is dual to a/c. To be more precise, I have 

Theorem 39. If I represents ascending i\ through i r , 

a 1 = {-VpU^i-^aic a N \ (280) 

Notice that this includes Eq. (|257j) as a special case when I has only one index i. 

Proof. To prove this, Pil calculate a I a^. 

a I a^ — a 1 J ajv 

= (8 ll A--'Afl v )Ja J v 

= (a ll A---Afl ! - 1 )J(o v Ja JV ). (281) 

To calculate a %T J ün I use Eq. ([77)) : 

a %r J ün = a %T J (ai A • • • A a n ) 

= ^(-1)- J ~ 1 (a lr J Oj)ai A- •• A&j A- •• Aa„ 

3 



Now Eq. (|28T|) becomes 



= (-1) 1 '' ai A • • • A a h , A • • ■ A a n 
= [-ïy^aic. (282) 

a J a N = (-l) i '- 1 (a il A---Aa i ''- 1 )Ja i e 

= (-l) i «- 1 (o <1 A---Ao < - a )J(o <r - 1 Joic). (283) 



When I evaluate a 1 ' -1 J using Eq. ([77|) again, I remove ai r _ 1 and multiply by (— l) ir - 1-1 . 
(This is why I put the indices of I in ascending order; i r is later than so cii r _ 1 is still in 

position i r _i in a^.) Thus each step removes a factor from cln and multiplies by (— 1)*^ , 
with the final result 

o J o A r = (-l) 5: J'= l(i '~ 1) o/.. (284) 

Now I'll take care of the special cases on each extreme: I is empty and I = N. When I is empty, 
Eq. ([28D|) reduces to a 1 = ojva^v = •"■) w hich is correct, and when / = N, I c is empty so a/c = 1, 
so Eq. (|280p becomes 



a 



N 



(_l)E? =I (i-i) a -i 
= (-i)"(' l_1 )/ 2 ajv 

«Af 



|ajv| 2: 



(285) 



where I used Eq. (|171[) in the last step. Since this matches Eq. (|278p . it too is correct and the 
theorem is proved. □ 
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To wrap up this part, let me consider the special case where {ai} is an orthonormal frame, which I'll 
denote {ei}. Then the reciprocal frame is a lot easier to find: it's clear on inspection that e l = = e^~ 2 ej 
fits the bill. Since the frame vectors are normalized, ef = ±1, so é 1 = ie^. For any I, let I m be the number 
of elements in the product for e/ that have negative square; then e 1 = (— l) Im ej, so the multivector basis 
and its reciprocal differ at most by signs. On such a basis, Eq. (|272[) for the scalar product becomes 

B *C = ^2(-l) Im Bi Ci = Y^i-^B 1 ^. (286) 
i i 

If the space of vectors is Euclidean, then I m = for any ƒ , so the magnitude is positive definite, so the entirc 
geometrie algebra is also a Euclidean space under the scalar product. If the space of vectors is non-Euclidean, 
then the algebra has (very) mixed signature under the scalar product. 

8.3. Orthogonal projections using frames 

In traditional vector algebra, the orthogonal projection of a vector into a subspace is given as a sum of 
projections onto a basis for the subspace. Although we don't need to do that in geometrie algebra, we still 
can. Let A r = ai A • • • A a r ; then 

= a\(ai A • • • Aa^A^ 1 

r 

= ^(-l) J ' _1 aJ a J (a i A---Aa i A---Aa r )A r : 1 . (287) 

i=i 

Now A r is a volume element for its subspace, so comparing with Eq. (I257[) shows me that the vectors in the 
sum above are the reciprocal frame to {et.,}, or 

r 

) =^(<zj aj)a>. (288) 
i=i 

Since the reciprocal frame volume element A r = a 1 A ■ ■ ■ A a r equals A r /\A r \ 2 (cf. Eq. (I285p b projection 
using either volume element gives the same result; had I used A r , Fd have ended up with Eq. (|288[) with a,j 
and a J interchanged. 

I can do the same thing with any s- vector B s . Let A r be defined as before; then using Eq. (|95p from 
TheoremfHlI find that 

P A r{B s )=B s \A r A- 1 

= B s J (ai A ■ • • A a^A' 1 

= J^i-l^^^^iB, J a h A- • • A a ls )(a l3+1 A • • • A a^A' 1 

= ( _x)-C— i)/a ^(_x)5:; =1 (i,-i)( Ba j a n A • • ■ A a t J(a l3+1 A • • • A a^A; 1 

= ^(- 1 ) E?=l(i "" 1) ( B J J <Hi A • • • A aJK +1 A • • • A a lr )Ar\ (289) 

where in the next to last line I used Eq. (|102|) . If I now let I be the sequence i\ through i s and use Eq. (|280|) 
for the reciprocal multivector basis, I find 

PaABs) = J2( B 1 J ai )(-l)ZU(y-D ai . A- 1 
i 

= ^(B s *a I )a I . (290) 
i 

This expression is still true if I let the sum run over increasing sequences of any length, since all the additional 
terms vanish. Thercfore for any multivector B 

P aiA -Aa r (B) = J2{B*a I )a I , (291) 
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where the sum runs over all increasing sequences, and the expression is still true if the bases {aj} and {a 1 } 
are interchanged. 

9. Linear algebra 

Now that I've said so much about linear spaces, let's take the next step and put some linear functions on 
them. If U and V are vector spaces with the same set of scalars, a function F : U — > V is said to be linear if 
F(au + j3v) = aF{u) + /3F(v), so linear functions respect linear combinations. Linear functions have a very 
well-developed theory, and they're important all over applied mathcmatics; in fact, when a function isn't 
linear, one of the fust things we do is consider its local linear approximation, the derivative. 

In this section Pil hit the highlights of linear algebra using the tools and perspective of geometrie algebra. 
I'll start by reviewing some basic properties of linear functions, and then Hl introducé the adjoint. Hl use 
it to describe three special types of functions: symmetrie, skew symmetrie, and orthogonal, each of which 
relates to its adjoint in a certain way All three have special forms in geometrie algebra, which Hl consider 
in detail. After that, Hl take a giant step into geometrie algebra proper by showing how to take a linear 
function on vectors and extend it in a very natural way to every multivector in the whole algebra. This is 
where geometrie algebra really starts to shine, because it Iets me see old topics in new and useful ways. For 
example, our old friend the eigenvector will be joined by eigenplanes, eigenvolumes, and more, and Hl show 
how to use them to describe linear functions. (It's immediately clear that a rotation has an eigenplane with 
cigenvalue 1, for example.) Hl also give a very easy and intuitive definition of the determinant, and Hl show 
how easy determinants are to calculate in geometrie algebra. 

I'm going to focus on functions that take vectors to vectors, and their extensions to the whole algebra 
will be grade-preserving. To some of you that might seem rather timid; since a geometrie algebra is itself a 
(big) vector space, why not just jump in with both feet and go right for linear functions from multivectors 
to multivectors, whether they preserve grade or not? Well, of course you can, and we already have; duality 
does that, for example. General tensors will also do that, and Hl consider them in due course. 

9.1. Preliminaries 

If F and G are linear, then so are their linear combinations, and so are their inverses if they exist. If 
F : U — » V, then U and V are called the domain and codomain of F respectively. (Some authors call the 
codomain the range.) F singles out two special subspaces: the kemel of F, or Kei(F), is a subspace of the 
domain consisting of all the vectors that F maps to 0, and the range of F, or Range(i 7 '), is a subspace of the 
codomain containing all the vectors that F maps something to. (The range is sometimes called the image, 
presumably by the same folks who've already used the word range to mean the codomain.) It's suggestive 
to think of Ker(F) as F -1 (0) and Range(F) as F(U). F is one-to-one iff Ker(F) = {0}, and F is onto 
iff Range(.F') = V. The dimension of Ker(F) is called the nullity of F, or nul^F), and the dimension of 
Range(-F) is called the rank of F, or rank(_F). If the dimension of the domain U is finite, the rank-nullity 
theorem says 

rank(F) + null(F) = dim U. (292) 

So if both domain and codomain have the same finite dimension, F is one-to-one iff it's also onto. Therefore to 
show invertibility, you only have to show either one-to-one or onto, and the other part follows automatically. 
If W is a subspace of the domain of F, then the restriction of F to If is well-defined and also linear; 
it's denoted Fw- Since blades represent subspaces, Hl sometimes write Fa for the restriction of F to the 
subspace A represents. 

Oftcn we care specifically about linear functions from U to itself, which Hl call linear transformations 
or linear operators. A pretty popular linear operator on any space is the identity; I denote the identity on 
C/byldc/. 

Since our subject is geometrie algebra, I will assume that all vector spaces have inner products and 
belong to geometrie algebras. With that, recall that Theorem [37] in Section [8] shows that any multivector is 
determined uniquely by its scalar products with all multivectors. Looking only at vectors, that means u is 
uniquely fixed if one knows u * v (or equivalently u J v) for all v. This has two useful consequences. The first 
Iets me reconstruct linear functions. 
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Theorem 40. A linear function F : U —ï V is completely determined by knowledge of F(u) * v for all u £ U 
and v . 



Proof. In the finite-dimensional case, F can be constructed explicitly: F(u) — J2i * F(u)] a i> 
where {ai}i=i,.. n is any frame in V and {a?}j=i,...,n is lts reciprocal frame. Since we know all 
the a l * F(u), we know F(u). We can use any frame for this construction since the set of all inner 
products determines F(u) uniquely for each u by Theorem [571 

In the hifmite-dimensional case, it's not obvious we can perform this construction, but for all 
applications I know of, the space is rigged in some way to allow something like this to be done. 
So Fll assume I can do it. □ 

Any linear transformation F defines a new bilinear product on vectors by aj F(b). The second consequence 
of Theorem |37j Iets me go the other way: I start with the product and define F. 

Theorem 41. If o is a bilinear function from vectors to scalars, there's unique a linear transformation F 
such that u o v = u \ F(v) for all u and v. 

Proof. Again, in finite dimensions the proof is by construction: Let F(u) = X^( a * ° u ) a i- 
This defines a linear function because the product o is bilinear, and it satisfies 



i\F(v)=u\ (^2(a l ov)a)j 



{a l o v) (u J ai) 



(u J ai) a l o v 



= uov. (293) 

Since F is determined by u \ F(v) — u o v for all u and v, this function is unique. 

In the infinite-dimensional case, I will assume that whatever structure is needed to make this 
result true has been added. □ 

This is useful, because any given vector space can support many different inner products, but only one inner 
product at a time can be encoded into a geometrie algebra. This theorem tells me I have a way to use 
the other products if I decide I need to. Also, there are bilinear products that I want to use that can't be 
inner products because they aren't symmetrie. This theorem Iets me include those products too, although 
of course the corresponding function F will have different properties. What those properties are will be the 
subject of future sections. 

9.2. The adjoint 

If F : U — > V is linear, then its adjoint is the unique linear function F : V — > U dcfined by 

F(v) *u:=v* F(u) for all u G U and v eV. (294) 

Notice that F switches domain and codomain compared to F and that * is interchangeable with either J 
or [ in this definition. The adjoint of the identity is pretty easy: Idu — Idu. Theorem I4Ü1 tells me how to 
construct F explicitly in the finite-dimensional case: 

F(v) = Yl[ aUT ( v )] a * 

i 

= [ F ( al ) *v\ ai . (295) 
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(Notice that for this to make sense, v e V while the frame {a^} C U .) The bilinearity of the inner product 
shows that taking the adjoint is itself a linear operation: 



aF + (3G = aF + /3G. (296) 

The relationship between F and F is symmetrie, so each is the adjoint of the other, or equivalcntly 

T = F. (297) 

Suppose F : U -> V and G : V -> W, and let GF : U ^ W denote the composition of F and G. Then 

w * GF{u) = G{w) * F(u) 

= FG(w) * u, (298) 

which tells me that 

GF = FG. (299) 
A special case of this arises if F is a operator on U, in which case F™ is defined for any n and 

F^=(F)™. (300) 

Now suppose F : U -> F is invcrtible, so there's an F -1 : V -> [7 such that F -1 F = Idu and FF -1 = Idy. 
In that case, since FF -1 = F -1 F = Idu = Idy, and similarly F -1 F = Idy, 



(F) -1 = F^ 1 . (301) 
The special subspaces defined by a linear function and its adjoint are related in interesting ways. 
Theorem 42. If F : U -> V is linear, 

Ker(F) = Range(F)- 1 -. (302) 

ƒƒ zn addition U is finite-dimensional, 

rank(F) = rank(F). (303) 

Proof. For the first part, 

u e Ker(F) iff F(u) = 

iff F(u) \v = for all w e V 
iff ujF(v)=0 forallueV 
iff u e Range(F)- 1 -. 

For the second part, we start with the rank-nullity theorem and the result of the first part: 

rank(F) = dim U - null(F) 

= dimt/ - dimKcr(F) 

= dimt/ — dimRange(F) ± . 

Now Range(F) and Range(F)- 1 are duals, so their dimensions add up to the dimension of U; so 
picking up where I left off, 

rank(F) = dimt/ - dimRange(F) 1 - 
= dimRange(F) 

= rank(F). (304) 

□ 
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9.3. Normal operators 

In the next few sections I'll be considering operators that commute with their adjoints: FF — FF. These 
are called normal operators, and they have properties that I'll describe here so I can use them later. 

Theorem 43. Any power of a normal operator is also normal. 

Proof. If F is normal, then F n (F) n is easily transformed to (F) n F n by moving all the F 
factors past all the factors of F. □ 

Theorem 44. F is normal iff F(u) \ F(v) = F(u) J F(v) for any u and v. 

Proof. First assume F is normal. Then 

F(u)\F(v) = FF{v)\v 
= FF(u) J v 

= F{u)\F(v). (305) 

Now assume the relation holds. Then 

FF(u) J v = F(u) J F(v) 
= F{u)\F(v) 

= FF(<u)\v, (306) 
so by Theorem gÖ] FF = ~FF and F is normal. □ 
Theorem 45. F is normal iff F(u) 2 — F{u) 2 for all u. 
Proof. Since 

F(u) J F(v) = - [F(u + vf - F{uf - F(v) 2 ] , (307) 

F(u) 2 = F(u) 2 for all u implies F(u) \ F(v) — F(u) \ F(v) for all u and v. On the other hand, 
F(u) J F(v) — F(u) J F(v) for all u and v implies F(u) 2 = F(u) 2 for all u just by considering the 
case u = v. So squares are equal iff inner products are equal, which takes us back to the previous 
theorem. □ 

Theorem 46. If F is normal and the inner product is nondegenerate on both Range(-F) and Range(_F) 7 
then Ker(i 7 ') = Ker(_F). If in addition the domain of F is finite-dimensional, Range(_F) = Range(-F'). 

Proof. For the first part, 

u € Kcr(F) => F(u) = 
FF(u) = 
=ï FF(u) = 

=> F(u) G Ker(F). (308) 

By Theorem Ker(F) = Range(F)- 1 -, so F(u) € Range(F) x . But wait a second: F(u) € 
Range(F) by definition, and Rangc(F) is nondegenerate, so it must be that F(u) — 0, so u £ 
Ker(F). Therefore Ker(F) C Kei(F). The same argument with F and F interchanged shows 
Ker(F) c Ker(F), so Kei(F) = Ker(F). 

For the second part, since the domain of F is finite-dimensional any subspace and its or- 
thogonal complement are duals, so each is the orthogonal complement of the other. That and 
Theorem 1421 teil me that Range(F) = Ker(F) J - and Range(F) = Ker(F)- 1 . But I just showed 
that Ker(F) = Ker(F), so Range(F) = Range(F) too. □ 

So if the conditions of this theorem are satisfied, normal F and F are both one-to-one (or onto) or neither 
one is. 



64 



9.4. Symmetrie and skew symmetrie operators 

A linear operator is symmetrie if it equals its adjoint, F — F, and skew symmetrie or skew if it is the 
negative of its adjoint, F = —F. The names come from this theorem. 

Theorem 47. The bilinear product a o b :— a \ F(b) is (anti) symmetrie iff F is (skew) symmetrie. 

Proof. Since a J F(b) = b J F(a), it follows that aob~boaiSF = F and a o b = — b o a iff 
F=-F. □ 

Recall that every bilinear product has this form for some F (Theorem |4"T|) . 

Both types of operator are normal, so all the results of Section 19.31 apply to them. (They're all pretty 
trivial in these cases, I have to admit.) Further, every linear operator is the sum of a symmetrie and a skew 
symmetrie operator, because 

F=l(F + F) + ±(F-F). (309) 

Also, for any linear operator F, both FF and FF are symmetrie. 

Powers of symmetrie and skew symmetrie operators are themselves symmetrie or skew symmetrie. 

Theorem 48. Any power of a symmetrie operator is symmetrie. Any even power of a skew symmetrie 
operator is symmetrie, and any odd power is skew symmetrie. 

Proof. Since F n = (F) n , F = F implies F n — F n , so F n is symmetrie also, and F = —F 
implies F n = (— l) ïl F™, so F n is symmetrie or skew as n is even or odd. □ 

The spectral theorem says that every symmetrie F has a frame {ai} of eigenvectors with eigenvalues 
{Ai}, which means F(a) for any a is given by 

F(a) =^A. i (a*a> l . (310) 

i 

Conversely, every F of this form is symmetrie. Analogous results hold on infinite-dimensional spaces with 
various additional restrictions. 

Skew symmetrie operators also have a canonical form, which is expressed very nicely in geometrie algebra. 
As motivation, notice that if F is skew, then a\ F(a) — 0, so F maps any vector to an orthogonal vector. 
Well, I know something else that does that: taking the dual by a bivector. In fact, the function F(a) ~ a J A 2 
for any bivector A2 is skew, because the resulting bilinear product is antisymmetric: 

ao b = a\(b\ A 2 ) 
= (aAb)\ A 2 
= -(&Aa)J A 2 
= -b\{a\A 2 ) 

= -boa. (311) 

It turns out all skew functions are of this form. 

Theorem 49. F is skew iff F(a) = a\ A 2 for a unique bivector A 2 . 

Proof. I just fmished showing that any F of this form is skew. Knowing F, I can reconstruct 
A 2 uniquely using any frame {a,} and Theorem 1381 when r = 2: 

A 2 = i^a l A(a,J A 2 ) 

i 

= \j2 alAF ^- ( 312 ) 
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Now assume F is skew and let A2 be defined as above. I find that for any a, 
a\A 2 = 1 -a\ ^ fl 'AF(«,)J 

i i 
i i 

= \F (j2(a J a l ) a)j + ± £(J-(a) j a,) a* 
= lF(a) + \F{a) 

= F(a). (313) 

□ 

Therefore every antisymmetric bilinear product is of the form a Ab J A2 — a J A2 \_b for some Ai . 
9.5. Isometries and orthogonal transformations 

The final special linear operator is an isometry, which preserves inner products: F(u) J F(v) = u J v. (Equiv- 
alcntly, isometries preserve squares of vectors.) Isometries are always one-to-one, because 

F(u) = Q=> F{u) J F(v) = for all v 
=> u J v = for all v 

=^u = 0. (314) 

So in finite dimensions, isometries are also onto and thus invertible. An invertible isometry is called an 
orthogonal transformation. The two are distinct only on infinite-dimensional spaces, but most of the results 
Til show don't actually require invertibility, so I'll continue to make the distinction. 

Any power of an isometry is also an isometry, as is clear from the definition. An isometry satisfies 

u\ v = F(u)\ F{v) = FF(u)\ v, (315) 

so if F is an isometry then FF = ld. If F is also invertible, then its inverse has to be F, so we also have 
FF = ld. Therefore orthogonal transformations satisfy F = F^ 1 and are also normal, and as a bonus F^ 1 
is orthogonal too. And as with isometries, any power of an orthogonal transformation is also orthogonal. 

9.5.1. Isometries and versors 

So far I've described three isometries: the parity operation (which was extended to the whole algebra as grade 
involution in Section lÜTTj) , reflections, and rotations (both in Section[7]). Now a rotation is two reflections, 
and as I showed in Section 17.2.21 the parity operation is reflection in a volume element, which amounts 
to n reflections in succession. So every isometry I've shown so far is a composition of reflections. That's 
no accident: the Cartan-Dieudonné theorem shows that every isometry in an n-dimensional space is the 
composition of at most n reflections along axes. That's fantastic news, because reflections are easy to do in 
geometrie algebra; so now we have powerful tools to perform and analyze any isometry at all. 

So what does a general isometry look like? Remembering Eq. (I237[) , I find that the isometry F that takes 
vector u and reflects it along axes 01, 02, . . . , a r in succession is 

F(u) = (— l) r (a r ■ • • ü2 ai)u(a^ 1 a 2 " 1 • • • a,^ 1 ) 
= (a r ■ ■ ■ ü2 a\)u* r {a r • • • 02 ai) 1 
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= A r u* r A- x 



(316) 



where A r = a r • • • ai a\ is an invertible r-versor. Thus a general isometry in finite dimensions is grade 
involution followed by conjugation with an invertible versor. (This is why I defined versors in the first place, 
and it's also why I've been proving so many results not just for blades but for versors in general.) Now this 
looks a lot like Eq. (I239P for reflecting a vector in a subspace; in fact, Eq. (I239P is just a special case of this 
result, since a blade is a special type of versor. Therefore this operation extends to the whole algebra the 
same way reflection in subspaces did in Section 17.2.21 a general isometry on multivectors takes the form 

F{B) = A r B* r A; 1 (317) 

and it reduces to reflection in a subspace iff the versor A r is an r-blade. This also makes it clear that the 
isometries generated by A r and A~ x are inverses of each othcr. 

Even though every versor is associated with an isometry, the association isn't exactly one-to-one. After 
all, A r and \A r generate the same isometry for any A ^ 0. (The ultimate reason for this is that a and Xa 
represent the same axis, and thus the same reflection.) We can eliminate most of that ambiguity, however, 
by composing our versors out of unit vectors; in that case, A r is a unit versor. That doesn't eliminate the 
sign ambiguity, but we can live with that. Is there any further ambiguity? Amazingly, no. IT1 show this in 
two steps. First, basically the same argument used to derivc Eq. (|247l) shows that 

A r r r a~ x = (-iyi, (3i8) 

so isometries dividc into two classes: even isometries, which leave I alone and are represented by even 
versors, and odd isometries, which change the sign of I and are represented by odd versors. This also shows 
that an isometry is odd iff it's the composition of an even isometry and one reflection. As I said I would 
back in Section [21 I'll now start referring to any even invertible versor as a rotor, so rotors represent even 
isometries. 

Now I'll prove the result. 

Theorem 50. Versors A r and B s represent the same isometry iff A r = XB S for some X =/= 0. 

So if we consider only unit versors, the association of versors to isometries is exactly two-to-one. 

Proof. If A r is a nonzero multiple of B s , we know they represent the same isometry, so let's 
prove it the other way. Suppose A r u* r ^4~ 1 — B s u* s B~ 1 for all u. Since r and s are both even 
or both odd, I can drop the grade involutions and I'm left with A r uA~ l = B s uB~ l , which can 
be sneakily rewritten 

B^A r xu = for all u. (319) 

Then Theorem[25]tells me that B~ 1 A r = X+fj, (J)_. But B~ x A r is even, so its odd part vanishcs 
and I'm left with Bj 1 A r = X. Now if A = 0, by Theorem [24] both A r and B s would be null, 
which they arcn't since I've been inverting both of them. Therefore A ^ and A r = XB S . □ 

You might be tempted at this point to associate an isometry with a unique sequence of reflections, but you 
can't. That's because the factorization of a versor into vectors isn't unique. For example, suppose a and b 
are orthogonal Euclidean unit vectors; then 

fa - b\ fa + b\ , 

° i= br)br)- (320) 

so in this case two different sequences of reflections give the same isometry. In fact, reflections along any 
orthonormal basis for a subspace will result in reflection in that subspace, so in that case infinitely many 
reflection sequences produce the same isometry. But that's only because they all produce the same versor 
to within a sign. 
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9.5.2. Rotors and biversors 

Nothing I've done in this section up to now has made any assumptions about the scalars, but for this last 
part I assume the scalars are real. Every rotor is a product of biversors, so I want to take a moment to 
examine them. Consider ab where a and b are unit vectors; it represents the composition of reflections in 
the b and a directions in succession, so the resulting isometry acts in the a Ab plane. My plan to analyze ab 
is to expand itasaJö + aAfe and figure out each piece separately. I'U do that by starting with Eq. ((5|) from 
back in Section 11.11 

a 2 b 2 = (oj b) 2 ~ (aAb) 2 . (321) 

Since a and b are unit vectors, the left hand side of Eq. (|321[) is ±1. In what follows, I will set fi = a\b, and 
I will set aAb = XB where B is a 2-blade. In the cases when (a A b) 2 ^ 0, Hl choose A so B is a unit blade; 
othcrwisc Hl come up with some other way to choose A. 

Before I get into the general cases, Hl handle a special case that Hl need to refer back to later: let a 
and b be orthogonal. Since a\ b = 0, (et A b) 2 = ±1. So ab = aAb is already a unit blade B. To see what 
isometry it generates, let u lie in B; then its product with B is an inner product, so they anticommute, so 

BuB 1 = uBB 1 = — u. (322) 

So versor B generates a reflection in the plane it represents. (Which we already knew from Section f7. 2. II ) 
The result in the general case depends on the sign of B 2 . 

1. Suppose first that B 2 — — 1. Then Eq. p75]) whe n n = 2 tells me that the number of negative-square 
vectors in a frame for B is cithcr or 2, so the inner product on B is either positive definite or negative 
dcfinitc, which I call the Euclidean or elliptic case. Then a 2 and b 2 have the same sign, so a 2 b 2 = 1. 
Putting all this in Eq. (|32Tj) , I find 

l = /i 2 + A 2 . (323) 
Thcrefore /i = cos(0/2) and A = — sin(0/2) for some 9, so 

ab = a\b + aAb 

= cos(0/2) -Bsin(0/2) 

= exp(-B0/2) (324) 

where the exponential is defined by its power series. You may remember this from the end of Section 
11.21 it's a rotation through angle in the plane defined by B. When = it, I recover the special case 
I solved above: a rotation by 7r in a Euclidean plane equals a reflection in the plane. 

2. Now suppose B 2 = 1. In this case the inner product is indefinite, which is called the hyperbolle case. 
Now I have to give some thought to a 2 and b 2 . First let them have the same sign, so Eq. (|321[) becomes 

l = /i 2 -A 2 . (325) 

Thcrefore fi = ±cosh(</>/2) and A = =p sinh(</>/2) for some <f), so 

ab = a\b + aAb 

= ± (cosh(</»/2) - B sinh(</»/2)) 

= ±exp(-Bct>/2) (326) 

where again the exponential is defined by its power series. (This isometry, by the way, is a rotation in 
the hyperbolic plane, and in special relativity it's a boost to velocity ctanh</>.) This time I couldn't 
absorb the sign of a J b into a choice for the parameter, because cosh is always positive; but that affects 
only the rotor, not the corresponding isometry. So aside from that, the rotors for the last two cases 
have the same polar form, and the difference in their expansions as scalar plus bivector is due to the 
different behaviors of the area element. 
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3. Sticking with B 2 = 1, now I consider a 2 = —b 2 . Eq. (I321j) becomes 



-1 = u 2 — A 2 , (327) 

so jU and A change roles: u — =p sinh(</>/2) and A = ± cosh(</>/2) for some (/>, so for the rotor I get 

ab = a\b + aAb 

= =F sinh(0/2) ± B cosh(0/2) 
= ±B (cosh(0/2) - Bsinh(0/2)) 

= ±Bcxp{-B(j)/2) . (328) 

This rotor is the product of the previous one and the area element. As I showed in the special case 
above, this extra factor generates a reflection in the plane. Why is it showing up as a separate factor? 
Because unlike the Euclidean case, there is no hyperbolic rotation that performs a reflection in the 
plane, so it has to be included separately. 

Finally, suppose B 2 = 0, which means by Theorem [23] that the inner product is degenerate. (This 
also has a name: the parabolic case.) Now Eq. (|321[) reduces to ±1 = /i 2 . This doesn't make sense if 
the left hand side can be —1, so let me show that it can't. The inner product may be degenerate, but 
it can't be identically zero, or every vector would be null and there would be no axes in the plane to 
rcflcct along. Therefore there's a non-null vector somewhere in there which I'll call v, and the direction 
orthogonal to it is a null vector which together with v spans the plane. Because of this, the length 
squared of any vector in the plane is just v 2 times the square of its component along v, so they all have 
the same sign. Thus a 2 b 2 = 1 and u = ±1. Let B be any 2-blade that's convenient to use to represent 
the plane; then aAb — ^f\B for some A, so the rotor becomes 

ab = a J b + aAb 

= ±exp(-.BA/2) (329) 
where once again the exponential is defined by its power series. 

You may wonder what this rotor does. Its inverse is ±(1 + ^ B), so for any vector u 



A A 2 

= u + ~(uB - Bu) BuB 

2 4 

A 2 

= u + Xu\B BuB. (330) 

Each term on the right hand side is a vector, and you can directly verify that the square of the wholc 
thing really is u 2 . The verification is an interesting exercise; you find that u\ B is orthogonal to u, 
BuB is orthogonal to u\ B, BuB is null, and the inner product of u and BuB cancels the square of 
u\B. 

Putting all this together, I've shown that a general even isometry on a real vector space consists of any 
number of rotations in planes (elliptic, hyperbolic, or parabolic) and reflections in hyperbolic planes. An 
odd isometry is the same thing plus one reflection along an axis. All the rotations can be represented in the 
same polar form; the properties of the different area elements produce different types of rotations. If the 
whole space is Euclidean, things simplify further because there are no hyperbolic or parabolic planes: every 
isometry is a sequence of rotations in planes, preceded (or followed) by one reflection if it's odd. 

The set of all isometries on a fmite-dimensional real vector space forms a group called the orthogonal 
group on that space. All of this analysis tells us two main things about such groups: 
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1. The subset of even isometries forms a group of its own (it includes the identity and is closed undcr 
products), and the subset of odd isometries is a one-to-one copy of the even subgroup. A reflection 
along any axis provides the relation between the two subsets. 

2. Aside from reflections in hyperbolic planes, all elements of the even subgroup are functions of param- 
eters that can be continuously varied to zero, which results in the identity transformation. 

These properties teil me that a space's orthogonal group is an example of a Lie group. This is a group which 
can be divided into a finite number of isomorphic subsets, and within each subset the elements can be labeled 
with a finite number of continuously- variable parameters. Whichever subset is lucky enough to contain the 
identity is a subgroup in its own right, so a Lie group is a continuously-parameterised group together with 
some isomorphic copies. (In our case, reflections along axes and in hyperbolic planes move us back and forth 
between the copies.) After I learn how to use geometrie algebra to study Lie groups in general, I believe I'll 
be showing that they form even subalgebras which are generated by exponentiating bivectors. But this is 
more than enough on isometries for now. 

9.6. Extending linear functions to the whole algebra 

So far, I have considered only linear functions dcfincd on vectors; but functions on vectors have an obvious 
extension to the whole algebra. For example, consider the r-dimensional space spanned by {ai}; this is 
mapped by linear function F to the space spanned by {F(üi)}. Since blades represent subspaces, it seems 
very natural to define F not just on vectors but on blades too; I set 

F(ai A • ■ • A Or) := F(a x ) A • • • A F{a r ). (331) 

If I then require this formal extension of F to be linear over multivectors, I get F(A A B) = F(A) A F(B) for 
any A and B. Wcll, almost any A and B; my picture doesn't teil me what F should do to scalars. Assuming 
I figure that out, then I've found a way to naturally extend any linear function on vectors to a function on 
the whole geometrie algebra that is not only linear but also respects outer products. As a matter of fact, 
I've actually done this already no fewer than six times. Four of the extensions were orthogonal projections, 
orthogonal rejections, reflections, and rotations in Section[71 and the fifth was general isometries in Section 
19.51 I defined all five on vectors to start with, and I extended them in exactly the way I just suggested: I 
had them respect outer products. The sixth extension was grade involution in Section [5.1) I started with the 
parity operation u — > —u, but I extended it by making it respect not outer products but the full geometrie 
product. (It ended up respecting outer products too, as the third equality in Eqs. (|120j) shows.) In retrospect, 
that looks a little daring; after all, a product of vectors contains terms of many different grades, and the rules 
I imposed on, say, three-fold and fivc-fold products could have put conflicting requirements on trivectors. So 
far, though, it looks like everything worked out. Whew. While my other five extensions look safer, because 
I didn't have different grades crossing over, they do raisc a question: how do they act on products? Do 
you operate on each factor separately and then multiply them back together? That looks like it'd work for 
general isometries (and thus reflections and rotations) because the internal factors would cancel out, but I'm 
not too sure about projections and rejections. And then there's one property that all six extensions have 
in common: they leave scalars alone. What's up with that? Grade involution does it by definition, but the 
others were found to do so after their final definitions were stated. How come? 

9.6.1. Outermorphisms 

To answer these questions, I want to lay some groundwork by describing something a little more general. 
A linear function on geometrie algebras that preserves outer products is called an outermorphism. It's not 
too hard to show that the composition of outermorphisms is also an outermorphism and that the inverse of 
an outermorphism, if it exists, is an outermorphism too. However, a linear combination of outermorphisms 
is not an outermorphism. To show this, let F = aF\ + where F\ and Tv are outermorphisms, and 
try evaluating both F(AAB) and F(A) A F(B). You'11 see the problem pretty quickly. For this reason, 
compositions and inverses of outermorphisms pop up frequently, but linear combinations don't. 

Next, outermorphisms are very restricted in how they handle scalars. First linearity has a say: an 
outermorphism F has to be linear over everything, not just vectors, which means F{\) — F{\\) = XF(1). 
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So .F(l) determincs T for all scalars. But scalar multiplication is also an outer product, so 1 = 1 A 1. Thus 
the outermorphism property requires J"(l) = J r (l)AJ r (l). So like 1, .F(l) must equal its outcr product 
with itself. Something that equals its own square is said to be idempotent; since we have several different 
products, we have several different types of idempotency. I've pointed out that 1 is an outer idempotent, and 
thanks to the outermorphism property, .F(l) has to be an outer idempotent too. It turns out there aren't 
too many of those. 

Theorem 51. IfA = AAA, then A = or 1. 

Proof. Let A — J2 r Ar', then A = A A A becomes 

A = Y / (A s A t ) s+t . (332) 

s,t 

I'll look at this one grade at a time. The grade-r part of this expression is 

r 

A r = J2 (A s A r - s ) r . (333) 

s=0 

When r = 0, this becomes 

A = A 2 , (334) 

so A is either or 1. When r = 1, I find 

A 1 = 2A A 1 . (335) 

Whether A is or 1, this equation requires A\ = 0. 

To show that all remaining A r vanish, I proceed by induction. Suppose it's true for r — 1; 
then most of the terms in the sum for A r drop out, leaving 

A r = 2A A r . (336) 

Whether Aq is or 1, this gives me A r = 0, and that completes the proof. □ 

So J"(l) = or 1. That means outermorphisms can do only two things to scalars. 

Theorem 52. If T is an outermorphism, then either T = or J"(A) = A. 

This is why all six extensions left scalars alone; they had to. 

Proof. Either J r (l) = 0orJ r (l) = l. In the former case, for any A 

T(A) = F{IA) 
= J(1AA) 
= J r (l)AJ(A) 

= 0, (337) 

so T = 0. In the latter case, J"(A) = A for all A by linearity. □ 

Ncxt, PU take a passing look at adjoints. Linear functions on geometrie algebras have adjoints just as 
they do on any other vector spaces: if T : Q 1 — > Q 2 is linear, T : Q 2 — > Q 1 is a linear function given by 

T{B) * A := B * T(A) for all A e Q 1 and B € Q 2 . (338) 

Outermorphisms are linear, so they have adjoints which are linear. I can't say more than that, though, until 
I consider a special class of outermorphisms, which is what's next. 
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9.6.2. Outermorphism extensions 

Now, back to the reason we're here. My goal is to start with a linear function from vectors to vectors, dcfinc 
it on scalars by having it leave them alone (as Theorem [SD says I have to do), and then extend it to the rest 
of the algebra by linearity and respecting outer products. The result is an outermorphism that matches the 
original function on vectors. Can I start with any linear function and do this? Why sure; in fact, I can do 
it exactly one way. 

Theorem 53. Let F be a nonzero linear function that maps vectors to vectors; then there exists a unique 
outermorphism [F] that reproduces F when applied to vectors. 

Proof. Existence is obvious, because definition on scalars and vectors, plus the outermorphism 
property, plus linearity is enough to define [F] on any multivector. Uniqueness follows for the 
same reason. □ 

I had to specify F ^ because technically if F = then Theorem [SD allows two extensions: one is 
[F] = 0, and the other is [F] (A) = (A). Both extensions vanish on all vectors and all blades; they diffcr 
only in the scalar option from Theorem 1521 thev use. 

It's clear that [F] maps r-blades into r-blades, so 

([F](A)) r = [F]({A) r ). (339) 

So by requiring [F] to preserve grade 1, I find that it preserves all grades. It also follows that 

[F] {A)* = [F] (A*) 
[F] (A)t = [F] (At) 

[F] (A)t = [F] (A*). (340) 

If Gi is the vector space of G, then it's also pretty clear that 

[Id 5l ]=Id a . (341) 

Now let's see how outermorphism extensions behave under composition, inverses, and adjoints; I'U show 
compositions and inverses first. 

Theorem 54. If U, V, and W are the vector spaces of G 1 , G 2 , and G 3 respectively, and F : U —ï V and 
G : V — ¥ W are linear, then 

(a) [GF] = [G] [F]. 

(b) if F is mvertible, [F]^ 1 = [i^ 1 ]. 

Proof. For part (a), there's no problem with the action on scalars or vectors, so all I need to 
check is the outermorphism property. I will check it for the product of vectors, which covers all 
the higher-grade cases too. 

[G] [[F] (ui A • • • A u r )] = [G] [F( Ul )A • • ■ A F(u r )} 
= GF{ Ul )A---AGF(u r ) 

= [GF] (ui A • • • A u r ). (342) 

Part (b) follows from part (a): 

[F- 1 ] [F] = [F-'F] 
= [Kt/] 

= Id e i, (343) 
and similarly [F] [F^ 1 ] = Id e 2, so [F" 1 ] is the inverse of [F]. □ 
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As for adjoints, outermorphism extensions obey these relations, which (for the moment) use an amazing 
number of brackets. 



Theorem 55. 



A\{{F](B)} = {F\({[F] (A)}\B) 
{[F] (A)}[B=[F] (A[{[F](B)}) 



(344) 



Therefore [F] = [F] . 



Because this theorem is hard to read, I'm going to state it in words. The purpose of the adjoint is to 
let you move F from one side of the scalar product to the other, as long as you change F to F along the 
way. This theorem says you can also move [F] from the "high side" of the inner product to the "low side," 
as long as (a) you change [F] to [F] and (b) you then act on the whole thing with [F]. 

Take a look at two extreme cases. If the multivector on the low side is a scalar, the inner products 
become products, [F] leaves the scalar alone, and it factors out. That's why there's an extra [F] (or [F]) 
acting on the whole thing. On the other extreme, if both sides have the same grade then the inner products 
are scalars and the extra [F] or [F] drops out because it leaves scalars alone. With a little tweaking, that 
gets us back to the definition of the adjoint in Eq. (|338[) . I'll do that tweaking in the proof. 

Proof. I'll prove only the first relation; the second relation is the reverse of the first with a 
few substitutions. The result is true for general A and B if it's true for A r and B s . If r > s both 
sides vanish identically, so let r < s. If r = 0, then A r — A and both sides reduce to A [F] (B s ). 
If s = 0, then r = so we're back to the previous case. For the remaining cases, I consider blades 
A r and B s . Next I'll prove r = 1 and any s > 1. It's true for s = 1, so assume it's true for s — 1, 
let B s = öA.B s _i, and consider 

a\ [F] (B s ) = a\ [F] (&AB a _i) 

= a\ [F(b) A [F] (£? s _i)] 

_!)-F(fe)A[ a J [F](B, 



= [a\F(b)} [F] (B a _0-F(6)A[aJ [F] (B.^)] 
[F(a) j b] [F] (JJ._0 - F(b) A [F] (F(a) j B a _0 
i]- [F] [&A(F(a)j £?._!)] 
!-fcA(F(a)J£? s 
F\ [F(o)J(6AB._i)] 



F [(F(a)\b)B 8 



F\ [(F(a)\b)B s 



l)] 



= [F](F(a)]B e ). 

Now for general r < s. Fix s and assume the result is true for r 
and consider 



1; then let A r = A T 



(345) 



A a 



A r \ [F] (B a ) = (A r _!Aa)J [F] (B B ) 
= A r _ 1 \(a\ [F] (B s )) 
= A r _ 1 \ [F] (F(o)jB.) 
= [F] [[F] (A r _ 1 )j(F(a)J J B s )] 
= [F][([F] (A r _ 1 )AF(a))\B s ] 
= [F] ([F] (A r _!Aa)JB s ) 
= [F] ([F] (A r )JB B ). 



(346) 



And that takes care of all cases. 
Finally, I'll show that [F] = [F] . 



A*[F](B) = (Al\{\F](B)}) 
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= ({[F] (At)}\B) 
= ({[F] (A)} t J b) 
= [F] (A) * B. 

Thus [F] satisfies Eq. ([555]) with F = [F] , so [F] = [F] . 



(347) 
□ 



Given the uniqueness of the outermorphism extension and its good behavior under composition, inverses, 
and adjoints, I will now happily drop the [F] notation and let F refer either to the linear function on vectors 
or the resulting outermorphism. That certainly makes Theorem [551 easier to read: I'll take A\ F(B) = 
F (F(A) J B) any day. 

Recall that the restriction of F to the subspace represented by A r is denoted F\ r ■ This restriction and 
F(A r ) are related in an important way 

Theorem 56. F{A r ) = iff rank(i 7 A 7 .) < T, which is true iff FA r is not one-to-one. 

Proof. F(A r ) — iff F maps A r to a subspace of dimension smaller than r, which means 
rank(i ?1 J 4 I ,) < r. Since A r is finite-dimensional, by the rank-nullity theorem this is true iff 
mJl(FA r ) > 0, so F is not one-to-one. □ 

So A r e KerF iff Ker(FA 7 .) 7^ {0}. Note the two different meanings of F in these two statements. 

Now only one mystery remains unsolved: why did grade involution turn out to be an outermorphism 
even though I made it preserve geometrie products instead of outer products? Because the parity operation 
on which it's based has a special property. 

Theorem 57. The following conditions on F are equivalent. 

1. F is an isometry. 

2. F(AB) = F(A) F(B) for all A and B. 

3. F{A \ B) = F(A)\ F(B) and F(A[B) = F(A) [F(B) for all A and B. 

Proof. I first assume F is an isometry. That means F(a) J F(b) — a\ b for all vectors; but 
a\ bis ascalar, so F(a\ b) = a\ b, so F(a\ b) = F{a)\F{b). Since F(aAb) = F(a) A F(b) by the 
outermorphism property, it follows that F(ab) — F(a)F(b). I can extend this to F(a\a2 • • • a r ) — 
F{a\)F{a2) • ■ • F(a r ) by induction: the result is true for r = 2, so assume it's true for r — 1 and 
consider 



F(a 1 a 2 ■ ■ ■ a r ) = F[ai J (a 2 • • • a r )} + F[ai A(a 2 • • • a r )} 
— F > ( — t y ~"a-\ I a, . ■ ■ ■ d, ■ ■ • • a r 



J=2 



3=2 



+ F( ai )AF[a 2 ---a r 



■ F{a 3 ) ■ ■ ■ F{a r ) + F( ai ) A{F{a 2 ) ■ ■ ■ F(a r )} 



^(-l)-'- 2 F(a 1 ) J F(a 3 )F(a 2 ) ■ ■ ■ F{a 3 ) ■ ■ ■ F(a r ) + F(a x ) A{F(a 2 ) ■ ■ ■ F(a r )} 
F( ai )\ {F(a 2 ) ■ ■ ■ F(a r )} + F(oi) A{F(a 2 ) ■ ■ ■ F(a r )} 

F( ai )F(a 2 )---F(a r ). (348) 
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Now let's look at F(AB); by linearity it's a sum of terms of the form F(A r B s ). If r = or 
s = then F(A r B s ) = F(A r ) F(B S ) by linearity and Theorem I5"2l so let A r — a\a 2 ■ ■ ■ a r and 
B s = 6i&2' • • b s \ then Eq. (I348P Iets me show 

F(A r B s ) = F(a\a 2 ■ ■ ■ a r b\b 2 • • • b s ) 

= F( ai )F{a 2 ) ■ ■ ■ F(a r )F{bi)F(b 2 ) ■ ■ ■ F(b s ) 
= F(aia 2 ■ ■ ■ a r )F(bib 2 ■■■b s ) 

= F{A r )F(B s ). (349) 

Therefore F(AB) = F{A) F(B). Notice that the products among the equal outer products, 
and the same is true of products among the bj ; it is the product of a r and b\ that does not equal an 
outer product, and this is the reason I needed F(ab) — F(a)F(b), not just F(aAb) = F(a) A F(b), 
for the proof. 

Next, assume F(AB) = F(A) F(B) for all A and B. Then, since F commutes with all grade 
operators, 

F(A r \B s )=F((A r B s ) s _ r ) 
= (F(A r B s )) s _ r 
= (F(A r )F(B s )) s _ r 

= F(A r )\F(B s ), (350) 

with the last line following because F preserves grades. Thus 

F(A\B)=F(A)\F(B) (351) 

in general. Replacing s — r with r — s in the proof yields the samc rcsult for the right inner 
product. 

Finally, assume F{A\ B) = F(A)\ F(B) for all A and B; then F(a)J F(b) =F(a\b) = a\b 
for any two vectors a and b, where the last equality follows because a J b is a scalar. Therefore F 
is an isometry, and all three results listed above are equivalent. 

If I had assumed F was not just an isometry but orthogonal, I could have used F = F^ 1 
and Theorem [55] to prove the third part, but this way is better because it assumes less (on 
innnite-dimensional spaces at least). □ 

Thus an isometry may be extended by having it respect either outer products or geometrie products, 
with the same result. That won't work for anything else, though. This is why grade involution came out 
fine, and that's why I could indeed have extended reflections and rotations by respecting products instead 
of outer products. Projections and rejections had to be done the way I did them, however. 



9.7. Eigenblades and invariant subspaces 

Vector a is an eigenvector of F with eigenvalue A if F(a) — Xa, or equivalcntly a 6 Ker(_F — Aid). In this 
dcfinition a can't be zero, but A can. Now suppose a\ and a 2 are eigenvectors with eigenvalues Ai and A2 
and let A — a\ A a 2 ; then 

F(A) = F( ai Aa 2 ) 

= F( ai )AF{a 2 ) 
= AiA 2 ai A a 2 

= AxA 2 A. (352) 

So A is an eigenblade of F; it's mapped by F to a multiple of itself. 

An eigenblade is defined generally by F(A) — XA, regardless of what the factors of A do. Rotation 
operators nicely illustrate the different ways eigenblades can arise, and how they are related (or not) to 
eigenvectors. Let F(A) = exp(—B6/2)Aexp(B8/2), so F is a rotation in plane B through angle 9. Then 

F(B) = exp(-B6i/2)Bexp(BÖ/2) 
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= exp(-J36>/2) exp{B6/2)B 
= B, 



(353) 



so B is an eigenplane of the rotation operator with eigenvalue 1. However, since every vector in B gets 
rotated, in general none of them are eigenvectors. (The exception is B 2 = — 1 and 9 = tt, which is a 
refiection: every vector in the plane is an eigenvector with eigenvalue —1.) You can check that B 1 - is also 
an eigenblade with eigenvalue 1, but that's for a different reason: every vector in B 1 - is left alone by the 
rotation. So F(A) — A is consistent with F(a) — a for every vector in A, but it's consistent with many 
other things too. 

If A = 0, then Theorem [56l tells me that F is not one-to-one on A, but that's all it tells me. On the other 
hand, A =^ tells me all sorts of things. First, Theorem [56] says F is one-to-one on A, so Fa is invertible. 
But it also hints at what Range(FA) is. In fact, it's hard to see how F could map A to a multiple of itself 
unless it also mapped all members of A back into A. That would make A an invariant subspace of F; that 
is, a subspace that is mapped to itself by F. Put that together with Theorem [561 and you get this result. 

Theorem 58. Any eigenblade of F with nonzero eigenvalue is an invariant subspace of F on which F is 
invertible. 

Proof. Suppose F(A) — XA and A 7^ 0. We already know that F is invertible on A, so Fll 
prove the first part. If a lies in A, then A A a = 0, in which case 

AAF(a) = \~ 1 \A A F(a) 

= \- 1 F(A) / \F(a) 
= X~ 1 F(AAa) 

= 0, (354) 

so F(a) lies in A too. □ 

So F maps A invertibly onto A. Now this is true for any A 7^ 0; but what does the actual value of A teil us? 
Well, if F(A) = XA then |F(A)| 2 = A 2 |A| 2 , so the value of A determines how much the norm squared of A 
changcs. If the scalars are real, I can interpret this further. Recalling Eq. (|173[) for the weight of a blade, I 
find that 

weight(F(A)) = |A| weight(A), (355) 

so F multiplies the weight associated with A by |A|; and F changes the orientation of A if A is negative. 
This suggests to me that A is actually the determinant of Fa, since it seems to be the factor by which the 
volume of A changes. The idea that the determinant of a linear transformation is actually an eigenvalue, not 
of the original transformation but of its outermorphism extension, is worth following up on and has general 
validity, even if the scalars aren't real. So I'll do that next. 



9.8. The determinant 

The determinant of a linear transformation is the factor by which it multiplies the volume element of the 
space it acts on. Finding it in geometrie algebra is easy; we consider F(I). This is an n-blade, so it has to 
be a multiple of I. (Put another way, I is an eigenblade of all linear transformations.) That multiple is the 
determinant, or 

F{I) =: det(F)I. (356) 

An equivalent defmition is 

det(F) = F{I)^. (357) 

This way of dcfining the determinant is very intuitive and also very easy to use in calculations, as I'll show. 
First, it's obvious that det(Id) = 1. Now let F and G be linear transformations; since 

dct(FG)I = FG(I) 

= F(det(G)I) 
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I find with minimum fuss that 
Therefore if F is invertible, 



so 



= det(G)F(I) 
= det(G) det(F)J, 

det(FG) = det(F) det(G). 

det(F _1 )det(F) = det(F~ x F) 
= det(Id) 
= 1 

det^- 1 ) =det(F) _1 . 



(358) 
(359) 



(360) 
(361) 



That tells me that det(F) ^ if F is invertible. I'll use that later. 
Now for adjoints. From the definition in Eq. ()338[) . 

F(I) *I = 1* F{I) 
det(F)I * I = det(F)I * I 
det (F) = det (F). 

It's easy when you know how. 

Next I'll calculate the determinants of some specific operators. 

• From Eq. (1235)) I get the determinants of orthogonal projections and rejections: 

det (Pa,) = S rn 
det (Fa,) = S r o- 



(362) 



(363) 



Let F be a symmetrie operator with a frame {ai} of eigenvectors with eigenvalues {A^}. Then 
ai A • • • A a„ is a volume element, so 



det (F) ai A • • • A a n = F(a 1 A • • • A a„) 

= F{a 1 )h---hF{a n ) 
= Ai • • • A n ai A • • • A 



so 



(364) 
(365) 



det (F) = Ai ■••An. 

• If F is orthogonal, then F(v) = ArV^A^ 1 for some rotor A r , so Eq. (|318l) shows that 

dct(F) = (-l) r . (366) 

In real matrix algebra, there's a well-known relationship between determinants, invertibility, and adjoints: 
a matrix is invertible iff its determinant is nonzero, in which case its inverse is its adjugate (adjoint of the 
cofactor matrix) divided by its determinant. The geometrie algebra equivalent, as you would expect, applies 
to the linear operator itself, not its matrix representation on some basis. On top of that, it's valid for any 
multivector. 

Theorem 59. F is invertible iff det (F) ^ 0, and for any multivector A 

F(^- ± ) ± 



F-\A) 



det (F) 



(367) 
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Recall that is the duality transform and is its inverse. 

Proof. I've already proven that if F is invertible, det(F) 7^ 0, so let's go the other way. 
Suppose det(F) ^ and let G be defined by G(A) = T(A^ ± )- L . Then I use the result A\ F{B) = 
F(F(A)\ B) from Theorem [55] to get 

GF(A) = G[F(A)] 

= F[F(,4)- L ]- L 
= F[F(A)J Z]/- 1 
= [AjF(J)]/- 1 
= det(F)(AJ I)/- 1 
= det(F)A//" 1 

= det(F)A (368) 
A similar argument shows FG(A) = det(F)A, so F -1 = G/det(F). □ 

You can see that you need the full apparatus of geometrie algebra to do this: I take the dual, which is a 
geometrie product with a volume element, and I need to use the outermorphism extension of F, because 
calculating F~ 1 (a) involves calculating F(a^ ± ), and a~ ± ' is not a vector. (Unless we're in a two-dimensional 
space, I suppose.) 



10. Applications 

10.1. Classical partiele mechanics 

The most obvious place to apply geometrie algebra is classical mechanics, since it relies heavily on vector 
algebra already. In this section only I'U adopt the notational conventions of classical mechanics, so vectors are 
denoted a, 6, and so on, the magnitude of vector a is denoted a, unit vectors are indicated with an overhat 
~, and the derivative of any quantity with respect to time is indicated by an overdot ' . Since boldface means 
something else, I will not use boldface for blades in this section. The material in this section is largely drawn 
from [3] and g]. 



10.1.1. Angular momentum as a bivector 

As a partiele moves over time, its position vector r sweeps out area at a rate that a picture easily shows to 
be 

A=-rAv (369) 

where v — r is the particle's velocity vector. Unsurprisingly, the rate at which area is swept out is a 2-blade. 
This blade is proportional to the dynamical quantity 

L := r Ap = mr Av = 2mA, (370) 

callcd the angular momentum. In Standard vector algebra, angular momentum is defined to be the vector 
L = r x p, the cross product of r and p; the definition given here is the dual of that vector (see Eq. (12081) ). 
which is more natural given the association with areas. Nonetheless, since the algebraic properties of the 
outer and cross products are so similar, much of what one knows from the Standard treatment holds without 
change; for example, 

L = mv A v + mr A v 
= mr A v 

^ r AF, (371) 
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so L is conserved iff the force F is central (parallel or antiparallel to r). Since L is conserved ifï A — 0, I 
have Kepler's Second Law: the position vector of a partiele subject to central forces sweeps out equal areas 
in equal times. Further, the plane in which central force motion takes place is L itself. Writing r = rf, 
which implies v = ff + rf, I find 

I = rAp = mr A v 
= mrf A(rf + rf) 

= mr 2 fl\f (372) 

since f A f = 0. But I know a bit more than that; f is a unit vector, or f J f = 1, the time derivative of which 
is f J f = 0. This is just the familiar fact that a constant-length vector and its time derivative must always 
be perpendicular. (Incidentally, this shows that v = ff + rf is a decomposition into radial and tangential 
components.) If f J f = 0, then f A f = rf , so 

L = mr 2 ff . (373) 

Now this is nice because the geometrie product has better properties than the outer product, and this is the 
hrst algebraic feature of this treatment that is genuinely new. Since L is a bivector, 

L = -Ü = -mr 2 ff (374) 

so the scalar l, the magnitude of L, is given by 

l 2 :=\L\ 2 = iïL = -L 2 = m 2 r i f 2 . (375) 

Notice that l equals the magnitude of the angular momentum vector from Standard treatments. It should, 
of course, since bivector L and the angular momentum vector are duals. 



10.1.2. The Kepler problem 

The Kepler problem is to determine the motion of a point partiele of mass m moving in a potential of the 
form V = —k/r, where r is the particle's distance from some fixed origin. The partiele experiences a force 

F = - Af (376) 

where the constant k is positive for an attractive force and negative for a repulsive force, so the particle's 
acceleration is given by 

v = -=■#. (377) 

mr 

Now take a look at Eqs. (|373[) and (|377[) . One is proportional to r 2 while the other is inversely proportional 
to r 2 , so their product is independent of r. In fact, let me calculate the product: 

Li) = ( —mr 2 ff] ( —^f 



= kf, 

and since L is conserved and k is a constant, this implies 

4- (Lv - kf) = 0. (378) 
at 

Well, look at that: another constant of motion. The second term in the constant, kf, is clearly a vector, and 
the hrst term can bc writtcn 



Lv = L [v + L f\v 
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= L [ v + mr At) Ad 

= L[v. (379) 

This is hardly a surprise; v is a vector in the plane defined by L, so by Theorems 1161 and 1171 Lv = L\_v is 
a nonzero vector in the plane of L perpendicular to v. Thus the conserved quantity is a vector in the plane 
of motion; it is often called the "Laplace-Runge-Lenz vector," and in traditional vector algebra treatments 
of the Kepler problem it typically appears at the end as the result of a great deal of work. Here it was the 
first thing I found. 

I would actually prefer to define a dimensionless conserved vector, and this quantity clearly has dimensions 
of k, so I define a conserved vector e by 

Lv 

e:= — -f. 380 

k 

Fd like to use this equation to get further expressions describing the motion of the partiele; first is the polar 
equation, r as a function of direction. Since the expression for e has r in it and rr = r, it follows that I can 
get an equation for r by multiplying Eq. (|380j) by kr, with the result 



The left hand side equals 



m 

while rr = r and er = er cos 9 + e Ar, so putting it all together 



Lvr = k(rr + er). (381) 



Lvr = L(v \r + v Ar) 

T2 

= (r\v)L 

m 

l 2 

+ {r\v)L (382) 



ï 1 

h (r J v)L = k(r + ercosO + e Ar), (383) 



or on separating the scalar and bivector parts, 

l 2 

— = kir + er cos 9) 
m 

(r\v)L = eAr. (384) 
The scalar equation can be solved for r with the result 

l 2 /mk 



1 + e cos 9 ' 



(385) 



which is the equation for a conic section with eccentricity e and one focus at the origin. Since the length of 
e is the eccentricity of the orbit, e is naturally called the eccentricity vector, which is the name I'll use for 
it henceforth. 

The direction of e also has a geometrical meaning, but it's different in the attractive and repulsive cases, 
so I'll do one at a time. First I assume k > and I note that r equals its minimum and maximum values 
when 9 — and tt respectively, which means that e points toward the particle's point of closest approach, 
called its periapsis, and away from its point of farthest retreat, called the apoapsis. (Fun fact: these two 
points are called the perigee and apogee if you're orbiting the earth, the perihelion and aphelion if you're 
orbiting the sun, and the pericynthion and apocynthion if you're orbiting the moon. So now you know.) 

Now the repulsive case. If k < 0, we run into a problem: r has to be non-negative, so we have to have 
1 + ecosö < for at least some values of 9, and the orbit may include only those values. This is possible iff 
e > 1, with the result that the orbit is a hyperbola. In this case, r takes on its smallest value when 9 = tt, 
so in the repulsive case the eccentricity vector points away from the periapsis. 



80 



The motion in the Kepler problem is completely determined by two vectors, the initial position and 
velocity, which are themselves determined by six parameters. The conserved angular momentum supplies 
three parameters because it's a bivector, and the conserved eccentricity vector supplies two more (only two 
because the angular momentum fixes the plane of the motion), so the motion is completely determined by 
these two constants plus one further parameter, which may be taken to be the initial value of 6. If that's the 
case, then anything that doesn't depend on the starting point, such as any other constants of motion, should 
be a function of only L and e. Hl now show that this is the case for the energy by hnding the magnitude of 
the eccentricity vector. 

(Lv - kr) 2 = k 2 e 2 
(Lv) 2 - 2k(Lv) J f + k 2 = k 2 e 2 (386) 

Using the fact that a vector equals its own reverse, the first term on the lcft hand side can be calculated as 

(Lv) 2 = LvLv 

= (Lv)^Lv 
= vl?Lv 

= l 2 v 2 . (387) 

(Lv)\ r 
r 

(Lvr) 
r 

_ P_ 

mr ' 

where in the last line I used Eq. (|382|) . so now Eq. (|386j) becomes 

9kl 2 



The second term on the left is — 2k timcs 

(Lv)\ f 



(388) 



^ Q™ 2 - *) = *V - 1) (389) 

or 

S~(e 3 -1). (390) 

This gives the energy in terms of l and e. 

I have derived all the main results of the Kepler problem (except for the time evolution) a whole lot 
more easily than Standard treatments do. In fact, many textbooks don't even get to the eccentricity vector. 
Here geometrie algebra is clearly superior to Standard vector algebra both for solving the equations and for 
understanding the results. 



A. Summary of definitions and formulas 

A.l. Notation 

Q Geometrie algebra 

Gr Grade-r subspace of Q (space of r-vectors) 

Q n Geometrie algebra of an n-dimensional vector space 



A, B, etc. General multivector 

A, /j,, etc. Scalar 
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a, b, u, v, etc. Vector 

A r r-vector (sometimes grade-r part of A) 

A r r-blade 

A + Even-grade multivector 

A_ Odd-grade multivector 

I Volume element 

(A) r Grade-r part of A 

(A) Scalar (grade-0) part of A 

(A) + Even-grade part of A 

(A) _ Odd-grade part of A 

A^ 1 Inverse of A 

A* Grade involution of A 

A* r r timcs grade involuted A 

A^ Re ver se of A 

A$ Clifford conjugate of A 

\A\ 2 Squared norm of A 

A 1 - Dual of A 

A^ 1 - Inverse dual of A 

AB Geometrie product of A and B 

A\B Left inner product of A into B 

A [ B Right inner product of A by B 

AAB Outer product of A and B 

A* B Scalar product of A and B 

A x B Commutator of A and B 

PA r (B) Orthogonal projection of B into A r 

RA r (B) Orthogonal rejection of B from A r 

U, V, W, etc. Vector space 

F, G, etc. Linear function of vectors (or its outermorphism extension) 

F(j Rcstriction of F to subspace U 

Fa Restriction of F to subspace represented by A 

ld Idcntity function 

Ker(F) Kernei of F 

Range(F) Range of F 

null(F) Nullity of F 

rank(F) Rank of F 

F Adjoint of F 

det(F) Determinant of F 

T Outermorphism 



A.2. Axioms 

A geometrie algebra Q is a set with two composition laws, addition and multiplication, that satisfy these 
axioms. 

Axiom 1. Q is a ring with unit. The additive identity is called and the multiplicative identity is called 1. 
Axiom 2. Q contains a field Q of characteristic zero which includes and 1. 

Axiom 3. Q contains a subset Q\ closed under addition, and X e Go, v € Qi implies Xv = vX £ Q\. 
Axiom 4. The square of every vector is a scalar. 
Axiom 5. The inner product is nondegenerate. 
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Axiom 6. Iföo = Qi, then Q = Q . Otherwise, Q is the direct sum of all the Q r . 



A.3. Contents of a geometrie algebra 

An r-blade A r is the outer product of r vectors, ai A • • ■ A a r . It represents the subspace spanned by 
{a,j}j=i,..., r , with a weight and orientation if the scalars are real. 

A r = iff the dj are linearly dependent. 

A r and B r dcfine the same subspace iff A r = XB r . 

If A r is a proper subspace of A s , then A r can be factored out of A s from either the left (A s = A r A A s _ r ) 
or the right (A s = A s _ r A A r ). The grade-s — r factors in each case may be chosen to be the same except 
for at most a sign. 

a A A r — iff a lies in A r . 

a J A r = iff a is orthogonal to A r . 

The rcfiection of multivector B in subspace A r is A r B* r A~ x . 

A r , A*, Al, A\,, and A~ x (if it exists) represent the same subspace. 

A.4. The inner, outer, and geometrie products 



min{r,s} 




(1) 



j=0 



(A r B s ) 



r+s-2j 



. = (-1)""^ (B s A r ) 



r+s-2j 



(2) 



(AB) = (BA) 



(3) 
(4) 
(5) 
(6) 



= (A*B*) 
= (A^Bi) 
= (A*B*) 




(7) 




(8) 




(9) 



A r \ B s 
A r A B s 



(- 1 y(s-i) Bs[Ar 
(—l) rs B s AA r 



(10) 
(11) 



d\ A a 2 A • • • A a r = (a\ü2 • • • a r ) r 



(12) 
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a\A=^(aA-A*a) (13) 
aAA= -(aA + A*a) (14) 



A[a = -a\A* (15) 
AAa= aAA* (16) 

a A A Ab = —bAAAa (17) 



a\(AB) = (a\A)B + A*(a\B) (18) 

= (a A A)B — A* (a A B) (19) 

aA(AB) = (aAA)B- A*(a\ B) (20) 

= (aj A)B + A*(aAB) (21) 

a J (A A £>) = (a\ A) A B + A* A(a J £>) (22) 

aA(ALB) = (aAA) [B - A* [{a\ B) (23) 

a A(A \ B) = (a \ A) \ B + A* \ (a A B) (24) 



aj (ai A a 2 A • • • A a r ) = ^(-l) : ' _1 (a j aj) a\ A a 2 A • • • A &j A ■ ■ ■ A a r (25) 
i=i 

ai A(a 2 A • • • A a r ) = a\ A a 2 A • • • A a r (26) 

If r < s then 

B r J (ai A a 2 A • • • A a s ) = ^2(-l)^=^~ j) (B r J a h A a l2 A ■ ■ ■ A a ir ) a lr+1 A ■ ■ ■ A a ls (27) 

where the sum is performed over all possible choices of {a» }j=i...., r out of {aj}i=i,..., s , and in each term i\ 
through i r and i r+ i through i s separately are in ascending order. 

AA(BAC) = (AAB)AC (28) 

A\(B[C) = (A\B)[C (29) 

A \ (B \ C) = (A A B) \ C (30) 

A[(B AC) = (A[B)[C (31) 



If A = aia 2 ■ ■ ■ a r , then 
and 



AB.rf = {AB s A^) s (32) 
(ABAi)A(ACAt) = \A\ 2 A{B A C)A^ (33) 
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A.5. The geometrie meaning of the inner and outer products 

A r A B s = iff A r and B s share nonzero vectors. 

A r AB S , if nonzcro, represents the direct sum of A r and B s . 

A r J B s = iff A r contains a nonzero vector orthogonal to B s . 

A r \ B s , if nonzero, represents the orthogonal complement of A r in B s . 

If A r and B s are orthogonal, then A r B s = A r A B s . 

If A r is a subspace of B s , then = A r J 

The converses of the previous two statements are true if (1) r = 1 or s = 1 or (2) A r or B s is invertible. 

P A ^(B)=BjA r A T : 1 = (BjA r )JA- 1 (34) 
R Ar (B)=BAA r A; 1 =BAA r [A- 1 (35) 



A.6. Grade involution 



A* = A (36) 

a* = -a (37) 

(AB)* = A*B* (38) 

(A + B)*=A*+B* (39) 

= (-l)M r (40) 

A* = <A) + -<A)_ (41) 

= (42) 

(A) ± = ^±^*) (43) 

A** = A (44) 

(A- 1 )* = (A*)- 1 (45) 

(AJB)*=A*JB* (46) 

(A[B)* = A* [B* (47) 

(AAB)* = A* AB* (48) 

AI = IA< n -V (49) 
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A.7. Reversion 



A.8. Clifford conjugation 



At = A (50) 

a 1 = a (51) 

(AB)t = B^t (52) 

(4 + B) f = A f + ^ (53) 

4 = (-ïy^-WAr (54) 

v4 tt = A (55) 

(A-y = (At)- 1 (56) 

(AjB)t = flt [At (57) 

= BtjAt (58) 

(AAB)t = BtAit (59) 



A* = A (60) 

at = -a (61) 

{ABf = (62) 

(A + B)* = A* + (63) 

At = A*t = At* (64) 

At = (-l) r ('' +1 >/ 2 A r (65) 

Att = A (66) 

(A^)t = (At)- 1 (67) 

(A\B) i = B t [A i (68) 

{A[B) t = B i \A i (69) 

(AAB)* = B* AA* (70) 
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A.9. The scalar product and norm 



A*B = (A*B) (71) 

= (A*\B) (72) 

= (Ai[B) (73) 

A * B = ^ A r * B r (74) 

r 

= J2 A t\Br (75) 

r 

= J24lB r (76) 

r 



A * B = B * A 

= A* *B* 

= A t *B t 



(77) 
(78) 

(79) 
(80) 



A * (BC) = (B^A) * C (81) 

A*(B[C) = (B^[A)*C (82) 

A* (B\ O) = (B f AA) *C (83) 

A*(BAC) = (B t J,4)*C (84) 



Multivector A is uniquely determined by either of the following: 

1. A * B for every multivector B. 

2. (A) and aj A for every vector a. 



\A\ 2 = A*A (85) 
= \A*\ 2 (86) 
= l^l 2 (87) 
= \At\ 2 (88) 

If A = a\a,2 ■ ■ ■ a r , then 

(a) \A\ 2 = tfA = a\a\---a 2 r . 

(b) A- 1 cxists iff |A| 2 ^ 0, in which case A^ 1 = Aï/\A\ 2 and |^4. 1 j 2 = \A\- 2 . 

(c) (AB) * (AC) = (BA) * (CA) = \A\ 2 B * C. 
For any bladc A r , 

(a) \A r \ 2 — iff the inner product is dcgenerate on A r . 

(b) A- 1 cxists iff \A r \ 2 ^ 0, in which case A~ x = (-l)^" 1 )/ 2 ^/^! 2 = A r / A 2 . 
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A.10. The dual 

A volume element I is a unit n-bladc. 



A ± =A\I~ 1 (89) 
= AI' 1 (90) 



A~ ± = A\I (91) 
= AI (92) 
= I 2 A ± (93) 



A 1 - is the orthogonal complement of A. 



A.ll. The commutator 



(AB) 1 - = AB 1 - (94) 

(AAB)^ =A\B ± (95) 

(A\B) ± =AAB ± (96) 

(A ± )~ 1 = IA- 1 (97) 

(A*) 1 - = (-1)™(^)* (98) 

(At)^=[/t(^)t] ± (99) 

{A i ) ± = [I i (A ± ) i ] ± (100) 

* B 1 - = \I\~ 2 A* B (101) 

A^B< n ~ r \A^)- 1 = (A r B* r A^ 1 )* (102) 

ik(a)=iV(a) (103) 

A x B = -(AB -BA) (104) 

A x (BC) = (Ax B)C + B(A x C) (105) 

A x (B x C) + B x (C x A) + C x (A x B) = 0. (106) 

(A x B)* = A* x B* (107) 

(A x B) t =B t xA t (108) 

(A x Bf = B t xA t (109) 
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A x A = O 



(110) 



a x A = a\ (A) + + aA (A) _ (111) 
Ax a= (A) + [a+(A)_Aa (112) 

A 2 x A r = (A 2 A r ) r (113) 

A 2 x (B\ C) = (A 2 x B)\ C + B\{A 2 x C) (114) 
A 2 x(B[C) = {A 2 xB)[C + B[(A 2 xC) (115) 
A 2 x (BAC) = (A 2 x B)AC + BA{A 2 x C) (116) 

r 

A 2 x (ai A a 2 A ■ ■ ■ A a r ) = J^ai A a 2 A ■ ■ ■ A(A 2 [ aj) A ■ ■ ■ A a r (117) 

3=1 

A commutes with all multivectors iff A commutes with all vectors iff A = A + fi{I)_. 
A.12. Frames and bases 

If {Oi}i=l,.. 

jn is a frame with volume element ün — ai A • • • A a n , the reciprocal frame is given by 

a l = (-l) ï_1 (ai A a 2 A • • • A a t A ■ ■ ■ A a n )a^ 1 . (118) 

It satisfics 

o*Joj=(55. (119) 
Let I be an increasing string of indices i\, i 2 , . . . , i r ; then o/ and a 1 are 

aj — a il A a i2 A ■ ■ ■ A cn r (120) 
a^a^A^A'-Ao^. (121) 

Thcy satisfy 

a/ * a J = 5j, (122) 

and for any multivector A, 

A = ^2 A ' a i wncrc A 1 = A* a 1 (123) 
i 

= ^2 Aki 1 where A 7 =^*a/. (124) 

If I is increasing and I c is the increasing string of indices complementary to /, then 

a 1 = (_l)i:;=i(*i-i) 0je o-i. (125) 
A frame and its reciprocal satisfy these identities: 

a l dj J A r = a l A(aj J A r ) = rA r for any A r . (126) 

i i 

^2 a { a % J A r = ^ a i A ( fli J A r) = rA r for any A r . (127) 
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cii a l = a 1 ai = ra. 

i i 

The volume element of the reciprocal frame, a N = a 1 A • • • A a n , is also given by 

N a N 
ü = i 12- 

\aN\ 2 

I 

= ^(Bïa 1 )^ 
I 

A.13. The adjoint of a linear operator 

F(B) *A = B* F(A) 

A\F(B)=F{F{A)\B) 
F(A)[B = F(A[F(B)) (133) 

T = F (134) 
GF = FG (135) 
F' 1 =F Z=T (136) 

det(F) = dct(F) (137) 

A.14. Symmetrie and skew symmetrie operators 

F is symmetrie if F = F and skew symmetrie (or skew) if F = —F. 

F is (skew) symmetrie iff a\ F(b) is (anti)symmetric. 

F is symmetrie iff F(a) — J2i ^i( a * a *) a i f° r some frame {ai} of eigenvectors with eigenvalues {Ai}. 
F is skew iff F(a) = a\ A 2 for some bivector A 2 . 

A.15. Isometries and orthogonal transformations 

F is an isometry if F(u) J F(v) — u J v for all u and v. 

F is an isometry iff FF = ld. 

F is an orthogonal transformation if F is an invertible isometry. 
F is orthogonal iff F = F -1 . 

F is an isometry on a finite-dimensional space iff F(a) — A r a* r A~ 1 for some invertible r-versor A r . 
The extension of orthogonal F on a finite-dimensional space to all multivectors is F(B) = A r B* r A~ 1 . 
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(128) 

(129) 

(130) 
(131) 

(132) 



A.16. Eigenvalues, invariant subspaces, and determinants 

A ^ O is an eigenblade of F if F(A) = \A. 



An eigenblade of F is an invariant subspace on which F is invertible. The eigenvalue is det (Fa) ■ 



dctF = F(I) 1 - 



(138) 



F(Ar^) 
det(F) 



_L 



F~\A) = 



(139) 



B. Topics for future versions 

These are the subjects I plan to add to the notes next, in no particular order. The items with asterisks are 
most interesting to me at the moment. 

• Linear algebra 

— More on invariant subspaces and determinants* 

— Representing a general linear operator as a sequence of multivector multiplications* 

• Differcntial and integral calculus 

— The directed integral of a multivector 

— The derivative of a multivector- valued function defined in terms of the directed integral* 

— Recovering traditional vector calculus 

— The fundamental theorem of calculus and its corollaries (Gauss' theorem, Stokes' theorem, Green's 
theorem, etc.)* 

— Taylor series* 

— Generalizations of Cauchy's integral formula* 

— The invertibility of the derivative (cf. the exterior derivative) 

— Solutions to Standard ODEs and PDEs (simple harmonie oscillator, wave equation, etc.)* 

— Fourier analysis 

— Manifold theory 

— Lie groups and Lie algebras 

— Curvature* 

• Geometry 

— Meet and join of subspaces 

— Projective splits (e.g. Minkowski spacetime into any observer's space + time)* 

— Different models of space (Euclidean, projective, conformal) 

— Geometrie algebra on a vector space without an inner product 

• Physics 

— Rotational dynamics and the inertia tensor* 

— Relativistic partiele mechanics* 

— Electricity and magnetism in 3D and 4D* 

— Lagrangian and Hamiltonian mechanics* 

— Continuüm mechanics and elasticity theory* 

— The Dirac equation 

— General relativity 

— The Galilei and Lorentz groups and their Lie algebras 
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